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PREFACE  TO  THE   THIRD  EDITION. 


AT  the  request  of  the  Eev.  Charles  Cheyne,  the  father  of  the 
late  C.  H.  H.  Cheyne  the  author  of  this  work,  who  to  the 
great  regret  of  all  his  relatives  and  friends  died  at  Torquay 
on  the  1st  of  January,  1877,  I  have  undertaken  the  charge 
of  this  new  edition.  I  have  made  but  two  changes  in  the 
body  of  the  work,  and  have  introduced  but  few  notes.  Inas- 
much as  the  methods  of  Chapter  II.  are  peculiarly  the 
author's  own  I  have  not  ventured  to  alter  them  :  but  I  have 
instead  added  to  the  Appendix  the  Corollary  at  the  end  of 
Art.  5,  a  method  of  determining  the  variation  of  the  elements 
(Art.  6)  which  is  for  the  most  part  due  to  Pontecoulant  with 
a  slight  addition  of  my  own  to  make  it  complete,  together 
with  Art.  8  also  largely  due  to  the  same  author.  I  have 
also  given  some  additional  examples  and  have  indicated  the 
sources  from  which  they  were  obtained. 

If  a  new  treatise  on  the  Planetary  Theory  had  to  be 
written  it  would  be  necessary  to  digest  the  work  of  Jacobi, 
Hansen,  Leverrier  and  Newcombe,  to  mention  only  a  few 
eminent  analysts.  A  brief  account  of  a  variety  of  methods 
might  be  found  in  Matthieu's  Dynamique. 
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A  biographical  notice  of  the  late  author  may  be  read  in 
the  Monthly  Notices  of  the  Royal  Astronomical  Society, 
Vol.  xxxvu.  pp.  147,  148.  This  I  find  as  complete  as  is 
required  in  such  a  sketch  of  the  too  brief  life  of  a  diligent 
lover  of  Mathematical  and  Astronomical  Science. 

A.  FREEMAN. 


MURSTON  KECTORY,  KENT, 
October  31st,  1883. 
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PREFACE  TO  THE  FIEST  EDITION. 


IN  this  volume,  an  attempt  has  been  made  to  produce  a 
Treatise  on  the  Planetary  Theory,  which,  being  elementary 
in  character,  should  be  so  far  complete,  as  to  contain  all  that 
is  usually  required  by  students  in  this  University.  But  it 
is  not  without  diffidence  that  I  submit  my  volume  to  their 
notice.  In  the  earlier  part  of  it,  the  methods  which  have 
been  adopted  are  to  some  extent  original*,  and  the  general 
arrangement  of  the  second  Chapter  will,  it  is  believed,  be 
found  to  be  new.  Through  the  kindness  of  the  Publishers, 
a  portion  of  Pratt's  Mechanical  Philosophy  has  been  placed 
at  my  disposal.  Of  this  I  have  availed  myself,  particularly 
in  the  Chapter  on  the  Stability  of  the  Planetary  System ;  but, 
on  the  whole,  comparatively  little  has  been  reprinted  ver- 
batim from  that  work.  Among  other  sources  of  information, 
my  obligations  are  mainly  due  to  Pohtecoulant's  Theorie 
Analytique  du  Systeme  du  Monde,  Airy's  Mathematical 
Tracts,  and  Frost's  Planetary  Theory  in  the  Quarterly 
Journal  of  Mathematics:  but  I  have  also  referred  to  the 

*  Some  of  these  have  already  appeared  in  Mathematical  Journals. 
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Mfoanique  Celeste,  the  Mecanique  Analytique,  Mrs  Somer- 
ville's  Mechanism  of  the  Heavens,  (a  work  forming  a  complete 
Mathematical  Treatise  on  Physical  Astronomy,)  a  Memoir 
by  Prof.  Donkin  on  the  Differential  Equations  of  Dynamics, 
Phil.  Trans.  1855,  &c.  A  collection  of  Problems  has  been 
added,  taken  chiefly  from  the  Smith's  Prize  and  Senate- 
House  Examination  Papers  of  the  last  twenty  years.  In 
conclusion,  I  would  express  my  sincere  thanks  to  Messrs. 
A.  Freeman,  P.  T.  Main,  and  other  friends,  of  St  John's 
College,  for  the  valuable  assistance  which  they  have  afforded 
me,  and  would  venture  to  hope  that  the  work  will  be  found 
useful. 

C.  H.  H.  CHEYNE. 

ST  JOHN'S  COLLEGE, 

October,  1862. 


IN  the  SECOND  EDITION  comparatively  few  changes  have 
been  made.  The  work  has  been  revised,  and,  it  is  hoped, 
in  some  degree  improved.  The  Stability  of  the  Planetary 
System  has  been  rather  more  fully  treated,  and  an  elegant 
geometrical  explanation  of  the  formulae  for  the  secular 
variations  of  the  node  and  inclination  introduced,  for  which 
I  am  indebted  to  a  paper  by  Mr  H.  M.  Taylor,  Fellow  of 
Trinity  College,  in  the  Oxford,  Cambridge  and  Dublin 
Messenger  of  Mathematics. 

C.  H.  H.  C. 

1,  DEAN'S  YARD,  WESTMINSTER, 
September,  1870. 
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THE  PLANETARY  THEORY. 


CHAPTER  I. 

INTRODUCTION. 

1.  To  determine  the  motion  of  a  system  of  bodies  sub- 
ject only  to  their  mutual  attractions,  is  a  problem  the  mathe- 
matical difficulties  of  which  have  not  yet  been  overcome: 
hence  in  the  particular  cases  of  this  problem  which  Physical 
Astronomy  presents,  recourse  must  be  had  to  methods  of  ap- 
proximation.    Happily  the  arrangement  of  the  Solar  System 
renders   approximate   methods   possible,  and  in  the  skilful 
hands  of  the  Mathematicians  of  the  last  century,  they  have 
been  brought  to  a  high  state  of  perfection. 

2.  If  the  Sun  were  the  only  attracting  body,  the  planets 
would  describe  exact  ellipses,  agreeably  to  Kepler's  first  law ; 
but  in  consequence  of  the  attractions  of  the  planets  them- 
selves,  slight  deviations  from  elliptic  motion  are  produced. 
The  method  of  calculating  these  deviations,  to  which  our 
attention  will  chiefly  be  directed,  is  due  to  Euler,  but  was 
subsequently  extended  and  perfected  by  Lagrange :  it  con- 
sists in  supposing  the  planets  to  move  in  ellipses,  the  ele- 
ments (or  arbitrary  constants)  of  which  are  continually  though 
slowly  changing*. 

*  The  legitimacy  of  this  hypothesis  will  appear  when  we  come  to  treat 
of  the  equations  of  motion.     See  Arts.  21  and  22, 
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2  PLANETAEY  THEORY. 

3.  Now  the  elements  of  an  elliptic  orbit  are  (i)  the  mean 
distance,  or  semi-axis  major,  (ii)  the   eccentricity,   (iii)   the 
longitude  of  perihelion,  i.  e.  of  the  point  of  the  orbit  nearest 
to  the  Sun,  (iv)  the  longitude  of  the  epoch*,  or  mean  longi- 
tude at  the  epoch  from  which  the  time  is  reckoned,  (v).  the 
inclination  of  the  plane  in  which  the  orbit  lies  to  some  fixed 
plane  of  reference,  (vi)  the  longitude  of  the  ascending  node. 
Of  these  (i)  and  (ii)  determine  the  magnitude  of  the  orbit, 
(iii)  determines  its  position  in  its  own  plane,  (v)  and  (vi)  de- 
termine the  position  of  this  plane,  and  (iv)  has  reference  to 
the  position  of  the  body  itself  in  its  orbit. 

If  the  planets  moved  accurately  in  ellipses,  these  would  be 
constants  :  we  must  however  be  prepared  to  consider  them  as 
variable  quantities,  which  it  will  be  the  object  of  the  problem 
to  determine.  They  are  termed  the  elements  of  the  orbit. 

4.  But  further,  not  only  is  it  found  that  the  true  orbit 
of  a  planet  is  not  an  ellipse,  but  that  it  is  not  even  a  plane 
curve,  although  the  departure  of  the  planet  from  the  plane' 
in  which  it  is  at  any  instant  moving  is  extremely  slow.     We 
define  as  the  plane  of  the  orbit  the  plane  containing   the 
radius  vector  and  direction  of  motion  of  the  planet  at  the 
instant  under  consideration. 

5.  We  shall  suppose  the  Sun  and  planets  so  distant  from 
each  other  that  they  may  be  considered  to  attract  as  if  they 
were  condensed  into  their  respective  centres  of  gravity;  a  sup- 
position which  would  be  rigorously  true  if  these  bodies  were 
exactly  spherical,  and  either  of  uniform  density  or  composed  of 
concentric  spherical  shells,  the  density  of  each  shell  being  uni- 
form throughout.     The  errors,  however,  thus  introduced  into 
the  motions  of  translation  are  found  to  be  inappreciable  for  the 
planets,  though  not  in  the  case  of  their  satellites.     The  mo- 
tions of  rotation  will  not  be  considered  in  the  present  treatise. 

*  Also  briefly  termed  the  epoch. 
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6.  Moreover,  since  the  masses  of  the  planets  are  ex- 
tremely small  in  comparison  of  that  of  the  Sun,  it  follows 
that  in  cases  where  it  is  not  necessary  to  carry  the  approxi- 
mation beyond  the  first  order  of  these  masses,  we  are  permit- 
ted to  avail  ourselves  of  the  Principle  of  the  Superposition  of 
Small  Motions,  and  thus  to  reduce  the  problem  to  a  case  of 
that  of  the  Three  Bodies. 

7.  So  far  the  Theory  of  the  Planets  resembles  that  of 
the  Moon,  and  the  same  method  of  treatment  might  be  em- 
ployed in  both  cases.     But  they  differ  in  this  respect:  the 
ratio  of  the  distances  of  the  disturbed  and  disturbing  bodies 
from  the  central  one*  is  much  smaller  in  the  Lunar  than  in 
the  Planetary  Theory,  so  that  if  in  the  latter  theory  the  ap- 
proximation were  made  by  means  of  series  proceeding  by 
powers  of  this  ratio,  'it  would  be  necessary  to  retain  many 
more  terms  than  are  required  in  the  former.     On  the  other 
hand,  the  perturbations  of  the  Moon  are  far  larger  than  those 
of  the  planets,  since  in  the  former  case  the  Sun,  of  which 
the  mass  is  enormous,  and  the  distance  not  proportionately 
great,  is  one  of  the  disturbing  bodies.     For  these  reasons 
different  methods  of  calculation  are  employed. 

8.  To  find  an  expression  for  the  component  in  any  di- 
rection of  the  force  which  disturbs  the  motion  of  a  given  planet 
relatively  to  the  Sun. 

Let  M  denote  the  mass  of  the  Sun,  m,  m,  m",  &c.,  those 
of  the  planets,  and  suppose  the  relative  motion  of  m  re- 
quired. 

Let  x,  y,  z,  x ',  yt  /,  x",  y",  z",  &c.,  be  the  co-ordinates 
of  m,  m,  m".  &c.,  referred  to  any  system  of  rectangular  axes 

*  By  the  central  body  is  meant  that  whose  attraction  exercises  the 
greatest  influence  on  the  body  whose  motion  is  required ;  the  Sun,  for 
instance,  in  the  Theory  of  the  Planets,  and  the  Earth  in  that  of  the  Moon. 
All  the  other  attracting  bodies  are  called  disturbing  bodies. 

1—2 
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originating  in  the  centre  of  gravity  of  the  Sun ;  r,  r',  r",  &c., 
their  distances  from  the  origin ;  p,  p",  &c.,  the  distances  of 
m,  m",  &c.,  from  m. 


Now  if  to  every  body  of  the  system  we  apply  forces  equal 
and  opposite  to  those  which  act  upon  the  Sun,  we  shall  re- 
duce the  latter  to  rest  without  affecting  the  relative  motion. 
Hence,  considering  the  action  of  only  one  disturbing  planet 
m',  the  forces  acting  upon  m  will  be 

M  +  ra  .     ,.       .          ,r        i   •  n     P 
— 2 —  in  direction  mm,  or  briefly,  —z , 

m  .     ,.       .  , 

-TO  m  direction  mm , 
P 

m  .     ,.      _         ,,, 
-72  in  direction  m  M, 

of  which  the  last  two  constitute  the  disturbing  force. 

On  the  hypothesis  of  Art.  5,  the  components  parallel  to 
the  axes  of  the  forces  acting  on  m,  will  be 

__  JIM     m'  (xr  —  x)     m'x 
"^4     ~~^~      ~^r> 

_py  ^m '(y '-y)     m'y 

rfH         p'*          ~  r3   ' 
__  liz     m'  (z  —  z)     m'z 
~^^     ~p^~      "T7^3 
of  which  the  terms  containing  m  are  the  disturbing  forces. 
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Let  s  denote  the  length  of  the  arc  of  any  curve  measured 
from  some  fixed  point  up  to  m:  then  the  resolved  part  of  the 
disturbing  force  parallel  to  the  tangent  at  m  to  this  curve 
will  be 

m  (x  —  x)  dx     m'  (y  —  y]dy     m  (z  —  z)  dz 
"T75"^*4     ~  p3        5s        ~Y*        Ts 

fm'x  dx     my'  dy     mz'  dz\ 
~  \^~ds  +  ~r'*~~  ds  +  ~^  ds)' 


d 

or  -=- 

ds 


on  supposition  that  x,  y,  z,  are  alone  affected  by  the  process 

cJTt! 
of  differentiation,  which  may  be  written  —  y-  ,  if 


If  we  express  in  like  manner  the  disturbing  forces  due 
to  the  action  of  m",  m"',  &c.,  we  shall  have  for  the  whole 
component  in  this  direction 

air 


ds  "  ds 


, 
ds 

The  function  R  is  called  the  disturbing  function:  the 
arbitrary  curve  employed  above  will  be  termed  the  curve  of 
reference. 

9.     From  the  manner  in  which  -=-  has  been  introduced, 

ds 

it  appears  that  R  is  supposed  to  be  expressed  in  terms  of  s 
and  quantities  which  do  not  vary  with  s.     It  must,  however, 

be  borne  in  mind  that  in  -7-  the  variation  is  purely  hypo- 

thetical, and  has  nothing  whatever  to  do  with  the  actual 
variation  of  R  due  to  the  motion  of  the  planet. 
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For  example,  suppose  the  curve  of  reference  a  straight 
line  parallel  to  the  axis  of  x,  and  let  R  be  expressed  in 
terms  of  x,  y  and  z ;  then  in  this  case  x  only  will  vary,  and 
the  disturbing  force  parallel  to  the  axis  of  x  will  be  denoted 

by  -7- ,  y  and  z  being  considered  constant  in  the  differen- 
tiation. Similarly,  the  disturbing  forces  parallel  to  the  axes 
of  y  and  z  will  be  expressed  by  -y-  and  -y-  respectively,  the 

differential  coefficients  being  strictly  partial. 

Again,  suppose  the  carve  of  reference  a  circle  with  its 
plane  parallel  to  that  of  xy,  and  its  centre  in  the  axis  of  z, 
and  let  R  be  expressed  in  terms  of  the  polar  co-ordinates 
(rl}  6^)  of  the  projection  of  the  planet  on  the  plane  of  xy,  and 
its  distance  (z}  from  this  plane ;  then  in  this  case  Ol  only  will 
vary,  and  the  disturbing  force  perpendicular  to  the  projected 

/77? 

radius  vector  will  be  expressed  by  — -^- ,  rl  and  z  being  con- 
sidered constant  in  the  differentiation.  Similarly,  the  forces 
parallel  to  the  projected  radius  vector  and  to  the  axis  of  2, 

will  be  expressed  by  the  partial  differential  coefficients  -7- , 

arx 

dR  ..    i 

-T-  respectively. 

10.  The  disturbing  function  is  independent  of  any  par- 
ticular system  of  co-ordinates  that  may  be  employed.  For 

™     m      m  .     ,         f         ,. 
R  =  -^-^(xx  +  yy  +  zz) 

=  m      m'r  fxx      yy    |  zz\ 
p        r'2  \r  r       r  r'      rr'J 

m       m'r 

—  — T /«-  COS  ft), 

P        r 

if  o>  denote  the  inclination  of  r  to  r. 
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11.  To  express  E,'  in  terms  of  the  polar  co-ordinates  of 
the  projections  of  m  and  m'  on  a  fixed  plane,  and  of  their 
distances  from  it. 

Take  the  fixed  plane  for  tbat  of  xy  :  let  rlt  r/  be  the  pro- 
jections of  r,  r  upon  it,  and  0lt  #/  the  inclinations  of  rlt  rt'  to 
the  axis  of  x  ;  then 

x  —  rl  cos  Ol  ,       y—v^  sin  0lt 
x  =  r[  cos  #/,     y  —  r/  sin  #/; 
therefore    ##'  +  yy'  +  zz  —  rj\'  cos  (6l  —  0^)  +  a^', 


2  -  2r,r,'  cos  (^  -  ^' 


r   = 


Hence  by  substitution, 
ra' 


cos  (0l  -  0/)  +  (z  -  zj} 
m'r      cos  O- 


12.  In  a  subsequent  Chapter  we  shall  consider  the  de- 
velopment of  R  in  terms  of  the  time  and  the  elements  of  the 
orbit,  in  a  series  ascending  by  powers  and  products  of  the 
eccentricities  and  inclinations,  which  for  the  principal  planets 
are  very  small*.  At  present  we  shall  content  ourselves  with 
shewing  how  R  may  be  expressed  in  terms  of  these  quantities. 

*  For  the  smaller  planets  and  comets  this  is  not  the  case;  so  that 
different  methods  of  calculation  are  required  for  these  bodies. 
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We  shall  assume  that  the  equations  connecting  the  co-ordi- 
nates, the  time,  and  the  elements  in  an  elliptic  orbit,  hold 
also  when  the  motion  is  disturbed. 


13.  To  explain  how  R  may  be  expressed  in  terms  of  the 
time  and  the  elements  of  the  orbit. 

Let  r,  6  denote  the  radius  vector  and  longitude  of  the 
disturbed  planet,  the  latter  being  measured  on  a  fixed  plane 
of  reference  as  far  as  the  node,  and  thence  on  the  plane  of  the 
orbit  :  let  the  elements  be  a  the  mean  distance,  e  the  excen- 
tricity,  OT  the  longitude  of  perihelion,  e  the  longitude  of  the 
epoch,  (the  last  two  being  measured  in  the  same  way  as  #.) 
O  the  longitude  of  the  node  measured  on  the  plane  of  refer- 
ence, and  i  the  inclination  of  the  plane  of  the  orbit  to  the 
plane  of  reference.  Our  object  is  to  express  R  in  terms  of  t 
and  these  elements. 

Again,  let  #0,  OTO,  e0,  fi0  denote  the  longitudes  of  the 
planet,  of  perihelion,  of  the  epoch,  and  of  the  node,  measured 
entirely  on  the  plane  of  the  orbit. 

Let  a  sphere  be  described  with  its  centre  coinciding  with 
that  of  the  Sun,  and  its  radius  of  any  magnitude  :  let  the 
planes  of  reference  and  of  the  orbit  cut  it  in  the  great  circles 
NM,  NP,  then  the  lines  of  nodes  will  cut  it  in  N;  let  the 
radius  vector  of  the  planet  cut  it  in  P,  the  projection  of  this 
radius  on  the  plane  of  reference  in  M,  and  the  lines  from 
which  0,  Q0  are  measured  in  L,  0  respectively.  We  shall 
suppose  L  to  be  the  same  origin  as  that  from  which  0A  is 
measured  in  Art.  11. 


Then  in  the   figure   LM=0lt  LN+NP  =  6 
LN=£l,  the  angle  PNM=i,  and  P31  =  the  latitude  of  the 
planet  which  .we  shall  denote  by  A,. 
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Hence  from  the  right-angled  triangle  PNM, 

tan(^-H)  =  cosi  tan((9-  O) (1), . 

sin  X  =  snu  sin  (0- ft) (2): 

also  i\  =  r  cosX (3), 

z  —  r  sinX (4). 


Again,  from  the  formulae  of  elliptic  motion*, 
r  =  a  {1  +  Je2  —  e  cos(nt  +  e0  -  OTO)  -  Je2  cos2  (nt  +  e0  —  w0)  -.  .  .}, 
00  =  n£  -f  e0  +  2e  sin(7tf  +  e0  -  *TO)  +  ^e2  sin2  (72*  +  e0  -  -or0)  +.  .  .f: 
but  6-eQ  =  LN-  ON  =  e-e0  =  OT-^; 

therefore  0=00  +  e-e0,     e0  -  i!T0  =  e  —  -cr, 

and  our  formulae  become 


0=  w<  +  e  +  2e  sin  (ytf  +  e  -  «r)  +  |e2  sin  2  (nf  +e—  BT)  +  ......  (6). 

In  Art.  11  we  have  expressed  .R'  in  terms  of  rlt  6^  and  z\ 
hence  by  equations  (1)  to  (4)  it  may  be  expressed  as  a  func- 
tion of  r,  0,  O,  and  i  :  we  may  then  substitute  for  r  and  0 
from  equations  (5)  and  (6),  and  R'  will  be  expressed  in  terms 

*  See  Tait  and  Steele's  Dynamics,  Art.  162,  Fifth  Edition. 
t  n  is  termed  the  mean  motion,  and  is  connected  with  the  mean  distance 
by  the  equation  n2a3  =  ju,  where  /j.=M+m. 
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of  t  and  the  elements  of  the  orbit.  Similarly  R"  ',  jR'",  &c., 
and  therefore  R  may  be  expressed  in  terms  of  t  and  the 
elements. 

14.  We  proceed  to  investigate  certain  relations  which 
subsist  between  the  partial  differential  coefficients  of  R  with 
respect  to  the  co-ordinates  of  the  disturbed  planet,  and  its 
partial  differential  coefficients  with  respect  to  the  elements 
of  the  orbit.     These  will  be  useful  in  obtaining  the  formula? 
by  which  the  values  of  the  elements  are  calculated. 

We  premise  that  when  we  speak  of  the  partial  differ- 
ential coefficient  of  R  with  respect  to  one  of  the  elements, 
we  suppose  R  expressed  in  the  manner  indicated  in  the  last 
Article,  and  that  the  time  as  well  as  the  other  elements  are 
considered  constant  in  the  differentiation  :  when  we  speak  of 
the  partial  differential  coefficient  of  R  with  respect  to  r  or  0, 
we  suppose  R  expressed  in  terms  of  r,  6,  i  and  O,  which  may 
be  done  by  equations  (1)  to  (4). 

„  dR     dR     dR 

15.  To  shew  that  -^  =  •*-  +  -3—  . 

dO       de       dtxr 

Equations  (5)  and  (6)  of  Art.  13  may  be  written 

r=f(nt  +  e  —  «r), 
0  —  w  =  (f>  (nt  +  e  —  -sr), 

whence  it  follows  that 


d0     d&_:=l 
de      diz 

Now  since  e  and  «r  enter  into  R  only  through  r  and  6, 

dR^_dRdr     dR  dO 
de      dr    de      dd   de  ' 
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dRdr^     dR   <W  t 
dr    dm      de   d-rs  ' 


therefore,  by  addition, 


dR     dR  =  dR 

de      dtp      de  ' 


^   .  dR     dR     dR  ,  dR 
16.     To  shew  that  -^-  =  -j-  +  T-  +  TO  • 
d^      de       der      dii 

From  equations  (1)  to  (4)  of  Art.  13  we  obtain 

^  =  <f>  (r,  e  -  a,  i), 


where  (f>,  %,  -^r  are  symbols  of  functionality. 
It  follows  that 


Now  since  by  Art.  11,  R  is  a  function  of  rv  0lt  and 
dRdR  dt       dR  d        dR  dz 


dR^dR  dr,      dR  de,      dR  dz  % 
<m  ~  dr,  d&      de,  d&     dz   dll ' 

therefore,  by  addition, 

dR     dR  ^  dR 

7/1       I        Jf*\     ~         JA      ) 

ctu      all      at/* 
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whence,  by  the  last  article, 

dR_dR     dR     dR 

~  +  ~      ' 


1(_      ^     ,.   .    dR  .  ,dR       ,  dR 

I/.     To  obtain    ,    in  terms  of-*-  and  ,Q  . 
de  dr  d# 

If  u  denote  the  excentric  anomaly,  we  have* 

r  —  a  (I  —  e  cosu)  ....................  (1), 

6  -  >& 


(3), 


from  which  r  and  6  may  be  expressed  in  terms  of  t  and  the 
elements  by  eliminating  u.     Assuming  r  and  6  so  expressed, 

,  ,      ,  ,   .    dr       ,  dO 
we  proceed  to  obtain  -j-  and  -y-  . 
de         de 


-P,         f-,\    dr         f     .      du  \ 

±*  rom  (1),  -^-  =  a  f  e  sin  u  -= cos  w  J , 


and  from  (3),  -r-  (I  —  e  cosu)  —  shut  =  0  ..................  (4)  ; 


eliminating      -y-  ,  we  have 


C  e  —  cosu  ] 

—  a  }  i —         \t  r  > 


=  a 


—  e  cosu  —  e 


sn 


*  See  Tait  and  Steele's  Dynamics,  Arts.  160  and  161,  Fifth  Edition. 
This  use  of  the  auxiliary  u  was  first  suggested  by  Mr  W.  Pirie,  Camb.  Math. 
Journal,  1837,  Vol.  i.,  p.  47. 
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1  +  tan2 


2 

=  —  a  cos  (0  —  «•). 
Again,  differentiating  the  logarithms  of  equation  (2), 

__1__^=1  /_J_         1    \         1     eZu 
sin  (0-  w)  cfo      2  \1  +  e  +  1  -  e)  +  sin  u  de  ' 

eliminating  -=-  by  means  of  (4), 
de 


sn     —  vre        —  e         —  e  cos  w 


if  /i,2  =  /^a  (1  —  e2)  ;  therefore 

-T-  =  a  (^2  +  -  J  sin  (^  —  /S3")- 
de         \h?     r) 

Now  since  R  is  a  function  of  e  only  because  it  is  a  function 
of  r  and  6, 


7  7->  /  1  \ 

-  a  cos  (V  -  *r)  -r-  +  a  ( &  +  -  }  sin  (6  - 


N  dR         /LL      1\ 
—  «r)  -J7  +  a  (p  +  -  ) 

$/7  X/i1          /  / 


Since  0—  w=  #0  —  w0,  (see  Art.  13,)  this  equation  may 
be  written 

dR  dR         ffi      1\   .  ^ 

5-  =  -acos(00-^0)^  +  a^  +  -Jsm(00-OTo)^, 

under  which  form  it  will  be  useful  in  the  next  Chapter. 


CHAPTER  II. 

FORMULA   FOR  CALCULATING  THE  ELEMENTS   OF  THE  ORBIT. 

18.  WE  now  proceed  to  form  equations  of  motion, 
taking  the  Sun's  centre  for  the  origin  of  co-ordinates,  the 
radius  vector  of  the  planet  for  the  axis  of  sc,  a  perpendicular 
to  it  in  the  plane  of  the  orbit  for  the  axis  of  y,  and  a  normal 
to  this  plane  for  the  axis  of  z.  With  this  system,  it  will 
be  shewn  that  two  of  the  resulting  equations  can  be  ex- 
pressed in  the  same  forms  as  if  the  planet  moved  in  one 
plane. 

Let  #,  y,  z  be  the  co-ordinates,  u,  v,  w  the  velocities  of 
the  planet  with  reference  to  three  rectangular  axes  originat- 
ing in  the  Sun's  centre,  and  moving  with  angular  velocities 
c^,  <£2,  <£3  about  their  instantaneous  positions  :  let  X,  Y,  Z\>e 
the  accelerations  due  to  the  impressed  forces  in  the  directions 
of  the  axes.  Then  (Routh's  Rigid  Dynamics,  Arts.  244,  245 
Third  Edition), 

_dx_  1 


dz 

w^-fa-vh 

and  the  equations  of  motion  are 
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du 


.(2). 


dv 


„     dw 

~W 


In  these  equations  c^,  <f>2,  <p3  are  arbitrary;  we  propose 
so  to  determine  them  that  the  axis  of  x  may  coincide  with 
the  radius  vector  of  the  planet,  and  the  plane  of  xy  with  the 
plane  of  the  orbit. 

In  order  that  the  axis  of  x  may  coincide  with  the  radius 
vector  of  the  planet,  we  must  have 

always ;  and  therefore 

dx  _dr     dy  _         dz  _ 
~dt~dt'    dt''     dt~ 

and  in  order  that  the  plane  of  xy  may  coincide  with  the  plane 
of  the  orbit,  we  must  have 

always ;  and  therefore 

dw 


Hence  equations  (1)  give 

=  di' 
and  equations  (2)  become 


dr 
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7/1 

-y 


7/1 

Now   let    cj>3  =  -yj°  :   thus   for   the  equations   of  motion 


we  have 


cFr        AWA* 
~d('     '  \~dtl  ' 


rdt       dt 


j' 


of  which  the  first  two  are  the  same  in  form  as  if  the  plane 
of  the  orbit  were  at  rest. 

19.  If  we   measure   on   this   plane  from   the   planet's 
radius   vector   in   a   direction   contrary   to  that  of  motion 
an  angle  equal  to  #0,  we  arrive  at  what  may  be  considered 
as  the  origin  from  which  00  is  measured.     Since  this  will  be 
a  point  having  no  angular  velocity  about  the  axis  of  z,  which 
is  normal  to  the  plane  of  the  orbit,  it  is  said  to  be  fixed  in 
the  plane  of  the  orbit*. 

20.  We  shall  for  the  present  confine  our  attention  to  .the 
first  two  of  the  above  equations. 

In  order  to  find  the  components  of  the  disturbing  force 
parallel  and  perpendicular  to  the  radius  vector  of  the  planet, 
let  us  take  first  the  radius  vector  as  the  curve  of  reference  ; 
then  s  =  r,  and  R  being  supposed  expressed  as  a  function  of 
r,  0,  O,  and  i  (see  Art.  13),  we  have 

dR^dRdr     dRdO     dRdn     dRdi^ 
ds      dr  ds      dQ  ds     d£l  ds       di  ds 

_dR 
"  dr' 

*  It  is  necessary  to  define  the  expression  fixed  in  the  plane  of  the  orbit, 
since  the  definition  of  Art.  4  is  not  sufficient  completely  to  regulate  the 
motion  of  this  plane,  though  affording  it  a  distinct  geometrical  position. 
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since  0,  H  and  i  do  not  vary  with  s.     Hence  the  disturbing 

force  in  direction  of  the  radius  vector  =  -=— . 

dr 

Again,  let  us  take  as  the  curve  of  reference  a  circle  in  the 
plane  of  the  orbit,  with  its  centre  coinciding  with  that  of  the 
Sun ;  then  Ss  =  rW,  and  we  have 

dR_\dR 
ds  ~r  d0> 

since  r,  fi,  and  i  do  not  vary  with  s.     Hence  the  disturbing 

force  perpendicular  to  the  radius  vector  •—  — ^  . 

r  dd 

We  have  then 


r^  dr'         ~r  d0} 

and  the  equations  become 


df 


21.  These  equations  do  not  admit  of  rigorous  integra- 
tion, but  we  may  reduce  them  by  the  method  of  the  Varia- 
tion of  Parameters  to  a  system  of  differential  equations  of  the 
first  order.  The  principle  of  this  method  may  be  explained 
as  follows.  Suppose  it  required  to  integrate  the  equations 

.   /         .  dx     dy     dfx 

Ai*w*;i5'  ^'  ~/P» 

\  QjU         \Mi          (Mi 

dx     dy     d*x 
V'  '~dt}    c5'    ~dt*'    df 

where  Px,  P2  are  functions  of  t.     The  solution  of  these  equa- 
tions can  be  made  to  depend  upon  that  of  the  equations 
C.  P.  T.  2 


•(i), 
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<^i  =  0,  <£2  =  0.     Suppose  the  four  first  integrals  of  these  last 
equations  to  be 

..  dx 


dx 


dy\ 

w 


•(ii), 


dx 


where  c1?  c2,  cg,  c4  are  arbitrary  constants  or  parameters. 
The  method  of  the  Variation  of  Parameters  consists  in  so 
determining  ct,  c2,  C3  and  c4  as  functions  of  ^  that  these  inte- 
grals f  and  therefore  the  two  final  integrals  of  the  equations 
<jf>i  =0,  </>2  =  0,  which  can  be  obtained  from  equations  (ii)  by 

Ci  ?*  fiii\ 

eliminating  -j-   and    •—}   shall   satisfy  equations  (i).     That 

cp  C2,  c3,  and  c4  can  be  so  determined,  may  be  seen  as  follows: 
by  the  solution  of  equations  (i),  values  of  x  and  y  and  there- 
fore of  -jz  and  —  can  be  found  as  functions  of  t  and  constant 
dt          dt 

quantities ;  if  these  be  substituted  in  equations  (ii)  the  requi- 
site values  of  c1?  c2,  c3  and  c4  will  be  obtained.  For  an  ex- 
ample of  the  application  of  this  method,  see  Boole's  Differen- 
tial Equations,  Chap.  IX.  Art.  11. 


22.     If  in  equations  (1)  and  (2)  of  Art.  20  we  put  R  =  0, 
and  then  integrate  them,  we  obtain 

1      fju 
T     h 

dr     lie  . 


•(3), 


•(4), 
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where  h,  e,  OTO  are  the  constants  of  integration.  Equation  (3) 
indicates  motion  in  an  ellipse,  of  which  e  is  the  excentricity, 
^  the  longitude  of  perihelion,  and  h  twice  the  area  described 
in  an  unit  of  time.  If  the  mean  distance  in  this  ellipse  be 
denoted  by  a,  we  have  in  addition 

A2  =  /m(l-e2)  ........................  (6). 

We  shall  assume  (in  accordance  with  the  principles  of  the 
method  of  the  Variation  of  Parameters)  the  first  and  second 
integrals  of  equations  (1)  and  (2),  together  with  equation  (6), 
to  retain  the  same  forms  when  R  is  not  zero  ;  h,  e,  ^Q)  and  a 
being  in  this  case  considered  variable*. 

The  values  of  these  elements  are  to  be  obtained  from  the 
condition  that  the  above  integrals  shall  satisfy  equations  (1) 
and  (2). 

23.  If  their  values  as  calculated  for  any  given  time  be 
substituted  in  equation  (3),  it  will  represent  an  ellipse  having 
a  contact  of  the  first  order  with  the  actual  orbit,  since  the 

values  of  -^~  and  -j-9  at  the  common  point  will  be  the  same 
at          at 

for  both  curves.  It  is  termed  the  instantaneous  ellipse,  since 
the  planet  may  for  an  infinitely  small  time  be  supposed  to 
move  in  it.  Moreover,  the  velocity  and  direction  of  motion 
of  the  planet  will  be  the  same  as  if  it  moved  in  this  ellipse, 
so  that  if  at  any  time  the  disturbing  force  were  to  cease,  the 
planet  would  continue  to  move  in  the  instantaneous  ellipse 
constructed  for  that  time.  This  is  accordingly  sometimes 
given  as  the  definition  of  the  instantaneous  ellipse. 

*  We  shall  also  for  convenience  suppose  the  equation  n2a3=/*  to  hold  in 
the  disturbed  orbit,  n  being  of  course  considered  variable. 

2—2 
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24.  To  obtain  formulae  for  calculating  the  elements  of  the 
instantaneous  ellipse  at  any  time. 

Suppose  the  value  of  c  required,  where  c  denotes  any  one 
of  the  elements.  From  equations  (3),  (4),  and  (6)  we  may, 
by  eliminating  the  other  elements,  obtain  c  as  a  function  of 

r,  00,  -77  ,  and  h*:  let  then 
at 

c  =f(r,  00,  r',  h\ 

where  r'  is  written  for  -7-  .     Differentiating,  we  have 
at 

dc^dfdr      dfdO,     d£  dr      df  dh 
~dt  ~~  dr  dt  +  d00  dt  +  dr'  dt  +  ~dh  dt  ' 

dr      tfr        (dO\z      fi      dR 
N°W         ==r--*+> 


d  /  9  ddA      dR    p  , .       ,_N 

=  fTt  V     "///  J  =  77^  '  from  equation  (2)  ; 

.       „        dc     d/  dr      d/  d00      df  ( 
therefore  -^  =  ?-  -^  +  ^r  -^  +  -rnr-jr ,    7j  ,      ^ 


d/  dr     dh  d6' 

But,  since  by  hypothesis,  if  ^  were  zero  and  c  constant, 
our  assumed  integrals  would  still  satisfy  the  differential  equa- 
tions, we  have  (making  R  zero  and  c  constant) 

dfdr      d/d00      df(    (Mtf      »} 
~  +  dO    dt  +  dr'         dt     ~ 


,  dc      df  dR     dfdR 

therefore  -j-  =  -f-,  -=-  +  ~  -j^  . 

d£      dr  dr      dA  dd 


We  retain  ^  for  convenience,  in  preference  to  replacing  it  by  r2  -r-°  . 
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Hence  in  obtaining  the  formulae  for  calculating  the  ele- 
ments of  the  orbit,  we  may  proceed  as  follows.  From  equa- 
tions (3),  (4),  and  (6)  we  may  express  the  element  required 

dr 

as  a  function  of  r,  00,  -r  ,  and  h.  We  may  then,  by  differen- 
tiating the  resulting  equation  with  respect  to  t  as  if  r  and  00 

dR  f     d*r         ,   dR  ,     dh 
were   constants,  writing   -y-  for  -^ ,   and    -^  for  -=-  ,   and 

dr 

eliminating,  if  necessary,  -7-  and  Ji  by  means  of  equations 

dt 

(4)  and  (6),  obtain  the  differential  coefficient  of  the  element 
required  in  terms  of  the  elements,  the  co-ordinates  of  the 
planet,  and  the  disturbing  force.  The  result,  however,  will 
in  every  case  admit  of  being  expressed  in  terms  of  the  ele- 
ments and  of  the  differential  coefficients  of  R  with  respect  to 
them,  a  form  under  which  it  is  very  convenient  of  applica- 
tion. 

25.  To  obtain  a  formula  for  calculating  the  mean  dis- 
tance. 

From  equations  (3),  (4)  and  (6),  if  e  and  -CTO  be  eliminated, 
we  shall  find 


a      \dt 

7  D 

Differentiating  as  if  r  were  constant,  and  writing  -r-  for 

cLr 

dR  „     dh 


dR  dO\ 

a*  dt~       drdt      dd   dt)' 
Now  since  r  =f  (nt  +  e  —  -cr), 

^o  -  OTo  =  0  -  v  =  ^  (n  t  +  e  -  «• ) , 

and  the  forms  of  -7-  and  -~  are  the  same  as  if  the  elements 
dt          dt 

were  constant,  we  have 
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dr        , ,        dr 


,    .    „    ,  dO,        d90        dO 

and  similarly,  -r  =  n  -~  =  n  -7-  : 

J  dt          de         de 


therefore  C  T  =  2n  f  ^  "  +  ^f  ^ 

a#  ae> 


c?a     2na*dR 

or  -y-  = =-  . 

at         yi6    ae 

26.  This  formula  may  also  be  obtained  as  follows.  If  5 
denote  an  arc  of  the  actual  path  of  the  planet  measured  from 
some  fixed  point  to  its  position  at  time  t,  we  have  the  equa- 
tion of  motion 

d*s         u,  dr     dR 


and  by  a  known  formula 


dt     "  r       a' 

Differentiating  the  latter  we  obtain 
ds  d?s  __     2/i,  dr     fi  da 


ds 
and,  multiplying  the  former  by  2  -5-  , 


=_  -. 

dt  d?~  "  r*  dt*     ds  dt'" 

.-.f  a  da      ^  dR  ds 

therefore  -^- 


where  —  ^  —  denotes  the  differential  coefficient  of  R  with 
at 
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respect  to  t,  only  so  far  as  R  involves  t  through  involving  the 
co-ordinates  and  elements  of  the  disturbed  planet.  Now  since 
the  velocity  of  the  planet  at  any  time  can  be  expressed  in 
terms  of  the  co-ordinates  and  elements  of  the  instantaneous 
ellipse  constructed  for  that  time  in  the  same  form  as  if  it 
moved  in  this  ellipse,  its  component  in  any  direction  can  also 
be  so  expressed.  Hence,  considering  R  as  a  function  of 

/i       /  ,  -i  \    1 1  r>  dr.    dd,    dz 

rlt  0t,  z  (see  Art.  11),  the  values  of  -~  ,  -~ ,  -y-,  and  there- 
of     dt    dt 

fore  of  — i-^  ,  may  be  expressed  in  the  same  forms  as  if  the 

elements  were  invariable.  Since,  then,  t  always  occurs  in  R 
coupled  with  e.in  the  expression  nt  +  e,  we  have 

d  (R)  dR  dR 

=  n 


therefore 


dt          d  (nt  +  e)         de  ' 
da     2na2  dR 


dt        /j,     de 


27.     To  obtain  a  formula  for  calculating  the  eccentricity. 
From  equations  (3)  and  (4),  if  «-0  be  eliminated,  we  obtain 

/^2  =  ^      //>_M2 
(dt)  ~~  hz       (r     h)  * 

7  I) 

Differentiating  as   if  r  were  constant,  and  writing  -T- 

dr 


^?  —  =  _  ,  _ 

dr  dt      h2  dt      (h?       \r  +  h2)  (r     h)}  dt 

_tfede      j*I-<?dR     h  dR 
=  + 


dt        ~W~     dd 


dh     dR 
= 
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therefore   &<*e_dRdr     dRM,   .^-e*)dR 
h?  dt~  dr  dt  +  dB   dt  '          h3        dO 

fdRdr     dRd0\_  ^  (1  -  ez]  dR 
n(drde+d0de)  h3         d0  ' 

f  de     ntfdR     l-ez/dR     dR\ 

therefore          -7-  =  —  -^  ---  ,  —  [  -j-  +  3—    ,     (Art.  15), 
dt     fSe  de         he    \  de      dvrj 


=  na(l-  e*}  dR     na  V(l  -  ez)  fdR     d 
fjie        de  pe         \de      d 

since  A2  =  /xa  (1  —  e2),  and  nV  =  p. 

28.     This  formula  may  also  be  deduced  from  that  of  the 
mean  distance,  by  means  of  the  equation 
#>  =  ^(l~e2). 

™-    i,         01  dh         ,.,        2,  da     ~        de 
We  have  2^=^(1-^-2^^, 

dR  dR      ,.        de 


-e*)  /dR     d 
'    --      - 


or  -      -  -       -^ 


therefore 


-r  ---  —  -- 
de  fie        \ 


dR     dR\ 

-y-  +  -7-    . 
de      dixj 


29.     To  obtain  a  formula  for  calculating  the  longitude  of 
perihelion. 

From  equations  (3)  and  (4),  if  e  be  eliminated,  we  obtain 
dr     ,  ,  h      IL 


Differentiating  as  if  r  and  00  were  constant,  and  writing 

,     d?r  ,,  . 

for  -j-s-  ,  we  obtain 
at 

dr  .  fc  ,  ,.          v  c?E      /I      At\  dh 
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,i_       /.  dr  ,n          .  dnrn 

therefore         -y-  cosec  (00  —  TTO)  —v^ 


but  from  equation  (4)  -=-  cosec  (0Q  —  OTO)  =  y  , 

fc0       h   dR     na*J(l-J)dR 
therefore          -°  =  -   -  —  =  -  —  /  —  . 

at      pea  de  pe         de 

Now  if  is  denote  the  longitude  of  perihelion  measured  on 
the  plane  of  reference  as  far  as  the  node,  and  thence  on  the 
plane  of  the  orbit,  H  the  longitude  of  the  node  on  the  plane 
of  reference,  H0  its  longitude  on  that  of  the  orbit,  we  have 

^-OTo  =  n-n0; 

,,       f  d™     dv?0     d£l 

therefore  --  =  -     +  -yr  - 


-     . 
dt       dt       dt       dt 

Now  -37  is  the  angular  velocity  of  the  line  of  nodes  on 
dt 

the  plane  of  reference,  —  ~  its  angular  velocity  on  the  plane 

of  the  orbit  ;  and  since  the  plane  of  reference  is  fixed,  the 
former  is  the  total  angular  velocity  :  hence 


r  * 

therefore  = 


or,  substituting  for  -=—  from  Art.  31  or  35, 
dt 


dt  ~  e         de 


2v  TT>       wa  tan  ^    7ri 
—e    dR  2  dR 
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To  obtain  formulae  for  calculating  the  longitude  of  the 
node,  and  the  inclination. 

30.     We  now  return  to  our  third  equation  of  motion, 

Hf!     *t-v...    ' 

or,  as  it  may  be  written, 

hfa  =  Zr  ........................  (1). 

We  have  seen  (Art.  18),  that  <£2  =  0;  hence  the  motion  of 
the  plane  of  xy,  which  coincides  with  the  plane  of  the  orbit, 
is  compounded  of  the  angular  velocities  fa  about  the  axis 
of  x,  and  <£3  about  the  axis  of  z.  Now  the  former  is  equiva- 
lent to  an  angular  velocity  fa  cos  (0  —  H)  about  the  line  of 
nodes,  and  an  angular  velocity  fa  sin  (6  —  O)  about  an  axis 
perpendicular  to  it  in  the  plane  of  the  orbit  :  but  the  angular 

velocities  of  the  plane  of  the  orbit  about  these  axes  are  -j-  and 
sin  i  -j-  respectively  ;  therefore 

fa  cos  (6  —  fl)  =  -77  ,     fa  sin  (6  —  fl)  =  sin  i  -7-  . 
at  at 

Hence,  by  equation  (1), 

n)  ..................  (2), 


(3). 


31.  In  order  to  determine  Z  we  must  suppose  the  curve 
of  reference  perpendicular  to  the  plane  of  the  orbit.  If  we 
denote  by  5  an  arc  of  this  curve  measured  from  some  fixed 

7  r> 

point  up  to  the  planet,  we  have  by  Art.  8,  Z—  -r-  .     Now 

in  the  same  way  that  the  position  of  the  planet  is  known 
when  we  know  r,  0,  i  and  fl,  the  position  of  any  point  on 
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the  curve  of  reference  may  be  determined  by  its  polar  co- 
ordinates r,  0  on  a  plane  passing  through  it  and  the  Sun,  the 
inclination  i  of  this  plane  to  the  plane  of  reference,  and  the 
longitude  fl  of  its  node.  Since,  however,  an  infinite  number 
of  planes  can  be  drawn  through  two  given  points,  we  must 
introduce  some  further  condition  to  fix  the  position  of  that 
on  which  r  and  6  are  measured*. 

Different  forms  of  expression  will  be  obtained  for  Z 
according  as  different  conditions  are  assigned.  First  sup- 
pose the  plane  to  pass  through  SN,  the  planet's  line  of 
nodes  :  let  P  be  the  position  of  the  planet,  draw  PC  per- 
pendicular to  SN,  and  take  for  the  curve  of  reference  a 
circle  AP  with  centre  C.  Through  SN  draw  a  plane  in- 
clined at  a  small  angle  Si  to  the  plane  of  the  orbit,  cutting 
the  circle  AP  in  p,  and  let  Pp  =  Ss ;  then  since  Pp  is  per- 
pendicular to  the  plane  SNP, 


or 
therefore 


Pp=CP.Si, 
$s  =  r  sin  (0  -  fl)  S 

di  _  I 

ds~  r  sin  (0  -  12)  ' 


See  Appendix,  Art.  7. 
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Now,  (see  Art.  13),  R  may  be  expressed  in  the  form 

in  which,  if  the  series  for  r  and  0  be  substituted,  R  will  be 
expressed  as  a  function  of  t  and  the  elements.     Hence 

dR_dRdr     dR  d(^     dR  dti     dRdi^ 
ds      dr  ds      d0  ds     d£l  ds       di  dsy 

dr^d0^dD,_d^_  1 

~7/I  ~~  ^>      ~J~~  ~ 


ds  ds  ds  ds     rsin(0-H)' 

therefore  -r  =  — : — ia      ~\  -j-r . 

as      r  sin  (0  —  O)   di 

On  substituting  this  value   for   Z  in   equation  (3),   we 
obtain 

7/""\  J  T) 

7    .      .  d,\L      aJt 
n  sin  i  —j-  =  — JT  , 

or  since  A2  =  pa  (1  —  e5),  and  w2a3  =  /x, 

6?O  na  dR 

dt      fju  v  (1  —  e2)  sin  i  di 

32.  By  substituting  the  value  of  Z  found  in  the  pre- 
ceding Article   in   equation    (2)   of  Art.  30,  we   might,   of 

di 
course,  obtain  an  expression  for  -=- ,  but  this  would  involve  0, 

and  consequently  be  in  a  form  inconvenient  for  calculation. 
We  proceed,  then,  to  obtain  an  expression  for  Z  by  means  of 
which  0  may  be  eliminated. 

33.  Suppose  the  plane  on  which  r  and  0  are  measured, 
instead  of  passing  through  SN,  to  pass  through  a  line  SC 
in  the  plane  of  the  orbit  perpendicular  to  £ZV":  let  P  be  the 
position  of  the  planet,  draw  PC  perpendicular  to  SC,  and 
take  for  the  curve  of  reference  a  circle  AP  with  centre  C. 
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Through  SC  draw  a  plane  inclined  at  a  small  angle  to  the 
plane  of  the  orbit,  cutting  the  sphere  in  the  great  circle  nmp, 
and  the  circle  AP  in  p.  Draw  Nm  perpendicular  to  np. 


Then  Nm  which  measures   the   inclination   of  the   two 
planes  =  —  SO  sin  i,  and 

<7P=rcos(<9-O); 
hence  if  Pp  =  Ss, 


therefore 
Again 

therefore 


=  -  r  cos  (0  -  O)  SO  sin  i; 
SO_  -1 

ds  ~  r  cos  (0  —  O)  sin  i  ' 


=  nm  —  nN 
=  -  SO  cos  i  +  SO  ; 
dO  .  dtl 


=  0         - 


Hence  -j-  =  -^  j-  +  -j^ 


_ 
"       S 
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On  substituting  this  value  for  Z  in  equation  (2),  we  have 
di  1       (dR  ..dR 


k  sin  i 


dR     «  .  9i  fdR  .  ^.u   ^Art  15^ 


+  " 


sin  i  d£l  2 

The  only  remaining  element  is  the  epoch,  but  before  pro- 
ceeding to  obtain  a  formula  for  its  calculation,  we  shall  give 
another  method  of  obtaining  the  results  of  Arts.  31  and  33. 

34.  To  obtain  a  formula  for  calculating  the  inclination. 
(Second  method.) 

•  If  the  motion  of  the  planet  be  referred  to  the  polar  co- 
ordinates of  its  projection  on  the  fixed  plane  of  reference  and 
its  distance  from  this  plane,  we  have  the  equation 


r^dt\l    dt)      rlddl 


Now  if  §A,  SA1  denote  the  vectorial  areas  swept  out  in 
the  time  $t  on  the  plane  of  the  orbit  and  the  plane  of  reference 
respectively,  we  have 


but  SA^r^,  $A= 

«  «  dO.       9  dOn        .      ^ 

therefore  r*  -^  =  r  -y*  cos  i  =  h  cos  i. 
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Hence  our  equation  of  motion  becomes 
d  .,         .      dR 


,    .    .  di  .dh     dR 

or  —  h  sm  i  -j-  +  cos  i  -j-  =  -™  ; 

dt  dt     dO^ 

.  di     dR  .  dh 

therefore  —  h  sin  i  -j-  —  -j^  —  cos  i  -=- 

dt     d&i  dt 

dR  .dR 


,  (Art.  16), 


,       di  1      (dR     ,-  ..  (dR     dR 

therefore  -^  =  —  -r— -. — H  -77^  +  (1  -  cos  i)   -=-  4-  -r- 

/V/"  n    Q1  Tl  Q     I  ft  \/  \    ft  G  SITT 


na         (  1     dR          i  (dR     dR\] 

—, o  <—. — -. 1-  tan-  [ 1 v . 

V  (1  —  e  )  (sin  i  all  2  \  ae      atjr/ ) 


35.     To  obtain  a  formula  for  calculating  the  longitude  of 
the  node.     (Second  method.) 

We  have  by  Art.  13, 

5=/(r,0,fl,0 (1), 

or,  since  evidently  6  —  fl  =  #<,  -  H0, 

-B=/(r,00  +  n-n0,n,;)* (2), 


*  We  have  made  this  transformation,  because,  although  the  value  of  —  ° 
is  the  same  in  form  as  if  the  elements  were  invariable,  this  is  not  the  case 
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Now  we  have  seen  (Art.  26)  that     ^       (denoting  by  it 

that  R  is  to  be  differentiated  with  respect  to  t  only  so  far 
as  it  involves  t  through  involving  the  position  of  the  dis- 
turbed planet)  may  be  expressed  in  the  same  form  as  if  the 
elements  were  invariable. 

We  have  then,  considering  the  elements  variable, 
d(R)_dRdr     dRd(00  +  Q-Q,9)     dRdti     dRdi 
~~dT  ~dr~dt  +  deo~     ~dt~        ^  d&~dt  +  ~di  dt' 

and,  considering  them  invariable, 

d(R)  _dRdr     dR  d^ 
~~di        ~d/r  ~dt  +  d00  ~dt  ' 

Equating  the  two  values  of    ^   '  ,  we  obtain 

Cit 

dR  AW!  __  d&0\  .d^d^.dRdi^r, 
dOQ\di  ""  dt  )     d£l  dt  +  di  dt~ 

dR  .     ,ON      dR  .    n,      dR  f  dR 
Now  m2=inl=+, 


and  (see  Art.  29),         =        cos  i; 

(dR  ,  (dR  ,  dR\\  dO,     dR  di 

therefore          +  (1  -  cos  .)          +  -     +  -         =  0. 


di 
Substituting  for  -    its  value, 


d£l  _  _  na  dR 

dt  ~~  j>      l  "~  g2  sm  *  ^* 


36.  To  obtain  a  formula  for  calculating  the  longitude  of 
the  epoch. 

If  R  be  expressed  in  terms  of  t  and  the  elements  (see 
Art.  13),  since  nt  +  e  always  occurs  as  one  symbol,  we  may 
write 

*>  a>  e>  ™>  &>  i}> 
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Differentiating,  the  elements  being  considered  variable, 
we  have 

d(R]  _dRd(nt  +  e)     dR  da     dR  de     dR  d™ 
dt     ~  de        dt  da  dt      de  dt     d-sr  dt 

dRdti     dRdi 


and  differentiating  as  if  the  elements  were  invariable,  which 
is  permissible  for  the  reason  explained  in  Art.  26, 

d(R)        dR 
dt        U  de  ' 

Equating  the  two  values  of      ,     , 


dR  =  dRf        .dn     de\ 
n  de  ~~deT  +    (Tt  +  dt) 


de\      dRda     dRde 

cHii  d±u     aJct  d/i 

*   TTO  ~Tj7~  ~t~  ~3j~"~  Tji  * 
ail  at       a*  cit 

Substituting  for  -y- ,    -=-    &c.,  their  values 
ttc      dt 

dR  /  cZn     ^e\      2na2  ^  c?,K     na  (1  -  e2)  dJS  a'.K 
=  -JT    *  -37  +  ^z    + •    7    -T-  + 


de  \   dt      dt  de  da  e         de   de 


-e*)  /dR     dR\  dR,     na*J(I-e2)  dR  dR 
\de      d-sj)  de  pe          de  dvr 


dRdR  na  dR  dR 


—  e*)  di  dm     /z  V(l  —  ^J  sin  i  di 
7  T-.  7  T-,         na  tan  - 


na 


_ 
YdB     dR\dR 

2\  1  ~j~   i    ~J~  I  ~jT  > 


fj,  V(l  -  e2)  sin  i  dti  di       p  V(l  -  0  V  de      dsrj  di 
C.  P.  T.  3 
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dR/dn      de\      2na*  dR  dR 
de  I    dt  +  dt    +    LU     de  da 


dRdR 


ILL  V(l  -  e2)  de    di  ' 

rJ  7? 
Dividing  every  term  by  -=—  ,  and  transposing,  we  obtain 


de  dn     Zna?  dR     na  V(l  -  e2)  f1 

=     t~~-{     —  - 


37.  Of  the  formulae  which  have  been  obtained  for  cal- 
culating the  elements  of  the  orbit,  that  of  the  preceding  Article 
is  the  only  one  which  contains  a  term  proportional  to  the 
time*.  It  may,  however,  be  replaced  by  one  in  which  no 
such  term  exists.  For,  let  f  denote  the  mean  longitude, 
then 

f  =nt -f  e; 

c  d^  .dn     de 

therefore  j7=w  +  ^^z  +  37» 

dt  dt      dt 

d  f*>     r   j*\     j.^n     de 
or    jjtf-M^^^jg-i 

Now  let  (  =  fndt  +  e, 

then  by  the  formula  of  the  last  Article 

de'         ZncfdR     na  V(l  -  e2)  f1  2  ,  <£R 

^  ~*      ~~~  "V(1"       'j 


*  The  reader,  if  acquainted  with  the  Lunar  Theory,  will  have  already  seen 
the  inconvenience  of  such  terms. 
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Since  in  the  elliptic  formulae  e  never  occurs  except  when 
coupled  with  nt,  in  the  expression  nt  +  e,  it  will  be  altogether 
eliminated  if  for  nt  +  e  we  write  fndt  +  e. 

Considered  as  replacing  the  element  e,  e'  is  called  the 
epoch*'.  Since  however  we  shall  never  have  occasion  to  em- 
ploy the  formula  of  the  preceding  Article,  the  accent  will  in 
future  be  omitted. 

fndt  is  termed  the  mean  motion  in  the  disturbed  orbit, 
and  is  denoted  by  f. 

38.     To  obtain  a  formula  for  calculating  f,  we  have 

d^  _dn 
di?~~di' 

and  differentiating  the  equation  nza3  =  p, 

dn      _.  9  „  da     « 
-T-  +  SnV-57-0: 

dt  dt 

dn        3n  da        3n2a  dR 


f 
inereiore 


or  ^  =  — 


77     ~~  ^    , 

dt  p     de 

3n*a  dR 


dtf  fi     de  ' 

39.     We  will  here  recapitulate  the  formulae  which  have 
been  obtained  for  calculating  the  elements  of  the  orbit : 

da  _  2na*  dR 

de  =  na(l-e*)dR     na^/^-e2) 
dt  ue        de  tie 


, 
...,     cZtsr     na</(l-<ndR  2 

ml  _  —  __  v  ^  _  £  _          I 

~  ~ 


*  It  may  be  noticed  that  if  n  were  constant,  e'  would  be  identical  with  e. 

3—2 
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<iv>    ^-^S  +  ^^f1-^1-^ 

an2    dR 

.  .       (Zfl  _  wa  cLR 

rfi      fj,»J(l—  e*)smi  di  ' 

di_  na         f  1    dR  t         i  (dR  ,  dR\\ 

~    •  7.^  ~r  tan  — 


We  have  also  (vii)  the  equation 


but  this  forms  no  new  relation,  since  it  has  been  deduced 
from  (i). 

40.  When  the  elements  have  been  calculated  by  means 
of  the  above  formulae,  the  position  of  the  planet  will  be  given 
by  the  equations 


Note.  Some  of  the  Author's  methods  are  similar  to  those  given  by 
Mr  M.  O'Brien  in  an  article  "On  certain  formula  in  Physical  Astronomy," 
published  in  the  Camb.  Math.  Journal,  1843,  Vol.  in.  pp.  249  —  259.  Com- 
pare sections  18,  29,  30,  8,  20,  and  31  of  this  work  with  sections  6,  7,  8,  10, 
11  of  the  article  referred  to.  [A.  F.] 


CHAPTER  III 

DEVELOPMENT   OF   THE   DISTURBING   FUNCTION. 

41.  IN  the  first  Chapter  we  have  obtained  equations  by 
means  of  which  R  may  be  expressed  in  terms  of  the  time 
and  the  elements  of  the  orbit ;  we  now  proceed  to  shew  how 
the  actual  development  may  be  effected  in  a  series  ascending 
by  powers  and  products  of  the  excentricities  and  the  tangents 
of  the  inclinations.     In  the  Planetary  Theory  these  are  ex- 
tremely  small,   and   the  series  will  converge  rapidly.     Ac- 
cordingly in  the  present  treatise  small  quantities  of  orders 
higher  than  the  second  will  be  neglected*. 

42.  If  we  recur  to  Art.  11,  it  will  be  seen  that,  consider- 
ing only  one  disturbing  planet, 

1 


cos  (^  _  0;  )  +  (z  _ 


+zz'~] 

- 


The  first  step  towards  the  required  development  will  be 
the  expansion  of  rlt  r/,  01}  #/,  z  and  /  in  terms  of  the  time 

*  We  may  remark  that  to  this  order  of  approximation  the  inclinations, 
their*  sines,  and  tangents  will  be  equal. 
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and  the  elements  of  the  orbit.  For  this  purpose  we  may 
employ  the  equations  which  have  already  been  obtained  in 
Art.  13,  viz.: 

tan  (6l  -  fl)  =  cos  i  tan  (0  -  Q), 
sin  X  =  sin  i  sin  (6  —  H), 
rl  =  r  cos  X,         z  =  r  sin  X, 

r=  a-ll  +  ^e*  -e  GQ$(nt  +  £-&)—-£:  e*c,QsZ(nt  +  e  —  is}—  ..\  , 


with  similar  equations  involving  the  co-ordinates  and  elements 
of  the  disturbing  planet. 


(i)     To  expand  rt.     We  have 

r  =  r  cos  X  =  r  (1  —  sin2  X 

**'£'•  -rfl-'rin'X-f.     ) 


=  r  |l  -  ^  sin*t  sin*  (^  -  H)  +. . . [ 
-  itan2t'sin2(0-O)+...i 


to  the  same  order  of  approximation, 


or  substituting  the  expansions  for  r  and  0, 

1 


i==a  Jl  +-e2— -  tan2i- 

H-  -  tan2  i  cos  2  (w*  +  e  -  H)  +  . . .  i 

=  a  (1  +  w),  suppose. 

Similarly,  r/  =  a'  (1  +  u'). 
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(ii)     To  expand  0t.     We  have 
tan  (0t  -  0)  =  tan  {(0,  -  O)  -  (0  -  II)} 


1+  tan  ^  -  ft)  tan  (0  -  Ii) 
(cos  i  -  1)  tan  (0  -  Q) 


-2sin2|tan(0-fl) 
1  -f  tan2  (0  -  O)  -  2  sin2 1  tan2  (0  -  fl) 


-sin2|sin2(0-n)  • 


=  -sm2^sin2(0-ft)-...; 
therefore      0,  -  0  =  -  sin2 |  sin  2  (0  -  flj  -  . . . 

* 

=  -  tan2*r  sin  2  (0  -  H)  -  ... , 

2J 

to  the  same  order  of  approximation ;    or,  substituting  the 
expansion  for  0, 


sin 


(nt  +  e  -  w)  4-  z  e2  sin  2  (w<  +  6  -  w) 


TC 

-  tan2  1  sin  2  (nt  +  e  -  fl)  -+'  ... 
e  +  #,  suppose. 
Similarly,  0/  =  w'«  +  e'  +  vr. 
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(iii)     To  expand  z.    We  have 

z  =  r  sin  X  =  r  sin  i  sin  (0  —  H) 

=  r  tan  i  sin  (0  —  H)  —  ... 

to  the  second  order  ;   or,  substituting  the  expansions  for  r 
and  0, 

z  =  a  (tan  i  sin  (nt  +  e  —  fl)  +  ...}, 

in  which  terms  of  the  second  order  will  not  be  required. 
A  similar  expression  may  be  found  for  /. 

43.  Having  obtained  the  expansions  of  rx,  rt',  0lt  0,', 
#,  /  we  must  now  substitute  them  in  the  expression  for  E. 
This  may  be  effected  as  follows. 

Let  E/  be  the  value  of  E  when  u,  u,  v,  v   are  severally 
zero  :  then,  writing  (j>  for  nt  +  e  —  (n'tf  -f  e')  we  have 
Rt  =  ra'  [{a2  +  a'2  -  2aa'  cos  $+(z-  /)2}  ~  *. 

-  (aaf  cos  0  +  ^')  (a"  +  s/2)~  f] 


=  m'  [(a2  +  a"2  -  2aa  cos  ^)"  -       cos 


-  m     ~  (a2  +  a/2  -  2aar  cos  0)"   (^  -  zj 

1      ,     3  a 


,2 


TVT  T^  T>       ,  /  .  J        f      f      ,  .    ,  ,. 

Now  E  =  E  t  +  -j-*-  au  +  j~,  au  +  -j-r'  (v-v) 
da  da  d<f>  v 


,,       ffR,     ,  ,,        ,. 
v  —  t  v  )  +  ,  /  7  .  aw  (v  —  v  ) 
da  d</> 
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44.     It  will  be  shewn  in  a  subsequent  Article  that 


can  be  expanded  in  a  series  of  the  form 


Assume  then 

(a2  +  a'2  -  2aa  cos  0)~^  =  -  (7  +  (7  cos  0  +  (72  cos  20+ 


os+,  cos 
22 

Thus 

B  =  «'  \l  C0  +  f(7,  -  4)  cos  *  +  C2  cos  2<£ 

(A          \         a  / 

+  w'att(i^  +  (^-4AcoS^+^c 
(2  da      \cw-aV  aa 

,  ,  ,  (1  ^(70      /d(7t     2a\  d(7, 

+  m  aV  \-  -  -•?  +    —  i  +  --3   cos  0  +  ^  cos  20  + 
[2  da       \da       a3/  «» 

-  m  (v  -  vr)  \(C,  -  4)  sin  0  +  2(72  sin  20  +  ...I 
l\.    •    •  /  ) 


2cos^ 


,  f  1  d*C0      f  d*C,       2  \ 
+  maauu  fe  j—  7^7+  [  T—  r^  +  —  , 
[Zdada      \dada      av 


-r-7 
dada 
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-  ro'aV  '(v  -  v)  (f^J  +  -^  sin  <f>  +  2  ~2  sin  2  A  +  .  .  .} 
'  " 


+  -     sn  <>  +    ~  sn       +  .  .  . 
(\da      a"/  da  J 


45.     By  Art.  42, 
w  =  ^  e2  —  -7  tan2  a  -  e  cos  (nt  +  e  —  «r)  —  ~  e2  cos  2  (n£  +  e  - 


+     tan2  -i'cos  2  (irf  +  c-  Q)  +  ... 
v  =  2e  sin  (n£  +  €  —  w)  +  7  e2  sin  2  (nt  +  e  —  -or) 

TD 

-  tan2-  sin  2  (nt  +  e  -  H)  +  ..., 

2 

£  =  a  (tan  i  sin  (w£  +  e  —  II)  +  .  .  .}, 
with  similar  expressions  for  u  ,  v',  z  . 
Hence  w2  —  e2  cos2  (nt  +  e  —  or)  +  .  .  . 

e2     e2 
=  jr  +  g-  cos  2  (rtf  4  e  -  «•)  +  .  .  ., 

(t;  -  v'Y  =  4e2  sin2  (nt  +  e  -  w)  +  4e"  sin2  (n'i  +  e'  -  w') 


—  See'  sin  (rz^  +  e  —  -or)  sin  (w'^  +  e'  —  OT')  +  .  . 
2  (e2  +  e'2)  -  2e2  cos  2  (nt  +  €  -  w)  -  2e/2  cos  2  (n  *  +  e  -  w') 
4ee'cos(<^—  CT+CT')  +-4ee'cos  {(w  +  ri)  t  +  e+  ef—  OT  —  'cr/}+... 
=  ee'  cos  (wi  +  e  —  «•)  cos  (n't  +  e'  —  w')  +  .  .  . 


=  -=-  cos  (d>  —  &  +  CT  ;  +  -^-  cos  \(n  +  n  )  t  +  e  +  e  --or  — 

—  >j 
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u  (v  —  v)  =  —  e2  sin  2  (nt  +  e  —  -cr) 


+  2ee'  cos  (rtf  +  e  —  OT)  sin  (n'£  +  e  —  «•')  +  .  .  . 
—  —  ez  sin  2  (nt  +  e  —  -or)  —  ee'  sin  (0  —  -or  +  w') 

+  ee'sin  {(w  +  w')  tf  +  e  +  e'-OT  —  'w'}  +  ...  , 
,N8     a2  tan2/     a/2tan2i'      a2  tan2/ 

-,2r)2=  -  g  -  +  --  £  ---  2  - 

-a      nZ  cos2 


+  aa'  tan  i  tan  i'  cos  {(w  +  w')  «  +  e  +  e'  -  12  -  f2'}  +... 
&c. 

46.  If  these  values  be  substituted  in  Art.  44,  it  will  be 
seen  that  cosines  will  be  multiplied  only  by  cosines,  and  sines 
by  sines.     Hence  the  series  will  consist  of  two  parts,  one  inde- 
pendent of  t  explicitly,  and  the  other  consisting  of  periodical 
terms  of  the  form 

Pcos{(pn±qri)t  +  Q], 

where  p  and  q  are  any  positive  integers  or  zero,  P  is  a  function 
of  the  mean  distances,  excentricities,  and  inclinations,  and  Q 
a  function  of  the  longitudes  of  perihelia,  nodes,  and  epochs. 
The  former  part  is  denoted  by  the  symbol  Fi  we  proceed  to 
determine  its  value  as  far.  as  the  second  order  of  small  quan- 
tities. 

47.  To  determine  that  part  of  E,  which  is  independent  of 
the  time  explicitly. 

If  those  terms  only  be  written  down  which  either  are,  or 
after  reduction  will  become,  independent  of  t,  we  have 


adCje* 

1       + 


0  tanV\ 

~  ~~  ~T 
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a2  d"  00  ea     a'2  d*  C0  ef*     D.  / 
4  da*  2  +  4  da'a  2       4  \ 


a?  tan2/     a*  tan  V 
2 
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o  \  -i «  J  cos  0  4ee  cos  [6  — 

2V  tt 


™      sin   ee/  sin  <  - 


sin  <  ee'  sin  <  —  -nr  + 


+  a'  (  ~T-T  4-  -73  )  si 

—  fl  +  H')  +  ...[. 


~  cos  ^f>  aa'  tan  i  tan  i'  cos 


Now  cos  ^>  cos  (<£  —  -cr  4-  OT')  and  sin  </>  sin  ((^  —  -cr  +  -cr')  con- 
tain the  term  J  cos  («r  -  tsr'),  cos  ^>  cos  (0  -  H  +  II')  contains 
the  term  J  cos  (O  —  II')  ;  hence 


+  T  aa'^  tan  ^  tan  T  cos  (li  -  Ox)  +  ...i  . 

We  shall  hereafter  be  able  to  simplify  this  expression. 

48.  We  have  seen  that  the  remaining  terms  of  R  are 
of  the  form  P  cos  {(pn  ±  qn)  t  +  Q]  :  if  then  values  of  p  and 
q  could  be  found  such  that  pn  ±  qn  =  0,  this  term,  being  in- 
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dependent  of  t  explicitly,  would  form  an  additional  term  in  F. 
No  instance  of  this,  however,  occurs  among  the  planets. 

49.  In  consequence  of  the  extreme  smallness  of  the  ex- 
centricities  and  inclinations   of  the  orbits  of  the  principal 
planets,  terms  in  R  of  orders  higher  than  the  second  may  in 
general  be  neglected  :  but  it  sometimes  happens,  as  in  the 
Lunar  Theory,  that  higher  terms  become  sensible   through 
the  process  of  integration.     This  we  shall  consider  in  a  sub- 
sequent Chapter,  but  the  following  proposition  has  an  impor- 
tant bearing  on  the  subject. 

50.  The  principal  part  of  the  coefficient  of  a  term  in  R 
of  the  form  P  cos  {(pn  —  qn')  t  +  QJ  is  of  the  order  p  ~  q. 

DEF.  By  the  principal  part  of  the  coefficient  is  meant 
that  part  of  P  which  is  of  lowest  dimensions  in  e,  e',  tan  i, 
tan  i'. 

If  we  return  to  the  expression  for  R  in  Art.  44,  it  will 
be  seen  that  in  order  to  obtain  the  general  term  it  will  be 
necessary  to  multiply  the  product  of  the  general  terms  of  the 
expansions  for  ua,  uf&,  vy,  v8,  £e,  z'^  by  cos  k$  or  sin  kef). 

Now  (i)  in  the  expansions  of  u,  u',  v,  v't  zt  z  the  follow- 
ing law  is  observed  to  hold : — The  number  which  multiplies 
nt  +  e  or  n't  +  e  in  the  argument  of  any  term  represents  the 
order  of  the  principal  part  of  the  coefficient  of  that  term. 

(2)  The  same  holds  good  in  any  power  of  u,  u',  v,  v', 
z,  or  z.  For  consider  a  term  Pcos  (pnt  +  q)  in  v?.  It  can 
only  have  arisen  in  the  following  ways;  partly  from  the 
multiplication  of  two  terms  in  u  of  which  the  arguments  are 
Int  +  X  and  mnt  +  //,,  where  I  +  m  —p ;  and  partly  from  such  as 
have  the  arguments  l'nt  +  \r  and  m'nt  +  p ,  where  I'  ~m=p. 
In  the  former  case  the  order  of  the  coefficient  will  be  I  +  m, 
which  equals  p,  in  the  latter  it  will  be  I'  +  m,  and  this  is 
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greater  than  p.     Hence  the  principal  part  of  the  coefficient 
of  a  term  P  cos  (put  +  q)  in  u*,  will  be  of  the  order  p. 

Since  then  the  law  holds  in  u*,  it  may  be  shewn  in  like 
manner  to  hold  in  the  product  of  u2  and  u,  i.e.  in  u3.  Thus 
it  may  be  proved  for  any  power  of  u.  In  like  manner  it 
may  be  shewn  to  hold  for  any  powers  of  u ' ,  v,  v',  z  or  z'. 

(3)  The  same  law  is  true  for  the  product  of  any  powers 
of  u,  v,  z;   and  likewise  for  the  product  of  any  powers  of 
u,  v',  z .     This  may  be  proved  by  a  method  similar  to  that 
of  (2). 

(4)  In  the  product  of  any  powers  of  w,  u,  v,  v',  z  and 
zf,  the  order  of  the  principal  part  of  the  coefficient  is  the 
arithmetical  sum  of  the  multipliers  of  nt  and  n't. 

For  let  us  consider  a  term  M cos  {(In  ±  I'ri)  t  +  N}.  Now 
this  must  evidently  have  arisen  from  the  multiplication  of 
L  cos  (Int  +  X)  with  Lf  cos  (I' n't  +  X'),  or  of  L  sin  (Int  +  X) 
with  L'  sin  (Int  -f  X'),  where  by  (3)  L  is  of  the  order  I  and 
L'  of  the  order  I'.  Hence  M  will  be  of  the  order  I  +  I'. 

Now  any  term  in  the  development  of  R  of  the  form 
P  cos  {( pn  —  qri)  t  +  Q]  must  have  arisen  partly  from  the 

f*OQ1 

multiplication  of  P1  .    k<f>,  or  as  it  may  be  written 


with  Pa  C°S  {[(p  -Qn-(q-k)  n'}  t  +  Q,}, 

sin 

and  partly  from  its  multiplication  with 

•^8     •      {[(P  +  &)  ^  —(#  +  &)  w/]  ^  +  Q3}> 

where  k  is  any  positive  integer  or  zero,  Pl  is  a  function  of  a 
and  a'  only,  and  P2,  Ps  are  functions  of  the  excentricities  and 
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inclinations,  such  that  the  orders  of  their  principal  parts  are 
given  by  law  (4).  Hence  the  order  of  the  principal  part  of 
P  will  be  equal  to  the  lesser  of  those  of  P2  and  P3. 

Now  the  order  of  the  principal  part  of  P2  will  be  the 
least  value  of  which  the  arithmetical  sum  of  p  ~  k  and  q  ~  k 
is  susceptible,  for  different  values  of  k. 

(i)     Suppose  k  intermediate  to  p  and  q\  then  this  sum 
=  p~k  +  q~k=p~q-, 

(ii)  Suppose  k  not  greater  than  the  smaller  of  p  and  q  ; 
then  this  sum  —p  -f  q  —  2k,  the  least  value  of  which  (by 
putting  k  equal  to  the  smaller  of  p  and  q)  =  p  ~  q  ; 

(iii)  Suppose  k  not  less  than  the  greater  of  p  and  q  ; 
then  this  sum  =  2k  —  p  —  q,  the  least  value  of  which  (by 
putting  k  equal  to  the  greater  of  p  and  g)  —p  ~  q. 

Thus  p  ~  q  is  the  order  of  the  principal  part  of  P2.  That 
of  Pa  will  be  the  least  value  of  which  p  +  k+q  +  k  is  sus- 
ceptible, i.  e.  p  +  q. 

Hence  it  appears  that  the  order  of  the  principal  part  of 
P  is  p  ~  q. 

51.  The  principal  part  of  the  coefficient  of  a  term  in  R, 
of  the  form  P  cos  {(pn  -f  qn')  t  +  Q}  is  of  the  order  p  +  q. 

This  term  arises  from  the  multiplication  of  such  terms  as 


with  P, 


and  P3       {[(p 

3  sm  IIA^ 
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and  as  in  the  last  Article,  the  order  of  the  principal  part 
of  P  will  be  equal  to  the  lesser  of  those  of  P2  and  P3. 

Now  the  order  of  the  principal  part  of  P2  will  be  the 
least  value  which  the  arithmetical  sum  of  p  ~  k  and  q  +  k 
can  assume,  for  different  values  of  k. 

(i)     Suppose  Ic  less  tban^>;  then  this  sum 
—p  —  k  +  q  +  k  =  p  +  q. 

(ii)     Suppose  k  not  less  than  p ;  then  this  sum 

=  k-p  +  q  +  k, 
the  least  value  of  which  (by  putting  k  equal  to  p)  —p  +  q. 

Similarly  it  may  be  shewn  that  p  +  q  will  be  the  order 
of  the  principal  part  of  P3. 

Hence  it  follows  that  p  +  q  will  be  the  order  of  the 
principal  part  of  P3. 

In  Art.  44  we  have  assumed  that  (a2  +  a'2  -  Zaa  cos  <f>)~8 
can  be  expanded  in  a  series  of  cosines  of  0  and  its  multiples, 
we  shall  now  give  a  proof  of  this  and  shew  how  the  coeffi- 
cients may  be  calculated. 

52.  To  shew  that  (a2  -f  a'2  -  2aa'  cos  <£)~'  can  be  expanded 
in  a  series  of  cosines  of  multiples  of  <f>. 

r 

Suppose  a  greater  than  a,  and  for  -  write  a;  then 

(a2  -f  a/2  -  2aa'  cos  <£)"'  =  of*  (1  +  a2  -  2a  cos  <£)-« 

=  a"2*  [1  +  a2  -  a  (e+^1  + 


r 

=  a  N 
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$v-i  , 

+  sze         + 


5(5  +  1)  («  +  2)      8,*V-1 

—   u.  6  *T"  •  •  • 


3 


5  (5+  1)      2.   -2*V-] 

^          ae 


g(*  +  l)(a  +  2) 

13 


1? 


5j5+J)     5  (5+1)  (5+2)  ]   //V     '  +  /^ 

"  1 2      fc  [3  -U  2         7 


where  the  coefficients  of  e  and  e   £          will  always  be 

equal.     Hence  we  may  write 

(a2  +  a2  —  2ad  cos  <j>)~9 
*=-A0  +  Aicus(f>  +  A2cos%j>  +  ...  +  Ak  cos^+  ..., 

where  A&  Alf  &c.,  are  functions  of  a  and  a'.  The  series 
which  they  represent  will  be  always  convergent  provided  a  is 
less  than  unity,  or  a  greater  than  af.  If  a  be  less  than  a',  we 
have  only  to  interchange  a  and  a'  in  the  above,  so  that  a.  will 
then  denote  the  ratio  of  a  to  a. 

c.  P.  T.  4 
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53.     To  calculate  C0  and  Cr 

In  the  preceding  Article,  let  s  =  -  ;   then 


1   1.3    3     1.3   1.3.5   5 

2'274a  +  274'27¥76&  H 

Unless  a  be  small,  these  series  will  converge  too  slowly 
to  be  practically  useful.  More  convergent  series  might  be 
obtained,  but  Ponte'coulant  remarks  (Systeme  du  Monde, 
Tome  in.  p.  81),  that  in  practice  it  is  more  convenient  to 
employ  elliptic  integrals  for  the  purpose,  in  the  manner  we 
proceed  to  explain.  We  have 


cos 


-f  ^  02  (cos  3<^>  +  cos  <£)  +  . 

Integrating  both  sides  of  these  equations  with  respect  to 
between  the  limits  0  and  STT,  we  obtain 

* 


7T 


_ 
" 


.  -f-  a2  —  2a  cos 
cos  <b  dd>  1  F*          cos 


These  integrals  may  be  reduced  to  the  standard  forms  of 
elliptic  functions  by  assuming 

sin  (6  -</>)=  a  sin  6 ,..(1), 

whence  tan  0  = -r-S- (2). 

CQS  </>  —  a 
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From  (1)       cos  (6  -  $)  (l  -  f  f)  =  a  cos  0; 
V        aoj 

therefore  &  =  cggff  -*)-«<**  g 

d#  cos  (0  -  <£) 

Now 
cos  (0  —  (f>)  —  a  cos  6=  cos  6  (cos  </>  —  a)  +  sin  0  sin  <£ 

+  Sin  *  sin      b      2 


sn 


sn 

=  V((cos  0  -  a)2  +  sin2  $,  by  (2), 
=  V(H-a2-2acos^>)   ..............  (3). 

Also  cos  (0  -  0)  =  V(l  -  a2  sin2  0)  .....................  (4). 

TT  dd>          //I  +  a2  —  2x  cos  d> 


Again,  from  equations  (3)  and  (4) 

V(l  +  a2  -  2x  cos  (/>)  =  V(l  -  a2  sin2  (9)  -  a  cos  0; 
therefore  1  +  a2  -  2a  cos  </>  =  1  -  a2  sin2  0 

+  a2  cos2  0  -  2a  cos  (9  V(l  -  a2  sin2  (9), 
2*  cos  </>  =  2a2  sin2  0  +  2a  cos  6  V(l  -  a2  sin2  0), 
or  cos  0  =  a  sin2  0  +  cos  0  \/(l  —  &2  sin2  0). 

Now  as  0  increases  from  0  up  to  2?r,  0  also  increases  from 
0  to  27r  ;  hence 


0     air  Jo 


V(l  +  a2  -  2a  cos 

i  r2-        do 


4—2 
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C  =  —  r  _       cos  0  d<fr 
1  ~  «7rJ0  V(l  +  a2  -  2a  cos  0) 

1    P«      a  sin2  Odd 


cos 


—  a  sm 


Hence  with  the  usual  notation  for  elliptic  integrals  (see 
Todhunter's  Integral  Calculus,  Art.  222), 


The  numerical  values  of  Flo.,-}  and  ^(a,  -J  may  be 
found  from  Legendre's  tables  of  elliptic  functions. 


54.     Given  Ck  and  Ck-l  to  obtain  Ck+1. 
We  have 


01 
z 

differentiating  with  respect  to  <£, 

aa'  sin  <£  (a2  +  a/2  —  2aa'  cos  <^>)~^  =  (7X  sin  <£  +  2  (72  sin  20  -h  ... 

+  kCk  sin  ^0  +  .  .  .  ; 

therefore  aa  sin  ^  f  -  C0  +  (7t  cos  $  4-  .  .  .  J 


sin  0+  202sin 
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equating  coefficients  of  sin  k(f>, 
\  aaf(Ck_l-Ch+l)=k( 


whence         0, 


2k     a?  +  a*          2k  -1 


, 

,..      -  -    —  —  0.  —  --.—  .,  C,  ,. 
2&  +  1     aa         *      2&  +  1     ^ 

55.     Given  Ck  cmcZ  Ck+1  ^o  obtain  Dk. 
As  in  the  last  Article,  we  have 
aa  sin  c£  (a2  +  a'2  —  2aar  cos  £)"*  =  C^  sin  ^>  +  2  (72  sin 


therefore  aa'  sin  (f>  (-  DQ  +  Dl  cos  0  +  D2  cos  20  +  .  .  .  J 

=  (71sin</>  +  2<72sin20  +  ...; 
equating  coefficients  of  sin  &<£, 


a'(A.-1-AJ  ..................  (1). 

writing  k+  1  for  ^r, 


Again, 

(a2  +  a/2-  2aa'  cos  0)~f  =  |  D0  +  Dx  cos  0  +  D2  cos  20  +  .  .  .  , 
and 

(a2  +  a"2  -  2aar  cos  0)  ~^  =  |  (70  +  (7,  cos  0  +  C9  cos  20  +  .  .  .  ; 

therefore  -  (70  +  (7,  cos  0  -f  .  .  . 

=  (a2  +  a/2  -  2aa/cos  0) 

equating  coefficients  of  cos  &</>, 

C,=  (a*  +  a'*)Dt-aa'(D^  +  DM}  .........  (3), 

writing  k  +  1  for  k, 
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Eliminating  Dk_1  between  (1)  and  (3), 

(2k  +  l)Ck  =  (a2  +  a"2}  Dk  -  2aa'Dk+1 (5). 

Eliminating  Dk+2  between  (2)  and  (4), 

(2k  +  1)  CM  =  -  (a2  +  a'2)  DM  +  2aa'Dt (6). 

Finally,  eliminating  Dk+l  between  (5)  and  (6), 

(2k  +  1)  {(a2  +  a'2)  Ck -  2aa'  Ck+1]  =  {(a2  +  a2)2  -  4aV2j  Dk, 

k         i  ,y2 '2\2  (\*~     ~          /        k  k+l)' 

56.     To  calculate  the  successive  differential  coefficients  of 
Ck  and  Dk  with  respect  to  a  and  a'. 

We  have 

2i 

+  Ck  cos  /t' 

differentiating  with  respect  to  a} 

—  (a  —  a' cos  </>)  (a2  +  a2  —  2aa  cos  </>)     =  ^ -~  +  - 


substituting  for  (a2  +  a/2  —  2aax  cos  £)~*  its  expression  in  series 


ldC0     dCl  dCk 

equating  coefficients  of  cos  k<f>, 
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By  giving  to  k  in  succession  the  values  1,  2,  3,  &c.,  those 
of  -y-1  ,    -j-*  ,  &c.  may  be   found,  the   right-hand   member 

being  calculated  by  the  formula  of  the  preceding  Article.    By 
equating  the  parts  independent  of  <£,  we  obtain 


da 
_,  - 


The  value  of    _,  -  may  be  found  by   differentiating   the 


expression  for  D    in  Art.  55,  and  substituting  for  —  =-*  and 

da 

-~  their  values  as  given  by  the  present  Article. 

The  successive  differential  coefficients  of  Ck  and  Dk  with 
respect  to  a  may  be  obtained  from  the  expressions  for  ~~ 

and    ,—  *  by  simple  differentiation  and  substitution. 

57.  We  might  determine  in  the  same  way  the  successive 
differential  coefficients  of  Ck  and  Dk  with  respect  to  a  ;  but 
when  those  with  respect  to  a  have  been  found,  the  former 
may  be  derived  from  them,  as  we  proceed  to  shew.  On  ex- 
amining the  expansion  of  (a2  -f  a'2  —  2aa'  cos  <£)~*  in  Art.  52, 
it  will  be  seen  that  Ak  is  a  homogeneous  function  of  a  and  of 
of  —2s  dimensions.  Hence  Gk  and  Dk  are  homogeneous 
functions  of  a  and  a,  the  former  of  —  1,  the  latter  of  —  3 

dimensions.     It  follows  that  —  ,—  ,  -~  *  will  be  homogeneous 

da      da 

functions  of  —  2  and  —  4  dimensions  respectively  ;  and  so  on. 
Now  by  a  known  property  of  such  functions 

,aCk         „ 


which  determines  -^-j  : 
da 
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da*          da  da'         "  da  ' 

.;.,,'    .     #ct 

which  determines  ^  —  ~  : 
dada 


~r~  7»  7  —  7  /  T~/  1 

da*         dada  da 


which  determines    7—  *  :  and  thus  all  the  differential  coeffi- 
aa 


cients  of  Ct  may  be  determined. 


In  like  manner  all  the  successive  differential  coefficients 
of  D  may  be  calculated. 

We  are  now  in  a  position  to  simplify  the  expression  for 
F.  We  have  (Art.  47) 


2 


cos  [«:  — 


+  7  ^'A  tan  ^'  tan  *'  cos  (ii  -  no  +  . 

.'4 

The  following  proposition  will  be  found  useful. 

d2C  d2C 

58.     To  shew  that  -—,0,=  -  Dx,  and  that        ,1/  =  -  Dr 
dada  dada 

We  have 

-CQ+  Clcos(f>+  (72cos  2<^>  +  ...  =  (a'  -I-  a'2-2aa'  cos  '^)"» 


0     dC,  dC2 

therefore          s  -i-2  +  -j-1  cos  0  +  -j-3  cos  2rf>  +  .  .  . 
2  da       da  da 

-a-a  cos  <)   a2  +  a'2  -  2aa7  cos  < 
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therefore 

1  CL  (A        d  L'1  i  f  z       "2     c\     ' 

~  -j — j^j  4-  -= — -,-1-,  cos  (f>  4- . . .  =  cos  6  (a  +  a   —  2aa  cos 

2  da  da      da  da '        r  r  ^ 


+  3  (a  —  a  cos  <£)  (a'  —  a  cos  0)  (a2  -f-  d*  -  2aa  cos 

=  cos  </>  (a2  +  a'2  —  2aa'  cos  0)~* 

+  3  [aaf  (1  +  cos2  0)  -  (V  +  a2)  cos  0}  (a2  +  a2  -  2aa'  cos  ^>) 
=  cos  <>  a2  +  a'2  —  2daf  cos  </>)~^ 

cos 


=  -  2  cos  0  (a2  +  a'2  -  2aa'  cos 

+  3aa'  sin2  <£  (a2  -!-  a'2  -  2aa  cos  </>)^. 

Now     (a2  +  a'2  -  2aar  cos  <£f  *  =  -  Z>0  +  Dl  cos  ^>  +  .  .  .  ; 

<L 

differentiating  with  respect  to  0 

Saa'  sin  </>  (a2  4  a/2  -  2aa  cos  ^ 

=  D1  sin  ^  +  2Z>2  sin  2<£  +  .  .  .  ; 

1  d*C0        tfC. 

therefore  -  -,  —  =2_,  +  -r—^-,  cos  </>+... 

2  aa  aa      aa  aa 

=  -  2  cos  <f>  (-  D0  +  Dx  cos  0  +  ... 
+  sin  </>  (Dj  sin  <£  +  2D2  sin  20  +  ...), 


whence,  equating  the  parts  independent  of  <£,  and  also  the 
coefficients  of  cos  < 


2  ' 


or 


da  da' 

and  j—J-' 

da 'da 
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59.     Since  —~  is  a  homogeneous  function  of  a  and  a'  of 
—  2  dimensions, 

*C         ffCt_  _    dcr.. 

Ct      7    o    T"  tfc    ~i 5-7  —  —  —   — y —  , 

aa  da  da  da 


therefore  a  -^  +  -  — ^  =  -  ~  -— ^ 

1 


a-    -i    i  ,dC0     azd2Cn     1      /_ 

Similarly,        fyg+^yf^*?* 

Hence  the  coefficients  of  e2  and  e'z  in  the  expression  for  F 
are  each  equal  to  ^  oa'  Da. 

Again,  since  Gl  is  a  homogeneous  function  of  a  and  a  of 
—  1  dimensions, 

dC,       ,dC. 

a-rA  +  a--r*  =  -Cl-> 
da          da 

hence  the  coefficient  of  ee''cos  («r  —  -cr') 


but  (Art.  55)          2kCk  =  aa  (Dk^  -  Dk+l) ; 
therefore,  making  Jc  =  1, 

hence  the  coefficient  of  ee  cos  (sr  —  w') 
Again,  (Art.  56) 
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therefore  azD0  +  a  ^  =  aa  Dt. 

da 

Similarly,  a2DQ  +  a  ^  =  aa  A. 

Hence  the  coefficients  of  tan2  i  and  tan2  if  are  each  equal 
to  —  -  aa'D*. 

o 
Finally,  the  expression  for  F  becomes 

F=  m  \^Q+l  aa  Dl  (e*  +  e2)  -  \  aa  D2  eef  cos  («•  -  w') 


—  -  aa'Dj  (tan2  i  +  tan2  i')  +  -r  aa  D^  tan  i  tan  i'  cos  (O  —  H') 
8  4 


CHAPTER  IV. 

SECULAR  VARIATIONS   OF  THE   ELEMENTS   OF  THE  ORBIT. 
STABILITY  OF  THE  PLANETARY   SYSTEM. 

60.  WE  have  seen  in  the  preceding  Chapter,  that  the 
disturbing  function,  when  developed,  consists  of  two  parts ; 
the  one  independent   of  the  time  explicitly,  the  other  in- 
volving it  under  a  periodical  form :   we  shall  consider  sepa- 
rately the  effects  of  these  two  parts.    In  the  present  Chapter 
our  attention  will  be  directed  to  the  first  or  non- periodical 
part  of  R,  which  we  have  denoted  by  F.     The  inequalities 
thus  produced  in  the  elements  of  the  orbit  are  termed  secular, 
in  consequence  of  their  very  slow  variation. 

61.  By  differentiating  the  expression  for  F  in  Art.  59, 
with  respect  to  the  elements,  we  obtain 

' ..?-*   '"/•'-    . 

dF     m'     /r.     ,  .    ,  ,. 

-j—  =  —  aa  Dzee  sm  (^  —  vr  ) , 
tttxr        4 

dF    m      f  -T.         tn      ,.„   ,        ,  , 

-j-  —  -y~  aa  D$  —  -r-  aa  D2e  cos  («•  —  OT  ), 

r  -r\  i       'i       •!   •    ff\      r\'\ 
Dt  tan  i  tan  i  sin  (H  — 11 ), 
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^=  _  ^'  aa'D,  tan  i  +  -.-  aa'D^  tan  i'  cos  (fl  -  fl'), 

ttl  4  4 

-T-  =  an  expression  similar  to  F. 

Substituting  these  in  the  formulae  of  Art.  39,  and  neg- 
lecting powers  of  e,  e,  tan  i,  tan  i'  higher  than  the  second, 
we  have 


=  .,s.  0_ 

ac          4/i 
tan          =  -  2?f  A  (tan  »  -  tan  f  cos  (O  - 


~  =  A  +  A^  (e2 -  tanS')  +  A2  (ez  -  tan2 i') 

+  Azee  cos  (-CT  —  OT')  +  J.4  tan  i  tan  i'  cos  (fl  —  O'), 

where  in  the  last  expression,  A,  Alt  &c.,  have  been  written 
to  denote  certain  functions  of  a  and  a. 

62.  To  calculate  approximately  the  secular  variations  of 
the  elements  of  a  planet's  orbit,  in  a  given  time. 

Let  a0,  e0,  vr0,  &c.,  be  the  values  of  the  elements  at  some 
given  epoch;  a0  +  Sa*,  e0+  Be,  -G70-t-  &JT,  &c.,  their  values  after 
an  interval  t:  then  8a,  Be,  Sur,  &c.,  are  the  required  variations. 
By  Maclaurin's  Theorem, 

*  It  will  be  shewn  in  Art.  64  that  Sa  is  always  zero. 
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which  may  be  carried  to  any  required  degree  of  accuracy,  but 
in  practice  the  first  two  terms  will  generally  be  sufficient. 

We  have  supposed  the  variations  of  the  elements  required 
at  a  time  t  after  the  epoch ;  if  they  be  required  at  a  time  t 
before  the  epoch,  we  have  only  to  change  the  sign  of  t  in  the 
above. 

We  may  remark  that  I-,-} ,  (-7-) ,  &c.  are  of  the  order 
\dtJ0'  \dtJ0' 

of  the  disturbing  force,  since  they  involve  the  first  power  of 

m  :  (-77)  ,  ( — 1-3  )  ,  &c.  are  of  the  second  order;  for,  since  the 
\dt  /0   \dt  ,/0 

expressions  for  ~^-,  -j   ,  &c.,  involve  elements,  their  differ- 
dt     at 

ential  coefficients  will  involve  the  differential  coefficients  of 
those  elements,  and  thus,  by  substitution,  m'2  will  be  intro- 
duced:  similarly  (-TT)  ?  ("33")  »  *^c>>  are  °^  ^ne  thivd  order, 

\at  /0    \at  JQ 

and  so  on. 

In  the  short  period  of  one  year  all  terms  after  the  first 
may  be  neglected,  so  that  putting  t  =  1,  we  have 

Be  =  (~},&c. 


Hence  the  coefficient  of  t  in  the  above  formulae  is  called 
the  annual  variation. 

63.     Since  the  elements  of  the  planetary  orbits  are  con- 
tinually changing,  it  will  be  interesting  to  shew   that   the 
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dimensions  of  these  orbits,  and  their  inclinations  to  the 
ecliptic,  nevertheless  fluctuate  between  very  narrow  limits. 
This  constitutes  what  is  termed  the  Stability  of  the  Plane- 
tary System  :  in  order  to  establish  it,  it  will  be  necessary 
to  prove  the  stability  (i)  of  the  mean  distances,  (ii)  of  the 
excentricities,  (iii)  of  the  inclinations. 

64.  To  prove  the  stability  of  the  mean  distances  of  the 
planets  from  the  Sun,  and  of  their  mean  motions. 

By  Art.  61  -j-  =  0,  so  that  a  is  constant.     Now  it  will  be 

shewn  in  a  subsequent  Chapter  (see  Art.  91),  that  to  the  first 
order  of  the  disturbing  force,  the  periodical  terms  of  R  can 
produce  only  periodical  variations;  consequently,  to  this 
order,  the  mean  distance  is  susceptible  of  no  permanent 
change*.  The  same  is  true  of  the  mean  motion  n,  since 

Jr 
it  =  — | ,  and  //,  does  not  alter.     We  are  hereby  assured  of  the 

cr 

impossibility  of  any  of  the  bodies  of  our  system  ever  leaving  it 
in  consequence  of  the  disturbances  it  may  experience  from 
the  other  bodies ;  and  this  secures  the  general  permanence  of 
the  whole,  by  keeping  the  mean  distances  and  periodic  times 
perpetually  fluctuating  between  certain  limits  (very  restricted 
ones)  which  they  can  never  exceed  or  fall  short  of. 

This  result  may  easily  be  extended  to  all  orders  of  the 
excentricities  and  inclinations  :  for  since  nt  +  e  always  occurs 
in  R  as  one  symbol,  e  cannot  occur  in  ^because  t  does  not, 

Al    A  dF  ,      da  . 

so  that  — ,  and  therefore  --=-  is  zero. 
lie  dt 

*  This  result  is  also  true  when  the  square  of  the  disturbing  force  is 
included:  for  the  demonstration  the  reader  is  referred  to  Pontecoulant's 
Sy steme  du  Monde,  Tome  i.  p.  395  (2nd  edit.). 
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65.     To  prove  the  stability  of  the  eccentricities  of  the  pla- 
netary orbits. 

We  will  first  consider  the  case  of  two  planets  only.     By 

Art.  61, 

de         m'na*a'  -^   ,  .    ,  /x 

3T      -^T^  "a  (---)• 

0.    .,    ,     de          mria"*a      ,     .     .    ,        . 
Similarly,  -rr  —  --  7  -  Dz  e  sin  (w  —  -or). 
CLu  fyjM 

Now  since  D2  is  the  coefficient  of  cos  2<£  in  the  develop- 
ment of 


an  expression  in  which  a  and  a   are  similarly  involved,  it 
follows  that 

£.'=•»,• 

TT  1^-      1      '  xl  -U  x'  1  m  m'        > 

Hence,  multiplying  the  above  equations  by  —  e,  —^-,  e  , 
respectively,  and  adding,  we  have 

m    de      m'    ,def 
-e-r.+  —r-,  e  -j-  =  0  ; 
na    at     na      at 

therefore,  since  a  experiences  no  secular  variation, 

™e'+^,«'2=a 

na         na 

A  similar  equation  holds  for  any  number  of  planets.     Re- 
placing for  convenience  —j  —  -  by  (a,  a'),  we  have 

-j-  =  —  m'na  («,  a)  &'  sin  (tzr  —  -sr') 


de 

-,-=  —  wwV  (ar,  a)  e  sin  (or'  — 

at 


—  m/rwa  (a,  a/x)  e"  sin  (^  —  -sr")  —  .. 

r'-OT) 

—  m"rid  (a'j  a")  e"  sin  («r'  —  «•")  —  .. 
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de" 

-j-  —  -  mri'a"  (a",  a)  e  sin  (&"-  -or) 

at 

-  mri'a"  (a",  a) e'  sin  (tsr" -  <*')  -  . . . 


Since  D2'  =  D2,  it  follows  that  (a,  a)  —  (a,  a):  hence  mul- 
tiplying  these   equations   by  -—  e,  -7—,  e,  &c.,  and  adding, 

1\0b          71 CL 


we  obtain 


m    de      m'    ,de       m"    t,de"  A 

-  e  -T:  +  -7-3  ,  e  -TT  +  -77-77  e  -Tr+...  =  0; 
na    at      na      at      n  a        at 

whence  by  integration 


f 
\ 


—       =  const. 
na 


Since  wa  =     /-  ,  this  may  be  written 


Now  observation  shews  that  all  the  planets  revolve  round 
the  Sun  in  the  same  direction,  so  that  the  mean  motions  n,  n, 
n",  &c.  are  of  uniform  sign.  Hence  all  the  terms  of  the  left- 
hand  member  of  the  above  equation  are  positive. 

We  learn  also  from  observation  that  the  excentricities  of 
the  planetary  orbits  are  at  present  very  small  indeed,  and 
in  the  case  of  the  Asteroids,  the  masses  are  very  small. 
Hence  the  constant  must  be  small.  Since,  then,  all  the  terms 
of  the  first  side  of  the  equation  are  positive,  and  their  sum 
always  equals  a  small  constant,  it  follows  that  every  term  is 
small  ;  and  therefore,  except  in  the  case  of  planets  of  small 
mass,  such  as  Mercury,  Mars,  Juno,  &c.,  that  the  excentricities 
must  remain  permanently  small. 

The  stability  of  the  excentricities,  however,  is  not  con- 
fined to  the  larger  planets  :  we  shall  give  another  proof  of 

C.  P.  T.  5 
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this  important  theorem  and  that  of  Art.  66  in  the  following 
Chapter. 

66.     To  prove  the  stability  of  the  inclinations  of  the  planes 
of  the  planetary  orbits. 


By  Art.  61      =  -D^  tan  i'  sin  (fl  -  fl'). 


„,.    .,    ,       di'  -^f.       .   .    ,_, 

Similarly,    -77  =  —  -  --  D/  tan  i  sm  (O 
dt         ^rfji 

As   in   the   preceding   Article,   it    may   be   shewn    that 
/  =  Dr     Hence,  multiplying  the  above  equations  by 


m  ,      .,    Vfi    ,       ., 
—  tan  i,  -7—,  tan  %  , 
na  na 


respectively,  and  adding,  we  have 

m  ,       .di      m'          .,  di' 
-  tan  i.-ji  +  -j-f  tan  i  -j  -  =  0, 
na          dt     na  dt 

or  to  the  same  order  of  approximation, 

m  ,       .d  (tan  i]      m'          .,  d  (tan  i') 
-  tan  i  —^  —  '-  4  -r-7  tan  *•  -^-5-  —  ^  =  (); 

rca  dt  na  dt 

^1  /•  Wi     ,  9    •  ^     ,  9   V 

therefore  —  tan%  +  -p-rtan  i  =  const., 

wa  ?ia 

or,  since  na  =  ./-  , 

m  *Ja  tan2  i  +  m'  */a'  tanV  =  (7. 

A  similar  equation  would  (as  in  the  case  of  the  excentri- 
cities)  be  true  for  any  number  of  planets.  Now  the  inclina- 
tions of  the  planetary  orbits  to  the  ecliptic  are  at  present  very 
small;  hence,  if  we  take  for  our  fixed  plane  of  reference  a 
plane  coinciding  with  the  present  position  of  the  ecliptic,  and 
except  the  case  of  planets  of  small  mass,  it  follows,  as  in 
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Art.  65,  that  their  inclinations  to  this  plane   must   always 
remain  very  small. 

67.  It  should  be  noticed  that  the  assumption  in  the  two 
preceding  Articles,  that  p  is  the  same  for  all  the  planets,  is 
equivalent  to  neglecting  the  square  of  the  disturbing  force  : 
for,  let  IJL  refer  to  the  planet  m  and  //  to  m;   then,  if  M 
denote   the  Sun's  mass,  /ji  =  M+m,  pf  —  M+m,  so  that  // 
differs  from  /z,  by  a  quantity  of  the  order  of  the  disturbing 
force  :    since,  then,  the  expressions  in  which  IJL  and  ///  occur 
are  themselves,  of  the  first  order,  it  follows  that  the  error 
introduced  by  supposing  //,  and  ///  equal   is  of  the   second 
order. 

68.  The  results  of  Arts.  65  and  66  may  also  be  obtained 
directly  from  the  equations  of  motion.     We  have  (Arts.  20 
and  15) 

dhdRdR     dR 


7  77T 

or,  replacing  R  by  F,  since  we  have  seen  that  -—  -  =  0, 


dt      dis ' 

Multiplying  both  sides  of  this  equation  by  m,  forming 
similar  equations  for  each  planet  of  the  system,  and  adding, 
we  have 

dF\ 
ra-r-    . 


Now  on  referring  to  the  expression  for  F  in  Art.  59,  since 

m 
OT  occurs  only  in  the  term    —  — -  aa'D2ee  cosfa  -  «•'),   it   is 

easily  seen  that  2  fm  ,— J  =  0;  hence  our  equation  becomes 


5—2 
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whence  by  integration 

2  (mh)  =  const.; 
or  since  A2  =  //a  (1  —  e2)  =  (M  +  m)  a  (1  —  e2), 

2   mV  l  +         +  ...      l  -      -!-  ...     =  const. 


Since  a  is  constant  as  regards  secular  variations,  and  our 
approximation  extends  only  to  the  second  order  of  the  excen- 
tricities,  we  have  to  the  first  order  of  the  disturbing  force 

2  (m  V«e2)  =  0, 
the  equation  of  Art.  65. 

Again,  by  referring  the  motion  to  the  fixed  plane  of  refer- 
ence, we  obtain 

d  ,,        ~      dR     dR     dR     dR 

dt  (k  cos  *}  =  at  -  *  +  s  +  35 

and  considering  the  whole   planetary  system,  we   have    as 
before 

^  f     d  /i.        -vl      v  f     fdF 

Z  \m  -j*  (h  cos  l)  r  =  2  sm    j 
v 


Now  on  referring  to  the  expression  for  Fy  it  is  easily  seen 
from  the  forms  under  which  OT  and  II  occur,  that 

dF  dF 


hence  our  equation  becomes 

2  jm  ^  (Aoos«)j  =0, 

whence  by  integration 

2  (mh  cos  i)  =  const.*, 

*  It  should  be  borne  in  mind  tbat  this  and  the  equation  S  (mh)  =  C,  which 
are  those  which  would  be  obtained  by  conservation  of  areas  were  it  lawful  to 
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or  2  \  mil  (1—2  sin2  - )  [•  =  const., 

I        \  35/J 

which,  from  above,  may  be  written 

2  f  7?zA  sin2  p  J  =  const., 

and  proceeding  as  before,  we  obtain  to  the  second  order  of 
the  excentricities  and  inclinations,  and  the  first  of  the  dis- 
turbing force 


(in  \fa  sin2  -  J  =  const., 


or,  since  to  this  order  of  approximation,  the  inclinations,  their 
sines  and  tangents  are  equal,  this  may  be  written 

2  (m  ija  tan2  i)  =  C, 
the  equation  of  Art.  66. 

We  may  remark  that  the  equation 
2  (mh  cos  i)  =  const. 

is  of  itself  sufficient  to  establish  the  stability  both  of  the 
excentricities  and  inclinations.  For,  proceeding  as  before, 
we  obtain 


which,  to  the  second  order  of  the  excentricities  and  inclina- 
tions and  the  first  of  the  disturbing  force,  gives 


or  to  the  same  order  of  approximation 

2  (m  fjae*)  +  2  (m  ^a  tan2  i)  =  C. 

Since  we  know  from   observation   that   all  the  planets 
revolve  round  the  Sun  in  the  same  direction,  all  the  radicals 

assume  the  principle,  are  true  only  as  regards  secular  variations.  To 
assume  their  actual  truth  would  be  to  neglect  terms  in  E  due  to  the  Sun's 
motion,  of  the  first  order  of  the  disturbing  force. 
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in  this  equation  must  be  taken  with  the  same  sign.  Also, 
since  the  excentricities  and  inclinations  are  at  present  very 
small,  the  constant  must  be  small.  Hence  it  follows,  as  in 
Arts.  65  and  66,  that  the  excentricities  and  inclinations  must 
always  remain  very  small. 

69.  From  the  preceding  analysis  we  draw  the  following 
remarkable  conclusion  :  The  fact  that  the  planets  revolve  about 
the  Sun  in   the   same  direction,  ensures  the  stability  of  the 
planetary  system.     The  converse  of  this  would  not  necessa- 
rily be  true,  as  we  shall  see  in  Art.  75  :  the  numerical  rela- 
tions of  the  dimensions  and  positions  of  the  orbits  of  the 
planets,  might  be  such  as  to  ensure  stability,  although  they 
revolved  in  opposite  directions.     But  the  above  is  indepen- 
dent of  particular  numerical  relations. 

70.  The  results  at  which  we  have  arrived  with  regard  to 
the  stability  of  the  planetary  system  are  of  especial  interest. 
In  consequence  of  the  changes  in  the  elements  it  might  have 
been   supposed  that  the  orbits  would  ultimately   undergo 
such  alterations  in  their  dimensions  as  to  bring  the  planets 
into  collision  or  hurry  them  into  boundless  space.     Or  even  if 
no  such  violent  catastrophe  occurred,  a  derangement  of  the 
seasons   might  seriously  have  interfered  with   the  physical 
comfort  of  man*.     But  our  analysis  shews,  (and  the  results 
are   confirmed  when  the  approximation  is  carried  further,) 
that   in   so  far  as  the  mutual  attractions  of  the  Sun  and 
planets  are  concerned,  the  dimensions  and  position   of  the 
orbits  will  for  ages  remain  nearly  the  same  as  they  are  at 
present,  i.  e.  nearly  circular  in  form,  and  but  little  inclined  to 
each  other,  thus  affording  a  beautiful  illustration  of  Gen.  viii. 
22 :    "  While    the   earth  remaineth,  seed-time  and  harvest, 
and  cold  and  heat,  and  summer  and  winter,  and   day  and 
night  shall  not  cease." 

*  See  Herschel's  Outlines  of  Astronomy. 


CHAPTER  V. 

SECULAR   VARIATIONS     OF     THE    ELEMENTS    CONTINUED. 
INTEGRATION    OF   THE   DIFFERENTIAL  EQUATIONS. 

71.  IN  Art.  62  we  have  given  a  method  of  calculating 
the  secular  variations  sufficiently  accurate  for  the  practical 
purposes  of  astronomy,  but  in  order  to  understand  their  real 
nature,  and  thus  to  examine  more  fully  into  the  important 
question  of  the   stability  of  the  excentricities  and  inclina- 
tions, it  is  necessary  to  proceed  to  the  actual  integration  of 
the  equations  of  Art.  61.     This  we  are  enabled  to  do  by  a 
method  due  to  Lagrange. 

72.  To  integrate  the  equations  for  the  eccentricity  and 
longitude  of  perihelion. 

We  have  (Art.  61)  for  the  planet  m 

de         m'ncfa'  ^  ,   .    ,  ,. 


dm 


with  similar  equations  for  the  planet  m. 


72  PLANETARY  THEORY. 

We  shall  be  able  to  reduce  these  to  a  system  of  linear 
differential  equations  if  we  assume 

u  =  e  sin  -or,          v  =  e  cos  CT, 
u'=  e'sin  -or',          v'  =  e  cos  -nr'  ; 

.,       ,,  du  dm       .       de 

therefore  -y-  =  e  cos  -or  -T-  -f  sin  w  -y  . 

iff  cfa  eft 

Substituting  the  values  of  -3-  and  -y-  ,   and  writing  a 


„  , 

for    —  -  -  ,  we  have 


du        .^  f.  ,          , 

--  =  a.  \P^  cos  w  —  l>2e  cos  OT 


Similarly,  ^T~a  (AX  ~~  DIU)- 

In  like  manner  for  the  planet  m',  writing  a!  for 
we  have 


mn'a'^a 


;      The  forms  of  these  equations  suggest  the  following  par- 
ticular integrals  : 

u  =  M  sin  (^  +  7),         v  =M  cos  (gt  +  <y), 
u'=M'sin  (gt  +  7),         v'  =  IT"  cos  (#£  +  7). 

Substituting  these  in  the  differential  equations,  we  obtain 
from  either  of  the  first  two 


and  from  either  of  the  last  two 
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eliminating  the  ratio  M  :  Mr 


or  g"  -  (a  +  a')  D,g  +  act'  (Dt2  -  Dfi  =  0  : 
and  the  roots  of  this  equation  will  be  real  and  unequal,  real 
and  equal,  or  impossible,  according  as 

(a  +  ayZ^-te'^-A2) 
is  positive,  zero,  or  negative.     Now 

(a  +  a)2  A2  ~  4aa'(A2  r  *>/)  =  («  ~  a7  A"  +  **a.'Vf, 
a  positive  quantity,  since  n,  n  and  therefore  a,  a7  are  of  like 
sign.  Hence  the  values  of  g  will  be  real  and  unequal  :  denote 
them  by  glt  ga,  and  let  7l,  7a;  Ml9  Mj  M^  M^\  be  the  cor- 
responding values  of  7,  M,  M'  respectively.  Then  the  com- 
plete solution  of  the  differential  equations  will  be 

u  =  M1  sin  (gf  +  7j)  +  M2  sin  (gf  +  72)> 
w  =  Ml  cos  (^  +  7l)  +  M9  cos  (^  +  7a), 
^  =  J//  sin  ^  +  7l)  +  Jl//  sin  (gjt  +  7a), 
vr  =  Jf/cos  (^  +  7l)  +  Mf  cos  (£,*+  7a). 

Of  the  constants  in  these  equations,  four  are  arbitrary  and 
must  be  determined  from  observation.     We  have 


tan  cr  =  -  = 


cos  {(ffl  -^)  *  +  7l  - 
sin  (^i^  +  %)  +  ^2  sin  (^  +  7,)  . 


v     Ml  cos  (gj  +  7l)  +  J/2  cos  (gj  +  7f)  ' 
with  similar  equations  for  ex  and  «•'. 

73.     Had  we  considered  a  system  of  several  planets,  we 
should  have  obtained  by  a  like  process 


cos        - 
sin 


tan  ^  = 

i  c°s  ft*+7i  +2  cos 

with  similar  equations  for  each  of  the  other  planets. 
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74.  We  may  hence  infer  the  stability  of  the  excen- 
tricities.  From  the  form  of  the  expression  for  e  it  appears 
that  ez  cannot  be  greater  than 


... 
and  therefore  that  e  cannot  exceed 


these  quantities  being  all  taken  with  the  same  sign. 

Thus,  by  determining  the  numerical  values  of  l/t,  l/2, 
&e.,  for  any  particular  planet,  we  may  assign  an  actual  limit 
which  its  excentricity  can  never  exceed.  For  the  principal 
planets,  Mlt  l/2,  &c.,  are  found  to  be  exceedingly  small,  so 
that  the  excentricity  must  always  remain  very  small. 

If  we  confine  our  attention  to  two  planets,  we  see 
from  the  expression  for  e  in  Art.  72,  that  the  excentricity 
fluctuates  between  the  limits  l/x  +  J/2  and  l/x  ~  J/2.  The 

period  of  these  changes  =  -  ,  and  is  the  same  for  each 

0*7& 
planet  :  it  appears  also  from  the  equation 

m  tjae2  +  m  >Jafe2  =  C, 

obtained  in  Art.  65,  that  the  maximum  of  each  excentricity 
takes  place  at  the  time  of  the  minimum  of  the  other. 

As  an  illustration,  take  the  case  of  Jupiter  and  Saturn. 
Sir  John  Herschel  finds  that 


<7t  =  21"-9905,     g9= 
Jfj  =  -  -01715,  l/2  =  -04321,  for  Jupiter  ; 
Jf/=     -04877,  M^=  -03532,  for  Saturn: 

the  year  1700  being  taken  as  the  epoch*.  Thus  we  obtain 
for  the  greatest  and  least  excentricities  that  Jupiter's  orbit 
can  attain,  '06036  and  '02606  respectively,  and  for  those  of 

*  Article  Physical  Astronomy  in  the  Encyclopaedia  Metropolitana. 
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Saturn,  '08409  and  *01345,  quantities  exceedingly  minute  ; 
while  the  period  of  these  changes  amounts  to  no  less  than 
70414  years. 

75.  It  appears  from  the  preceding  Article  that  the  sta- 
bility of  the  excentricities  is  a  consequence  of  the  periodical 
form  of  the  solution  of  the  differential  equations,  a  result 
which  depends  upon  the  fact  that  gl  and  gz  are  real  and 
unequal.  Now  we  have  seen  that  in  order  that  this  may  be 
the  case,  it  is  only  necessary  that 


shall  be  positive,  a  condition  which  might  be  satisfied  if  the 
signs  of  ny  n,  and  therefore  of  a,  a'  were  different.  In  this 
case,  then,  the  stability  would  still  subsist.  Let  us,  however, 
consider  what  would  be  the  effect  of  equal  or  impossible  roots 
to  the  quadratic  from  which  g  is  found.  In  the  former  case 
a  term  would  be  introduced  into  u}  u',  v,  and  v'  proportional 
to  the  time,  and  in  the  latter  the  periodical  terms  would  be 
replaced  by  exponentials.  Consequently  the  excentricities 
would  increase  indefinitely  with  the  time,  and  the  stability 
would  no  longer  subsist. 

76.  We  now  proceed  to  examine  the  expression  which 
has  been  obtained  in  Art.  72,  for  the  longitude  of  perihe- 
lion, viz. 


tan  v  =      *  sn        +       4       sn 


M1  cos  (gf  +  7,)+^  cos(g2t  +?2)  ' 


dt  -        M'+  M*  +  2MM  cos        - 


The  maxima  and  minima  values  of  -or,  if  such  exist,  will  be 

found  by  equating  -j-  to  zero.     Thus 
dt 
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If  this  (disregarding  sign)  be  not  greater  than  unity,  the 
perihelion  will  oscillate,  the  period  of  a  complete  oscillation 

being  the  same  as  that  of  the  excentricities,  viz.  -     — ;  but 

y i    i/2 
if,  as  is  the  case  with  Jupiter  and  Saturn,  this  be  greater 

than  unity,  the  longitude  of  perihelion  has  no  maximum  or 
minimum,  and  the  perihelion  moves  constantly  in  one  direc- 
tion. 

Again, 

--  ... 

W 


dt      2J^  +  J/2*+2^2cos  {(gi-gz)t  +  7l-ry2}      2 
-  Mf)     1 , 

-~+v(9i+y*)- 


Hence  when  e  is  a  maximum  or  minimum,  -=-  will  be 

ctt 

either  a  maximum  or  minimum,  and  the  apsidal  line  will  be 
moving  most  rapidly  or  most  slowly,  different  cases  occur- 
ring according  to  the  signs  and  magnitudes  of  the  quantities 
involved. 


77.     When  the  apsidal  line  oscillates,  to  find  the  extent 
and  periods  of  its  oscillations. 

We  have  (Art.  72) 
tan  OT  =  M  sin 


Ml  cos  (gj  +  7,)  +  M2  cos  (gf  +  yj  ' 

**•'./  N       tan  -sr  —  tan  (at  +  7l) 

therefore     tan  («r  —  a  ,t  —  7l)  =  -  —  ^  --  — 

Jl       7l/      1  -h  tan  «•  tan  ( 

—  If  sin  Ar 


COS 
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Also  by  the  last  Article,  if  r  be  the  least  positive  angle 
whose  cosine  is  .  --1—  1  —      - 


COS  T' 


Different  cases  will  occur  according  to  the  signs  of  Mlt 
Mz,  &c.  Suppose  Mlt  M2  of  like  sign,  gl  and  gz  positive,  and 

gl  greater  than  g2.     Then  ^  increases  as  t  increases,  and  ~^~ 

will  be   negative,   or   the   apsidal    line  will   regrede,  while 
•cos  ^r  +  cos  T   is   negative,   i/e.    so   long    as   ^    is    between 

(2n  —  1)  TT—  T  and  (2w  —  1)  TT  +  T:  —j-  will  be  positive,  or  the 

apsidal  line  will  progrede,  while  ^  is  between  (2»  —  1)  tr  +  r 
and  (2n  +  1)  TT  -  T. 

To  find  the  angle  through  which  the  apsidal  line  regredes 
and  the  period  of  the  regression.  Let  t't  t"  be  the  values 
of  t,  -or',  is"  the  values  of  OT  corresponding  to  the  values 
(2n  -  1)  TT  -  T  and  (2n  -  1)  TT  +  r  of  ^  :  then 


-      =     ?i  ~ 


tan  (^-^/-7i)=— -^ 

MI  —  Mz  cos  T 


/   //        ."        \ 
tan(w  -^  -% 


—          COS  T 


From  these  equations  the  values  of  t\  t",  -car',  -sr'7,  may  be 
found,  and  thus  «/  —  «•"  the  amount  of  regression  will  be 
known,  The  period  of  regression 

2r 
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In  like  manner  the  amount  and  period  of  the  progression 
may  be  obtained.  The  latter  will  be  found  to  be  -  — - . 

The  period  of  a  complete  oscillation  will  be  the  sum  of 
the  periods  of  the  regression  and  progression,  that  is  -  -  , 
which  agrees  with  the  preceding  Article. 

78.  The  motion  of  the  centre  of  the  instantaneous  ellipse 
in  consequence  of  the  secular  variations  of  e  and  OT  may  be 
exhibited  geometrically  as  follows. 

We  have,  by  Art.  72, 

e  cos  -57  =  M^  cos  (gf  +  yx)  +  M2  cos  (g£  +  y2), 
e  sin  to-  —  MI  sin  (gf  +  y4)  +  J/2  sin  (gf  +  ya). 

Let  a  circle  be  described  in  the  plane  of  the  orbit  with 
its  centre  8  coinciding  with  that  of  the  Sun,  and  its  radius 
equal  to  M^a,  where  a  is  the  mean  distance.  Let  a  point  P 
describe  this  circle  uniformly  with  a  velocity  glf  starting  from 
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0.  Again,  with  centre  P  and  radius  equal  to  M2a  let  another 
circle  be  described,  and  let  a  point  Q  describe  this  circle 
uniformly  with  a  velocity  <72,  starting  from  C.  Let  SL.  be 
the  line  from  which  longitudes  are  reckoned,  and  draw  PK 
parallel  to  it:  then  if  the  angle  OSN  be  equal  to  7^  and 
CPK  to  72,  the  angle  PSN  will  be  equal  to  gf  +  ylt  and 
QPK  to  gf  4-  72.  Produce  QP  to  meet  the  circle  again  in  R, 
and  draw  QN  perpendicular  to  SN.  Then,  supposing  J/x  and 
3/2  to  be  both  positive,  we  have 

SN  =  SP  cos  PSN  +  PQ  cos  QPK 

=  Mja  cos  (gj  +  7,)  +  M2a  cos  (gj  +  72) 

=  ae  cos  •or. 

Similarly,  it  may  be  shewn  that 
QN  =  ae  sin  OT. 

Hence,  the  apse  being  supposed  to  move  from  L  in  the 
direction  contrary  to  that  of  the  hands  of  a  watch,  Q  will  be 
the  centre  of  the  instantaneous  ellipse. 

If  Ml  be  positive  and  M2  negative,  it  may  be  shewn  in 
like  manner  that  the  centre  of  the  ellipse  will  be  R.  If 
Mt,  Mz  be  both  negative,  join  QS  and  produce  it  to  Q'  so 
that  SQ  =  8Q  :  then  the  centre  of  the  ellipse  will  be  Q'. 

A  similar  construction  will  of  course  apply  for  the  motion 
of  the  further  focus. 

It  is  easily  seen  from  the  above  that  the  excentricity  is 
least  when  Q  is  in  the  line  SP  and  greatest  when  Q  is  in  the 
line  SP  produced.  Hence  the  maximum  and  minimum 
values  of  the  excentricity  are  M^  +  Mz  and  Mv  ~  J/2  respec- 
tively. (Art.  74.) 
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79.     To  integrate  the  equations  for  the  inclination  and  lon- 
gitude of  the  node. 

We  have  (Art.  61)  for  the  planet  m 

i  __  mnaa  ^  ^  ^.,  g.Q  ,^  _  ^ 

dt  4//, 

.d£l         m'ncfa   -^  (l       .  .,        /^      ^/^ 

tan  i  -T-  = x>j  {tan  t  -  tan  i  cos  (il  -  fl )} ; 

with  similar  equations  for  the  planet  m. 

To  integrate  these,  assume 

p  —  tan  i  sin  O,      g  =  tan  i  cos  H, 
£/  =  tan  ^  sin  IV,   gr  =  tan  i'  cos  IT ; 

therefore  -~  =  tan  i  cos  II  -j-  +  sin  H  (1  -f  tan2^)  ^- . 

Substituting  the  expressions  for  -y-  and  -7, ,  and  writing 

at          at 

m'na?a'      .  « .  di ,    .         ,,  , ,     , ,  .   -,       -,  , 

a  for since  tan^  -y-  being  of  the  third  order  may  be 

^fju  at 

omitted,  we  have 

3-  =  a^  (tan  i'  cos  H'  —  tan  i  cos  H) 


Similarly,  ^  =  aDj  (p  —  p). 

Also  for  the  planet  m,  writing  a'  for 


4/4 
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The  forms  of  these  equations  suggest  the  following  par- 
ticular integrals  : 

p  =  JV  sin  (ht  +  8),         q  =  Ncos  (ht  +  8), 
8),        2'=JV"cos  (ht+  8). 


Substituting  these  in  the  differential  equations,  we  obtain 
from  either  of  the  first  two, 


and  from  either  of  the  last  two, 
hN'  =  oL'DL( 

eliminating  the  ratio  N  :  N't 


or  +  (<x  +  ax   =  0; 

therefore  h  =  —  (a  +  a')  Dv  or  h  =  0. 

Denote  the  former  by  \t  and  let  8t,  82,  Nlf  Nz,  J\T',  JV/, 
be  the  values  of  8,  ^V,  N'  corresponding  to  h  =  h1  and  h  =  0. 
Then  j\y  =  AT2,  and  the  complete  solution  of  the  differential 
equations  will  be 

p  =  Ni  sin  (^e  +  St)  +  JV2  sin  88> 
gr  =  JV;  cos  (\t  +  Sx)  +  ^2  cos  82, 
jp'=  J\V  sin  (Kf  +  aj  +  iY2  sin  S2, 
5'  =  JV/  cos  (hf  +  8,)  +  N2  cos  S2. 

Of  the  constants  in  these  equations,  four  are  arbitrary, 
and  must  be  determined  from  the  known  values  of  i  and  O  at 
some  given  epoch. 

We  have  then 
tan2  i  =p*  +  (f  =  N?  +  N;  +  2N^2  cos  (hf  +  Bl  -  Ba\ 

n  o  -ff  _  ^  Bin  (V  +  8,)  +  ^2  sin  S2 
?     JVt  cos  (^+  SJ  +  j^cosV 

with  similar  equations  for  i'  and  H'. 

c.  P.  T.  6 
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Had  we  considered  a  system  of  several  planets,  we  should 
have  obtained  a  result  precisely  similar  to  that  of  Art.  73. 

80.  From  the  form  of  the  expression  for  tan  i,  the  sta- 
bility of  the  inclinations  may  be  inferred.  For  it  may  be 
shewn,  as  in  Art.  74,  that  tan  i  can  never  exceed 


these  quantities  being  taken  with  the  same  sign.  Since 
then  JVj,  N2,  JV3  &c.  are  found  to  be  very  small,  it  follows 
that  the  inclinations  must  always  remain  exceedingly  small. 

In  the  case  of  two  mutually  disturbing  planets,  we  learn 
from  the  expression  in  Art.  79,  that  tan  i  fluctuates  between 
the  limits  N^  +  N~z  and  Nt~  Nz.  The  periods  of  these 

2?r 
changes  are  the  same  for  the  two  planets,  being  —7-  ;  and 

as  appears  from  the  equation  of  Art.  66,  the  maximum  of 
each  inclination  will  take  place  at  the  time  of  the  minimum 
of  the  other. 

In  the  case  of  Jupiter  and  Saturn,  the  period  is  50673 
years  ;  the  maximum  and  minimum  inclinations  of  Jupiter's 
orbit  to  the  ecliptic  are  2°  2'  30"  and  1°  17'  10",  those  of 
Saturn's  orbit  2°  32'  40"  and  0°47'. 

81.  We  now  proceed  to  examine  the  expression  which 
has  been  obtained  in  Art.  79  for  the  longitude  of  the  node. 
We  have 


-ZVj  cos  (hj  +  SJ  +  Nz  cos  S2 " 
The  maxima  and  minima  values  of  O,  if  such  exist,  will 
be  found  by  equating  -=-  to  zero.     Thus 

CLu 

cos  (h.t  +  8,  -  M  =  -  -TF  . 
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If  this  (disregarding  sign)  be  not  greater  than  unity,  the 
node  will  oscillate,  the  period  of  a  complete  oscillation  being 

the  same  as  that  of  the  inclinations,  viz.  —  -j-  .     But  if  it  be 

±fit 

greater  than  unity,  there  cannot  be  any  stationary  positions, 
and  the  node  will  move  continually  in  one  direction. 

It  may  be  shewn,  as  in  Art.  76,  that  the  motion  of  the 
node  will  be  fastest  or  slowest  whenever  the  inclination  is 
either  a  maximum  or  a  minimum. 

82.  When  the  line  of  nodes  oscillates,  to  find  the  extent 
and  periods  of  its  oscillations. 

It  may  be  shewn  as  in  Art.  77,  that  if  ^  be  written  for 
h1t+S1-B2,  and  T  denote  the  least  positive  angle  whose 

.       .    ^ 
cosine  is  W, 

IV2 

-^  sin 


cos  T  sn 


~~  1  +  COS  T  COS  i/r  ' 

and  tan2  i  -^-  =  h^N^  (cos  ty  +  cos  T). 

Different  cases  will  occur  according  to  the  signs  of  NI} 
Nz  and  1\.  Suppose  Nlt  N2  of  like  sign,  and  \  positive: 
then  ^  increases  as  t.  increases,  and  the  line  of  nodes  re- 
gredes  so  long  as  -v/r  is  between  (2n—  1)  TT  —  T  and  (2n  —  1)  TT+  r, 
and  progredes  so  long  as  ^  is  between  (2n  —  1)  TT  +  T  and 
(2n  +  1)  TT  -  T, 

Let  O',  fl"  be  the  values  of  O  corresponding  to  the  values 
(2n  —  1)  TT  -  T  and  (2n  —  1)  TT  +  T  of  ^r  ;  then 

tan  (H7  —  S2)  =  cot  T, 

tan(rr-S2)=-cotr; 

6—2 
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therefore  Ii'  —  8  =  mir  +  —  —  r, 


therefore  fl1  -  fl"  =  w  -  2r, 

which  is  the  angle  through  which  the  line  of  nodes  regredes. 
Also  the  period  of  this  regression  may  be  shewn  as  in  Art.  77 

to  be  j—.     Similarly,  the  angle  through  which  the  line  of 

"i 
nodes  progredes  may  be  shewn  to  be  TT  —  2r,  and  the  period 

UL  '          2(7T-T) 

of  the  progression  -  -—  -,  -  -  . 
"'i 

The  period  of  a  complete  oscillation  will  be  the  sum  of 

the  periods  of  the  regression  and  progression,  that  is  -y-  , 

"i 
which  agrees  with  the  preceding  Article. 

The  remaining  cases  corresponding  to  different  arrange- 
ments of  the  signs  of  Nlt  N2  and  \  may  be  treated  in  like 
manner. 

The  mean  value  of  O  is  mir  -f  82,  whatever  be  the  signs 
of  Nv  Nz  and  hiy  and  the  mean  value  coincides  with  the 
true  whenever  sin  ty  =  0.  Since,  then,  ^  is  the  same  both 
for  the  disturbed  and  disturbing  planet,  the  nodes  of  both 
orbits  will  arrive  simultaneously  at  their  mean  positions. 

In  the  case  of  Jupiter  and  Saturn,  N2  is  for  each  planet 
numerically  less  than  NI}  so  that  the  node  oscillates  ;  the 
extent  of  oscillation  being  13°  9'  40".  in  Jupiter's  orbit,  and 
31°  56'  20"  in  that  of  Saturn  on  either  side  of  their  rnean 
position,  the  ecliptic  being  taken  for  the  plane  of  reference, 
and  supposed  immoveable. 
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83.  To  shew  that  the  inclination  of  the  orbits  of  two 
mutually  disturbing  planets  to  each  other  is  approximately 
constant. 

If  7  denote  this  inclination,  we  have  by  Spherical  Trigo- 
nometry, 

cos  7  =  cos  i  cos  i'  +  sin  i  sin  i'  cos  (fl  —  O') 

=  cos  i  cos  i'  (1  +  tan  i  tan  i'  cos  (H  —  O')} 

=  (1  +  tanV)-*  (1  +  tanY)-*  {1  +  tan  i  tan  z"  cos  (Q,  -  IT)} 

=  l-o  {tanV  +  tanY  -  2  tan  i  tan  i"  cos  (£1  -  H')J, 

if  we  neglect  small  quantities  of  orders   higher  than  the 
second. 

Now     tansi  +  tanY  —  2  tan  i  tan  i'  cos  (H  —  IT) 


therefore  1  —  cos  7  —  «  (^i  "" 


sm    =     (~.;);   ;     _  . 

whence  it  follows  that  7  is  constant. 

84.  The  equations  which  give  the  secular  variations 
of  the  node  and  inclination  may  be  explained  geometrically 
as  follows*. 

The  equations  to  be  interpreted  are 

p**Nt  sin  (\t  +  SJ  +  N2  sin  S2, 
q  =  #i  cos  (\t  +  Sx)  +  N2  cos  S2, 
where  p  =  tan  i  sin  H,     ^  =  tan  z  cos  H. 

*  This  elegant  geometrical  explanation  is  due  to  Mr  H.  M.  Taylor,  M.A., 
Fellow  and  Tutor  of  Trinity  College,  Cambridge.  Oxford,  Cambridge  and 
Dublin  Messenger  of  Mathematics,  Vol.  in.  p.  189. 
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Let  NA,  NB  be  the  intersections  of  the  planes  of  refer-' 
ence  and  of  the  orbit  respectively  with  a  sphere  of  radius 
unity,  the  centre  of  which  coincides  with  that  of  the  Sun  : 
ABZP  the  great  circle  of  which  the  pole  is  N:  Z  and.  Pthe 
poles  of  the  great  circles  NA,  NB :  L  the  point  from  which 
longitude  is  measured. 

Then  LN=  O,  and  PZ=BA  =  i. 


Now  project  PZ  and  the  other  great  circles  by  radii 
drawn  from  the  centre  of  the  sphere  on  the  tangent  plane  at 
Z:  then  tan/ sin  H  and  tan /cos  II  are  the  Cartesian  co- 
ordinates of  the  projection  of  P  referred  to  the  projection 
of  LZ  as  axis  of  y,  and  a  line  at  right  angles  to  it  as  axis 
of  oc. 

Suppose  P'  the  projection  of  P,  then  if  the  co-ordinates 
of  Pr  be  x  and  y,  we  have 

x  =  Nt  cos  (hf  +  8,)  +  N2  cos  S2> 
y  =  Nj_  sin  (hj  +  SJ  +  N2  sin  82. 

These  equations  shew  us  that  P'  always  lies  on  a  circle  of 
which  the  centre  is  at  the  fixed  point  (JV"2cosS8,  j^asin88), 
and,  the  radius  is  N^.  also  that  P  describes  this  circle  with  a 
uniform  angular  velocity  ht. 
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We  may  hence   arrive   geometrically  at  the   results  of 
Arts.  81  and  82.     If  C  be  the  centre  of  the  circle,  we  have 


angle  CZx  =  S2 ,     angle  P'CO  =  l\t  +  8,  -  S2 ; 
and  if  P'Z  be  joined, 

P'Z=it  angle  P'Zx  =  O. 

Suppose  NI}  N2  both  positive;  and  first,  let  iYx  be  less 
than  Nt2. 

Now  the  only  time  when  the  node  will  be  stationary  will 
be  when  P'  is  moving  directly  towards  or  directly  from  Zt 
that  is  at  T,  T'  the  points  of  contact  of  tangents  to  the  circle 
drawn  from  Z.  As  P'  moves  from  Tr  to  T,  H  increases,  or 
the  node  progredes,  and  as  P'  moves  from  T  to  T\  fl  decreases, 
or  the  node  regredes. 

/         N\ 
If,  then,  T  denote  the  angle  TCZl  cos"1-^1!  ;the  period  of 

the  progression  will  be  the  time  P'  takes  to  move  from  T'  to 

T,  that  is  — ~ — — ;  the  period  of  the  regression  will  be  the 
ft, 

2r 
time  P'  takes  to  move  from  T  to  T",  that  is  -T-;  and  the  node 


88  PLANETARY  THEOEY. 

is  stationary  whenever  \t  +  S^  —  S2  =  (2n  +  1)  TT  ±  T.  Also 
the  angle  through  which  the  node  oscillates  is  TZT'  which 
=  7r-2r. 

Secondly,  suppose  N^  =  N9.  Then  Z  will  lie  on  the  cir- 
cumference of  the  circle  which  P'  traces  out.  In  this  case  T 
and  T'  coincide,  and  the  node  after  becoming  stationary 
begins  to  move  in  the  same  direction  as  before. 

Thirdly,  suppose  N^  greater  than  N2.  Then  ^lies  within 
the  circle,  and  no  tangents  can  be  drawn  from  it  to  the  circle  : 
from  this  we  see  that  the  node  is  never  stationary. 

It  is  easily  seen  that  in  all  three  cases  the  maximum 
and  minimum  values  of  the  inclination  are  ZO  and  ZO',  that 
is  Nt  +  N2  and  N^  •*  N2  respectively.  And  in  all  cases  the 
node  moves  fastest  or  slowest  when  P'  coincides  with  0  or  0', 
that  is,  whenever  the  inclination  is  a  maximum  or  minimum. 

The  above  geometrical  construction  also  affords  a  proof  of 
the  theorem  of  Art.  83. 

Since  in  the  case  of  two  mutually  disturbing  planets  the 
quantities  hlt  Slt  N9,  82  are  the  same  for  both,  it  follows  that 
the  point  P"  for  the  second  planet  traces  out  a  circle  con- 
centric with  that  traced  out  by  P',  and  with  the  same  uni- 
form angular  Velocity  ht  ;  also  that  the  two  points  P',  P" 
always  lie  in  a  common  radius  vector  through  C.  Now  the 
angle  subtended  by  P'P"  at  the  centre  of  the  sphere  is  the 
inclination  of  the  two  orbits;  this  inclination  is  therefore 
very  nearly  constant,  as  P'P"  is  small  and  constant  and  very 
near  Z. 

85.  To  integrate  the  equation  for  the  longitude  of  the 
epoch. 

We  have  (Art.  61) 


(e*  -  tan'i)  +  Az  (e'*  -  tanY  ) 
at 

+  Asee'  cos  (to-  —  -cr')  +  A4  tan  i  tan  i'  cos  (ft  — 
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Now  from  the  formulae  of  Art.  72,  we  obtain 

e~  =  M?  +  Ml  +  2  J/^  cos  {(gl  -  </2U  +  7,  -  7,}, 

ee  cos  (-55-  —  CT')  =  M^M{  +  MJtfi 

In  like  manner,  from  the  formulas  of  Art.  79, 

tan2 i  =  N*  +  N*  +  2Ayy"2 cos  [Jij  +  8,  -  SJ, 
tan2  iv= Ar1/2+  Ar22  +2N^N2  cos  {/i^  +  8t  —  &J, 

tan  i  tan  i'  cos  (O  —  H')  =  N1Ni'+  N* 

+  A"2  (N^  +N{]  cos  {\*  +  Sl  -  B2}. 

If  these  values  be  substituted  in  the  expression  for  -=-  , 
it  takes  the  form 

•j-  =Bn  +  B:  cos  {(gl  —  gz)  1 4-  7l  —  72}  +  -#2  cos  {/i^  +  St  —  Sj, 

where  5/?,  .5^  and  Ba  denote  certain  constants.    Integrating, 
we  have 

p 

e  =  e0  +  Bnt  -t-  —      —  sin  {(g^  —  g2)  t  +7!  — 72} 


We  may  omit  the  term  #ft£,  if  we  consider  it  as  furnishing 
a  correction  on  the  mean  motion  n,  which  thus  becomes 
(1  4-  B]  n.  With  this  understanding 


If  this  expression  be  developed,  we  may  again  omit  the 
term  involving  the  first  power  of  t,  and  consider  it  as  affording 
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a  further  correction  to  the  mean  motion*.     Thus  we  shall 
obtain  a  series  of  the  form 


86.  In  the  Theory  of  the  Planets  this  inequality  is  insen- 
sible, but  in  that  of  the  Moon  it  amounts  to  upwards  of  10 
seconds  in  a  century,  forming  what  is  termed  the  secular  ac- 
celeration of  the  Moon's  mean  motion.     Thus  it  appears  that 
this  inequality  does  not,  as  its  name  would  seem  to  imply, 
contradict  the  general  theorem  of  the  invariability  of  the  mean 
motions,  since  it  is  due  to  a  variation,  not  of  the  mean  motion, 
(as  we  have  employed  the  term  in  the  preceding  pages,)  but 
of  the  epoch.     If,  however,  as  in  the  Lunar  Theory,  the  epoch 
be  omitted,  any  variation  in  the  mean  longitude  will  of  neces- 
sity be  thrown  upon  the  mean  motion  ;  only  in  this  case,  n 
will  not  be  given  by  the  equation  n*a*  =  p. 

87.  We  have  hitherto  supposed  the  planetary  motions  to 
be  referred  to  a  fixed  plane,  but  have  left  the  particular  plane 
undetermined.    In  practice  it  is  usual  to  take  the  position  of 
the  ecliptic  at  some  given  epoch,  as  for  instance  the  year 
1800  ;  but  since  it  is  to  the  true  ecliptic  that  astronomers 
refer  the  celestial  motions,  we  will  now  obtain  formulae  for 
determining  relatively  to  the  plane  of  the  Earth's  orbit,  the 
position  of  that  of  any  other  planet. 

Let  then  ra,  m'  denote  the  masses  of  the  Earth  and  the 
planet  considered,  and  suppose  the  orbits  of  m  and  m'  but 
little  inclined  to  each  other  and  to  the  fixed  plane  of  reference. 
Let  \,  X'  denote  the  latitudes  of  points  in  these  orbits  corre- 
sponding to  the  same  longitude  0±;  then  (see  fig.  to  Art.  13) 

tan  X  =  tan  i  sin  (0X  —  fl),    tan  X'  =  tan  i'  sin  (0±  —  fl'). 

*  The  advantage  of  thus  disposing  of  these  terms  arises  from  the  fact 
that  the  mean  motion,  as  determined  by  observation,  is  the  complete  coeffi- 
cient of  i  in  the  expression  of  the  mean  longitude* 
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Now  sinc.e  i  and  i'  are  very  small,  we  may  replace  tan  X, 
tanX'  by  X,  X'  respectively:  thus 

X'  -  X  =  tan  i'  sin  (9l  -  Of)  -  tan  i  sin  (0l  -  O) 
=  (tan  i'  cos  IT  —  tan  i  cos  fl)  sin  61 

—  (tan  i'  sin  11'  —  tan  i  sin  II)  cos  #x  , 
or,  with  the  notation  of  Art.  79, 

X'  -  X  =  (q'-q]  sin  9l-(p'-p)  cos  0l  .........  (1). 

Now  let  7  denote  the  inclination,  v  the  longitude  of  the 
node  of  the  orbit  of  m'  relatively  to  that  of  m  ;  then  approxi- 
mately 

X'  —  X  =  tan  7  sin  (91  —  v) 

=  tan  7  cos  v  sin  Ol  —  tan  7  sin  v  cos  0t  ..........  (2). 

Hence  equating  coefficients  of  sin  6^  and  cos  0l  in  equations 
(1)  and  (2), 

tan  7  cosz>  =  q'  —  q,    tan  7  sin  v  =p'—  p\ 
whence  tan2  7  =  (p  —  pf  +  (q  —  q)z, 

and  tan  v  =  —  —  *-  . 

q  -q 

These  expressions  determine  the  position  of  the  orbit  of 
m  relatively  to  that  of  m,  when  the  values  of  p,  p  ,  q,  and  q 
are  known.     Differentiating  them,  and  neglecting  small  quan- 
tities of  orders  higher  than  the  second,  we  obtain 

dry      (dp      dp\   .  (dq      dq\ 

-i  =  LJ^  -  ~TL   sm  v  +  \-ji  ~  -ji]  cos  v> 
dt      \dt      dtj  \dt      dtj 

dv  _  i  dp'      dp\  cos  v      fdq      dq\  sin  v 
dt      \dt       dtjt&ny      \dt       dt 


If  the  values  of  ~  ,  -=2  &Ct  fa  substituted,  these  equa- 


, 
tions  give  the  variations  of  7  and  v. 
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88.  It  has  been  proposed  to  employ  the  invariable  plane 
of  the  Solar  System  as  a  plane  of  reference  ;  and  since  it 
remains  absolutely  fixed,  it  would  afford  a  means  of  com- 
paring together  observations  separated  by  long  intervals. 
But  the  want  of  a  convenient  fixed  line  upon  it,  from  which 
to  measure  longitudes'"",  stands  in  the  way  of  its  practical 
application  ;  so  that  it  must  be  looked  upon,  at  all  events  for 
the  present,  rather  as  theoretically  interesting  than  as  avail- 
able for  astronomical  purposes. 

*  See  Pontecoulant's  Systeme  du  Monde,  tome  i.  p.  469  (2nd  edit.)- 

Note.  Gauss  in  a  memoir  (Determinatio  Attractionis,  d*c.}  contributed  to 
the  Transactions  of  the  Royal  Society  of  Science  of  Gottingen,  1818,  has 
stated  without  proof  the  following  theorem.  The  secular  variations  which 
the  elements  of  a  planetary  orbit  experience  from  the  perturbation  produced 
by  another  planet  are  independent  of  the  position  of  this  planet  in  its  orbit, 
and  would  be  the  same  whether  the  disturbing  planet  moved  in  its  elliptic 
o'rbit  according  to  Kepler's  Laws,  or  whether  its  mass  be  conceived  to  be 
equably  distributed  throughout  the  orbit  in  such  manner  that  equal  parts  of 
the  mass  are  now  assigned  to  such  parts  of  the  orbit  as  were  described  in 
equal  intervals  of  time,  provided  that  the  times  of  revolution  of  the  disturbed 
and  disturbing  planets  are  not  commensurable.  Gauss  has  also  shewn  how 
to  reduce  the  expressions  for  the  attraction  of  such  a  mass-orbit  at  any 
assigned  mass-point  not  on  the  orbit  to  calculable  forms  of  elliptic  inte- 
grals. [A.  F.] 


CHAPTER  VI. 

PERIODIC   VARIATIONS   OF  THE   ELEMENTS   OF   THE   ORBIT. 

89.  WE  come  now  to  consider  the  variations  produced 
by  the  periodical  terms  of  R.     These  are  called  Periodical 
Variations,  as  opposed  to  the  Secular  Variations  produced  by 
the  non-periodic  terms.      We    have  seen,  indeed,  that  the 
latter  are  for  the  most  part  periodical  in  form,  but  in  the 
Planetary  Theory,  the  term  Periodical  Variations  is  restricted 
to  those  we  are  about  to  consider  in  the  present  Chapter. 

90.  We  have  seen  (Art.  46)  that  the  general  type  of  a 
periodical  term  is  P  cos  {(pn  ±  qn)  t  +  Q],  where  P  is  a  func- 
tion of  a,  e,  i  \  and  Q  is  a  function  of  -or,  e,  H.     Now  such  a 

.,,  .    .,  .    dR     dR        jdB, 

term  will  produce  a  similar  term  in  -=-  ,  -j-  ,  and  -j-r- :  but 

da       ae  at 

a    term   of  the   form  P  sin  {(pn  ±  qn'}  t  +  Q]   in   ^— ,   -y-  , 

jjD 

and    -^=  .     If,  then,  these  be  substituted  in  the  equations 

d\.L 

of  Art.  39,  they  will  take  the  forms 


94  PLANETARY  THEORY. 


dt 

-T7  =  P7  sin  {(pn  ±  qri)  t  +  Q], 

where  PI}  P2,  &c.  are  functions  of  the  elements  of  the  dis- 
turbed and  disturbing  planets,  and  involve  the  first  power  of 
the  disturbing  mass. 

91.  In  integrating  these  equations,  we  may  in  general 
consider  the  elements  which  enter  in  the  right-hand  members 
as  constant  and  equal  to  their  values  at  the  epoch  from 
which  the  time  is  reckoned*. 

Let  then  a,  e,  OT,  &c.  denote  the  values  of  the  elements 
at  epoch,  Sa,  &e,  SOT,  &c.  their  periodical  variations  after  an 
interval  t :  then  integrating  the  above  equations,  and  omit- 
ting the  constant  terms,  we  have 

Sa  =  -  ^TT^  cos  {(Pn  ±  VV  t  +  Q}> 


7  cos  {(pn  +  qri)  t  +  0}, 
sin  {(pn  +  qri)  t  +  0}, 


pn  ±  q 


*  This  is  equivalent  to  neglecting  the  square  of  the  disturbing  force  :  see 
Art.  95. 
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p 

-*—  ,  sin  {(pn  4-  qri)  t  +  0}, 
ri 


pn     q 


p 

s—7  sin  {(pn  +  qri)  t+Q], 
ri 


pn  ±  q 
p 


pn  ±  q 


,cos{(pn±qri) 


Hence  it  appears  that  the  variations  produced  by  the 
periodical  terms  of  R  are  all  periodical  in  form. 

92.  It  will  be  seen  that  all  the  expressions  of  the  last 
Article  involve  the  divisor  pn  ±  qri,  while  £f  involves  the 
divisor  (pn  ±  qri)*.  If  then  it  should  happen  that  either 
pn  -f  qn'  or  pn  -  qri  is  very  small,  a  term  in  R  containing 
(pn  ±  qri)  t  in  its  argument,  though  of  a  high  order,  may 
have  a  sensible  effect  on  the  elements  of  the  orbit.  Now 
since  p  and  q  are  either  positive  integers  or  zero,  pn  +  qri 
cannot  be  small  unless  n  and  ri  are  small,  a  case  which  does 
not  occur  with  any  of  the  planets  :  but  we  have  instances 
in  which  pn  ~  qn  is  small  *. 

Since  the  period  of  such  inequalities  is  very  great  f  being 
-,)  ,  they  are  called  inequalities  of  long  period,  or  long 


pn  ~  qn 
inequalities. 

93.     To  select  such  terms  in  R  as  will  produce  the  prin- 
cipal inequalities  of  long  period. 

*  If  in  any  case  the  mean  motions  of  two  planets  were  exactly  commen- 
surable and  in  the  ratio  of  p  to  q,  the  corresponding  term  of  E,  as  we  have 
already  remarked  (Art.  48),  would  cease  to  be  periodical  and  would  form 
a  part  of  F,  but  no  instance  of  this  occurs  among  the  planets. 
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We  have  seen  that  the  dimension  of  the  principal  part 
of  the  coefficient  of  a  term  containing  (pn  ~  qn)  t  in  its  argu- 
ment is  p  ~  q  (Art.  50) ;  hence  if  we  can  find  two  integers 
p  and  q  nearly  in  the  ratio  of  n  to  n,  and  having  a  small 
difference,  the  corresponding  term  of  R  will  produce  an  im- 
portant long  inequality  in  the  elements  of  each  planet. 

In  the  case  of  Jupiter  and  Saturn  n  :  n  : :  5  :  2  nearly, 
and  5—2  =  3;  hence  there  is  a  long  inequality  arising  from 
a  term  in  R  of  the  form  Pcos{(2r&—  5ri)  t  +  Q},  the  prin- 
cipal part  of  P  being  of  the  third  order.  This  inequality  is 
interesting  in  an  historical  point  of  view,  having  long  baffled 
the  labours  of  mathematicians  and  appeared  inexplicable  on 
the  hypothesis  of  gravitation.  It  was  at  last  successfully 
explained  by  Laplace. 

For  the  Earth  and  Venus,  n  :  n  ::  8  :  13  nearly,  so  that 
there  is  a  long  inequality  arising  from  a  term  in  R  of  the 
fifth  order.  The  discovery  of  this  inequality  is  due  to  the 
Astronomer  Royal. 

Finally,  in  the  case  of  Neptune  and  Uranus,  n  :  ri  ::  1  :  2 
nearly,  hence  there  is  a  long  inequality  arising  from  a  term 
in  R  which  is  of  the  first  order. 

94.  Between  corresponding  terms  of  the  long  inequali- 
ties in  the  mean  motions  of  two  planets,  arising  from  the 
near  commensurability  of  n  and  ri,  there  is  a  simple  approx- 
imate relation. 

Let  m,  m'  be  the  masses  of  the  two  planets,  R,  R  their 
disturbing  functions :  then  by  Art.  8,  considering  only  the 
mutual  action  of  m  and  m',  we  have 

m'     m'  , 
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„,     m     m,     ,  ,       .        ,x 
R=-j--,(xx+yy  +  zz). 

We  shall  distinguish  the  first  and  second  terms  of  R  and 
R'  as  the  symmetrical  and  unsymmetrical  parts  respectively, 
since  the  co-ordinates  of  m  and  m'  are  involved  symmetrir 
cally  in  the  former  but  not  in  the  latter. 

Since,  then,  the  symmetrical  parts  of  R  and  R'  differ  only 
in  having  m  and  m'  interchanged,  if 

m'if  cos  {(pn  —  qn)  t  4-  Q] 

be  any  term  in  the  symmetrical  part  of  R,  that  of  R'  will  con- 
tain the  term 

mMcos  {(pn  —  qri)  t  +  Q}. 
Confining  our  attention  to  these  terms,  we  have  (Art.  39) 

d*£  =     3  ro'q  dR  =     Bna  d  (R) 
dtz  p     de  p       dt 

^2:  m'M  sin  {(^n  —  gn)  t  +  Q}, 


therefore   8?  =  -  ,  /N2  sin  {(^n  -  ?n')  t  +  Q}. 

yit     (pn  — 

d-    -i    i      c>v/     Sn^dq       mM 
Similarly,  Sf=       7-^ 


'' 


mun'ap        mfjna 
Hence      ~  =  --     /2  /   =  ---  ^7 
o^          m/^n  a  c[        wifjun  a 

approximately,  since  qri  is  nearly  equal  to  pn  ;  therefore 


m 


or,  since  /&'  differs  from  /A  by  a  quantity  of  the  order  of  the 
c.  P.  T.  7 
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disturbing  force,  the   square  of  which  we  are   neglecting, 
we  have  to  the  first  order, 


m 


the  required  relation. 

The  same  relation  is  also  approximately  true  in  the  case 
of  terms  arising  from  the  unsymmetrical  parts  of  R  and  R*. 
For  denoting  these  by  R^  and  R{  respectively,  we  have 


Now  the  equations  of  motion  of  the  planet  m'  referred  to 
rectangular  axes  are 


_ 

d?       r*  ~  '  daf  '       ' 

and  since,  by  the  principles  of  the  method  of  the  variation  of 
elements,  theybrm  of  solution  of  these  equations  is  the  same 
as  it  would  be  if  R'  were  zero,  it  follows  that,  if  the  differ- 
ential coefficients  be  taken  as  if  the  elements  were  constant, 
we  shall  have 

x       _l_dW     y^     _LA'     *?__  _1^?' 
r'3"         //,'  d?  '   r's  ~      fif  df  '   7*~     ^'W'1 

and  therefore,  with  this  understanding,  that 


*  For  the  demonstration  of  this  we  are  mainly  indebted  to  The  Theory  of 
the  Long  Inequality  of  Uranus  and  Neptune:  an  essay  which  obtained  the 
Adams  Prize  for  the  year  1850.  By  E.  Pierson,  M.A. 
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O*         .t        1  TV, 

Sinnlarly,         R,  = 


Now  any  term  in  K1  containing  (pn  —  qn) t  in  its  argu- 
ment can  arise  only  from  the  combination  of  terms  in  x,  y, 

and  z,  containing  put  in  their  arguments  with  terms  in  -j-2-  , 

dt 

d*y  d^z' 

— - ,  and  -y-j  containing  qnt.     Suppose,  then,  x  and  x'  when 
dt  dt 

developed  in  terms  of  t  and  the  elements  to  contain  respec- 
tively the  terms 

L  cos  (put  +  I),     L'  cos  (qnt  +  1'). 
Hence  the  product  x  -^-  will  contain  the  term 

-  -  LL'qV cos  {(pn  -  qn')  t+l-  I'}, 
and  the  product  x'  -^  the  term 

—  -  LL'p*n* cos  {(pn  -qn')t  +  l-l'}: 

the  coefficients  are  in  the  ratio  qV  to  p*n*.     Similarly,  the 
coefficients   of  the   same  cosine   in  y—rf  and  z  -^  are  to 

those  in  y'  -~  and  z'  -^  in  the  same  ratio. 

Hence  if       — r  MqV  cos  {(pn  —  qn')  t  +  Q} 
be  any  term  in  Rlt  then  Rf  will  contain  the  term 

—  Mp*n* cos  {(pn  —  qn)  t  +  Q}. 

7—2 
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Confining  our  attention  to  these  terms,  we  have 
' 


s.n 
/-I/A       (pn  —  qn  )2 


Hence  ^>  =  --  >    =  -.  —  7-7  nearly, 

of          map         mud 


m 


m  */a 
the  square  of  the  disturbing  force  being  neglected. 

By  means  of  this  relation,  when  one  of  the  long  inequali- 
ties is  known,  the  other  may  be  calculated  :  it  may  also  be 
used  as  a  formula  of  verification. 

95.  We  have  remarked  that  in  integrating  the  equations 
of  Art.  90,  we  may  in  general  consider  the  elements  which 
enter  in  the  right-hand  members  as  constant  and  equal  to 
their  values  at  the  epoch  from  which  the  time  is  measured. 
In  the  case,  however,  of  inequalities  whose  periods  are  very 
long,  the  secular  variations  of  the  elements  in  the  interval 
produce  a  sensible  effect.  In  order  to  take  account  of  these, 
we  may  integrate  our  equations  by  parts,  considering  the 
elements  variable  ;  and  then  substitute  their  values  as  calcu- 
lated by  the  method  of  Art.  62.  For  example,  consider  the 
equation 


=  P  sin  X,  suppose. 

Integrating  by  parts,  and  remembering  that  n  is  constant 
with  regard  to  secular  variations,  we  have 
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1  dP    . 

--  sin  X 


- 

dt         pn  —  qn  (pn  —  qn) 

1          d*P 


therefore  S   =  -  ~          sin  X  - 


1         d*P  . 

+  7~  —  7T4-T75-S 

(pn  —  qn)*  dt* 

dP  1          cFP. 

-7   cos  X  +  ,  ,.4  -3-sin  X  +  .  .  . 


(pn  —  qn)3  dt  (pn  —  qn')* 

1         rf'P 

•f  r~   -  774-77  s 
(pn  —  qri)*  dt' 

|    . 

"  " 


2          cZP          I 
—  -,  -  —*  —rr  —  •••  f  cos  X. 
(pn-  qrif  dt  j 


In  this  equation  P,  -3-  ,  -^  ,  &c.  are  functions  of  the 
elements  ;  their  values  may  be  calculated  by  the  formulae  of 
Art.  62.  It  may  be  noticed  that  P  is  of  the  first  order,  -7-  of 

d2P 
the  second,  and  -p-  of  the  third  of  the  disturbing  force: 

for  -=-  ,  being  found  from  P  by  differentiation,  will  involve 
dt 

the  differential  coefficients  of  the  elements,  which  are  them- 

cfP 

selves  of  the  first  order  ;  and  similarly  for  -7-5-  . 

96.     Having  now  completed  our  account  of  the  methods 
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of  treating  the  secular  and  periodic  variations  of  the  elements 
of  the  orbit,  we  will  say  a  few  words  on  the  distinction  be- 
tween them.  In  the  first  place  we  may  observe  that  the 
periodic  variations  involve  the  mean  longitude  of  the  dis- 
turbed and  disturbing  planets,  and  therefore  depend  chiefly 
upon  the  configuration  of  the  planetary  system.  On  the  con- 
trary the  secular  variations  depend  solely  upon  the  values  of 
the  elements.  The  latter  class  of  variations  take  place  with 
extreme  slowness,  so  that  if  these  only  existed,  a  considerable 
time  must  elapse  before  the  deviation  of  the  planet  from 
elliptic  motion  became  appreciable.  On  the  other  hand,  the 
periodic  variations  (such  at  least  as  are  rapidly  periodic) 
"  are  in  their  nature  transient  and  temporary :  they  disappear 
in  short  periods,  and  leave  no  trace.  The  planet  is  tempo- 
rarily drawn  from  its  orbit  (its  slowly  varying  orbit),  but 
forthwith  returns  to  it,  to  deviate  presently  as  much  the  other 
way,  while  the  varied  orbit  accommodates  and  adjusts  itself 
to  the  average  of  these  excursions  on  either  side  of  it ;  and 
thus  continues  to  present,  for  a  succession  of  indefinite  ages, 
a  kind  of  medium  picture  of  all  that  the  planet  has  been 
doing  in  their  lapse,  in  which  the  expression  and  character 
is  preserved;  but  the  individual  features  are  merged  and 
lost*."  On  this  account  it  is  convenient  to  suppose  the 
planet  to  move  in  an  ellipse,  the  elements  of  which  are  cor- 
rected for  secular  variations  only,  and  to  take  account  of  the 
periodic  variations  by  applying  small  corrections  to  the  radius 
vector  and  longitude  as  calculated  from  the  elliptic  formulae. 

97.  We  will  accordingly  shew  how  by  means  of  the 
periodic  variations  of  the  elements,  the  corresponding  in- 
equalities in  the  radius  vector  and  longitude  may  be  calcu- 
lated. If  we  take  for  our  plane  of  reference  the  position  of 

*  Herscliel's  Outlines  of  Astronomy,  10th  edit.  Art.  656. 
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the  plane  of  the  orbit  of  the  disturbed  planet  at  the  epoch 
from  which  the  time  is  reckoned,  the  inclination  will  be  of 
the  order  of  the  disturbing  force,  and  therefore,  if  we  neglect 
the  square  of  the  latter,  we  may  also  neglect  the  square  of 
the  former. 

98.  To  calculate  the  periodic  variations  in  radius  vector. 

Let  Ba,  Be,  S-cr,  &c.  denote  the  periodic  variations  in 
a,  e,  <&,  &c.,  and  let  Br  be  the  corresponding  variation  in  r ; 
then 

£       dr  5,       dr  .        dr  ~         dr  .  „     dr ~ 
gr  =  ~-Ba  +  -j- 8e  + -j- tier  +  -^6$  +  jr  6e, 
da          de          d^  d£          de 

in  which  the  square  of  the  disturbing  force  is  neglected, 
since  this  would  be  introduced  by  the  squares  and  products 
of  Ba,  Be,  &c.  The  values  of  Ba,  Be,  &c.  have  been  found  in 

dr      dr 

Art.  91,  those  of  -j-  ,    -j- ,  &c.   may  be   obtained  from  the 
da     de 

equation  (Art.  40) 

r  =  a\l  +  ~e2  —  ecos(f-f  e  —  OT)  —  ^e2cos2  (f +e  —  OT)—  ...  -. 
I       *  *  J 

99.  To  calculate  the  periodic  variations  in  longitude. 

These  might  be  found  in  the  same  manner  as  the  varia- 
tions in  radius  vector,  but  they  may  also  be  deduced  from 
them :  we  proceed  to  obtain  a  formula  for  this  purpose. 

We  have  tl  =  p>     (Art  22)' 

and  0-00  =  ft-O0; 

therefore  -j-  =  —z  +  (1  —  cos  t)  -j-     (see  Art.  29) 

* 
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since  (1  —  cos  i)  —7-  ,  being  of  the  order  of  the  cube  of  the 

disturbing  force,  may  be  neglected. 

Let   Sr,  80,  and  Sh  be  corresponding  variations  in  r,  0 
and  h  ;  then 

d(0+$0)=  h  +  Sk 

dt 

d6     dW     h 

°r      S+^-? 

h 


neglecting  the  square  of  the  disturbing  force  ;  therefore 

Sh 


which  gives  the  variations  in  longitude.     The  value  of  $h 
may  be  found  from  the  formula 

dh^dR     dR 
dt  ~  de      dm' 

For  the  periodic  variations  in  latitude,  we  refer  to  Ponte- 
coulant's  Systeme  du  Monde,  Tome  I.  p.  492. 

100.     As  an  example  of  the  processes  of  this  Chapter,  we 
will  calculate  the  variations  in  radius  vector  and  longitude 
due  to  the  term  m'Me  cos  {(n  —  2ri)  t  +  e  —  2er  +  -cr}  in  E. 
Considering  this  term  only,  we  have 

R  =  m'Me  cos  {(n  -2n')t  +  e-  2e'  +  ^} 
=  m'Me  cos  X,  suppose. 

dR        ,dM  dR       ,,, 

Hence     -^—  =  m  -7-  e  cos  X,     -j-  =  mm  cos  X, 
da  da  de 

dR  ,„.-    .   .        dR  f^f    .    _ 

—-  =  —  m  Me  sin  X,      -j—  =  —  m  Me  sin  X, 
de  di* 


ao 
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Substituting  these  in  the  formulae  of  Art.  59,  and  neg- 
lecting small  quantities  of  orders  higher  than  the  first,  we 
have 


da  /71,     .    . 

-y-  =  --  m  Me  sm  X, 
dt          ft 

de     na    ,  ,,  .    . 
—  -  =  —  m  M  sm  X, 
dt       t 


'sr     na 

e  —rr=  — 

dt       jt, 


de         2na?    ,dM  I  na      ,,, 

-TT  =  --  -m  —j—  e  cos  X  +  „  —  em  If  cos  X, 

dt  da  2   j, 


•jz  =  --  m  Me  sm  X. 
at        /it 

By  integration  we  have 

2m'  M   na^e 

oa=  —     --  —  ,  cos  X, 
fji     n—2n 

~          m'M     na 
de  =  --    --  —,  cos  X, 
/-t    n  —  2n 

m'M     na       .    . 
--  —  ;  sm  X, 
fju    n  —  2n 

^       /I  m'M     2m  a  dM\     nae      .    . 

oe  =    ~  ----  -y-  -  —  —,  sin  X, 

\2   /^          p,     daJn—Zn 

<*,         3m'  M     n*ae 


s  —  --  •  7 

^     (n  — 

Also 


7\2  sm  X. 
z 
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therefore,  small   quantities  of  orders   higher  than  the  first 
being,  neglected, 


dr 

-j-  =  a  [e  -  cos  (f  +  e  -  w)  -  e  cos  2  (f  -f  e  —  «r)}, 


=  -  a  {e  sin  (f  +  e  -  w)  -f  e2  sin  2  (f  +  e  -  «r)}, 


T-  =  a0  sin  (f  +  e  -  «•). 


XT        *       dr  2       dr  s       dr  ~      ,  dr          dr 
Now  dr  =  -T-  6a  +  ^-de  +  T-  OOT  +  3 
aa          ae          ats- 


n  — 


,  cos  \ 


%f  ...  x  ^/cx 

-  s—  ,  cos  X  {e  —  cos  (£+  e  —  «r)  —  e  cos  2  (f  +  e 
u»     n  —  £n 


Since  we  are  neglecting  the  square  of  the  disturbing  force, 
the  elements  in  this  equation  may  be  considered  as  constant, 
and  therefore  nt  written  for  £:  we  have  then,  restoring  to  X 
its  value 

.,       m'M    na?  0  f/         ,^  .  n 

5~?  cos  2  {(n  -  n  )  <  +  e  -  e  } 


IJL     n  —  Z 
m'M   na?e 


f/       0  ,N  .  .         0  , 
cos  \(n  —  2n  )  t  +  e  -  2e 
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m'M   na?e  f,~        ^  f\  ,     ^       «  /        i 

r-  >  cos  {(3w  -  2n  )  £  -f  3e  -  2e  -  w), 


//,     n  —  2 
which  is  the  variation  in  radius  vector. 

101.     To  calculate  the  variation  in  longitude,  we  shall 
employ  the  equation 


~dt~  r*        r3    ' 
dR 


therefore  oh  =  --  ^-7  cos  X  ; 

71—  2ft 


8A,          2m'Me 

therefore  -g-  =  -=-7  —     T-TT  cos  X 
r        a  (n  —  tn) 


2m  M    n*ae          .,       0  ,.  ,  0  , 

~~'cosw~      +6~ 


Also         = 

m'M     n* 


5  m'M  (/       0  /\  ,  ,         0  / 

^  ---  s—  /  cos  f  (w  —  2w  )  ^  +  e  —  2e 
2    /t    n  —  2n 

5  m'lf    n*ae 


'  2    n    n-'. 
Hence  by  substitution 


r/0        0  /x  .  ,  0       0  , 
'  C°S  {  (      "       }    +      ~ 
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3m'  M 


5m'M  f/0       0  ,, 

~~COS^        ) 


By  integration 
tnfjf 


w  — 
3m'  Jf      n 


(n  -  2n 
5m'M 


0  ,.  /0          v 
/>t     (n  -  2n)  (3n  —  2rc  ) 

which  is  the  variation  in  longitude. 

In  the  case  of  Uranus  and  Neptune,  since  n  :  n  nearly 
as  2  :  1,  the  term  we  have  been  considering  is  important  in 
the  theory  of  the  long  inequality. 


CHAPTER   VII. 

DIRECT  METHOD   OF   CALCULATING  THE  INEQUALITIES  IN 
EADIUS  VECTOR,   LONGITUDE,  AND  LATITUDE. 

102.  IN  the  calculation  of  the  planetary  inequalities,  we 
have  hitherto  employed  exclusively  the  method  of  the  Varia- 
tion of  Elements,  but  there  is  another  method  of  solving  the 
problem,  which  demands  our  attention.     It  consists  in  obtain- 
ing equations  for  calculating  the  inequalities  in  radius  vector, 
longitude,  and  latitude  directly  from  the  equations  of  motion. 
This  method  is  indeed  the  simplest  to  employ  in  the  case  of 
periodic   variations  of  short  period,  that  of  the  preceding 
Chapters  being  the  most  convenient  for  the  calculation  of 
secular  variations  and  long   inequalities.      For,  since  these 
latter  take   place  with    extreme   slowness,  the  elliptic  ele- 
ments, when  once  corrected  for  them,  continue  for  a  consi- 
derable period  to  represent  the  actual  motion ;  while,  in  the 
former  case,  the  values  of  the  elements  change  rapidly,  and 
the   motion   cannot   for   long   be   represented  by  the  same 
ellipse.     We  proceed,  then,  to  the  direct  method  of  calcu- 
lation*. 

103.  If  riy  6lt  and  z  denote  the  projected  radius  vector, 
longitude,  and  distance  from  the  plane  of  reference,  of  the 
planet,  we  have  (see  Art.  9)  the  equations  of  motion 

*  The  two  methods  are  sometimes  distinguished  as  those  of  Lagrange 
and  Laplace,  but  in  the  Mecanique  Celeste  we  find  both  employed. 
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dR 
W' 


-  _ 
•S?  "     7^     dz  ' 

If  we  take  for  the  fixed  plane  of  reference  the  position  of 
the  plane  of  the  orbit  at  the  epoch  from  which  the  time  is 
reckoned,  the  inclination  (as  we  have  remarked  in  Art.  97) 
will  be  the  order  of  the  disturbing  force,  the  square  of  which 
will  be  neglected.  Now  it  will  be  seen  on  referring  to  Art. 
42,  that  rl  and  6^  differ  from  r  and  6  by  quantities  depending 
upon  the  square  of  the  inclination  :  hence  in  the  above  equa- 
tions, we  may  replace  r,  and  0:  by  r  and  0  respectively.  Also 
if  \  denote  the  latitude  of  the  planet,  we  have 

z  —  r  sin  \. 
Hence  our  equations  of  motion  become 


_ 
dt(r  dt)~dO  ......... 


dt*  dz 


104.  As  a  first  approximation,  let  values  of  r,  6  and  X 
be  obtained  from  these  equations  by  neglecting  the  disturbing 
force,  and  let  r  +  £r,  6  +  80,  X  +  8X  denote  the  true  values  of 
these  co-ordinates  ;  then  £r,  §6  and  B\  will  be  very  small 
quantities,  of  the  order  of  the  disturbing  force  :  they  are 
termed  the  perturbations  in  radius  vector,  longitude,  and 
latitude.  We  proceed  to  investigate  equations  by  means  of 
which  these  quantities  may  be  determined. 
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105.    To  obtain  the  equation  for  the  perturbation  in  radius 
vector. 

From  equations  (1)  and  (2)  of  Art.  103,  we  obtain 

HEFT'   \dt)  =^  +  2}(dr~di  +  d0  dt)  +  G 

t  +  C ....(4). 

Multiply  (1)  by  r  and  add  to  it  (4) :  thus 


r         dr         j    dt 


dt'         r        <r  dr         J     dt 

If  the  disturbing  force  be  neglected,  this  equation  be- 
comes * 


Let  a  value  of  r  be  obtained  from  this  equation,  and  let 
r  +  Sr  denote  the  true  radius  vector  :  then  if  we  agree  to 
neglect  the  square  of  the  disturbing  force  in  our  next  ap- 
proximation, it  will  be  sufficient  to  write  r  +  &r  for  r  in  those 
terms  of  (5)  which  do  not  involve  the  disturbing  force  :  also 
since  8r  is  itself  of  the  order  of  the  disturbing  force,  its  square 
may  be  neglected.  Hence  from  (5) 


Sr  dr         J     dt 

§r)     2p     2^  ,        . 
72  —-^   =  —  --  ~  Sr  -f  2r  - 

dt  df          r       r4  .  dr 


,  nd*(r§r)     2p     2^  ,        .    dR 

therefore      7;2  '  +  2  —-^   =  —  --  ~  Sr  -f  2r  -j- 

4 
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hence  by  (6), 


which  is  the  equation  for  the  perturbation  in  radius  vector. 
We  may  express  the  right-hand  member  in  a  more  convenient 
form,  for  since  (Art.  42) 

r  =  a  (1  +  u\ 

dR_dRdr^  .dRm 

da      dr  da  dr  ' 

dR        dR 

therefore  r-T~  =  a^r  ; 

dr         da 

d(R)        dR 

also  \.     =  n  -y-  . 

Hence  our  equation  becomes 

dz  (rBr)     A^  /  PS  x         dR*    rt     CdR  , 
— W-*  4-  —3  (for)  =  a  -j — H  2»  I  -r-  «^. 

^         r3V  da          J  de 

106.  To  obtain  the  equation  for  the  perturbation  in  longi- 
tude. 

We  have  from  equation  (1)  of  Art.  103, 

r  dt2     r3     r  dr 

As  before,  let  a  value  of  6  be  obtained  from  this  equation, 
the  disturbing  force  being  neglected,  and  let  0  +  80  denote 
the  true  value  of  the  longitude :  then  writing  r  +  Br  for  r 
and  0  4-  80  for  0,  and  neglecting  the  square  of  the  disturbing 
force,  we  have 


___,   _ 
dt  dt  ~r  de      r*  df     r4          r  dr  '' 

Zd0     ,         dR        dR 
but  r=h>    r-  =  a> 
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.,       -       07  -  dV     3/K*         dR 

therefore  2h  -j-  —  r  —,^  —  or  -rr2  — £  cr  —  a  -=- 
2  2        2 


d  /  d$r      ^   dr\      3^     ~          dR 
==TAr^Tr-°r-j-)--f~.rSr-a-T-. 
dt  \    dt  dtj      r3  da 

This  equation  will  become  integrable  if  we  eliminate  the 
term  — ^ .  r$r  by  means  of  the  equation  for  the  perturbation 
in  radius  vector.  We  have  from  that  equation 

df         r3  *  da          ]  de       ' 


therefore  by  addition, 


d/  dSr     ,  dr 


dt\jM       "dt)~    "da 

?dR 

de 

dR 


;)- 

r^T? 
dt 


^      r       ois-a     Od(r8r)        dSr     «,  dr          [ 
therefore  2/iS^  =  3  —~—^-  +  r—^  --  Sr  -y-  -  4a  I 
dt  dt          dt         J 


, 

-j-  dt 
da 


the  arbitrary  constant  being  considered  as  included  in  the 
sign  of  integration  ;  therefore 

oisa     Ad(r*r)     e>$  dr          (dRj*  HdRj# 

2hSv  —  4i    ^        —  28r  -=-  —  4a  I  -j~  dt  —  6n  1  1  -j-  dt  , 
dt  dt          J  da  JJ  de 

0d(rSr)      .   dr  (dR  ,  [fdR  72 

or         h§6  =  2  —V—  -  Sr-rr  -  2a     -r~  dt  -  3n      -j-  dt2, 
dt  dt          J  da  JJ  de 


which  determines  the  perturbation  in  longitude,  when  that  in 
radius  vector  is  known. 

c.  P.  T.  8 
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107.  To  obtain  the  equation  for  the  perturbation  in 
latitude. 

From  equation  (3)  of  Art.  103,  we  have 

(f  (r  sin  X)     p  (r  sin  X)  __  dR 
~f  ~rT~      ~~fa*   .. 

Since  the  plane  of  reference  is  supposed  to  coincide  with 
the  position  of  the  plane  of  the  orbit  at  the  epoch  from  which 
the  time  is  reckoned,  we  have  at  the  epoch.  X  =  0  :  hence, 
denoting  by  SX  the  latitude  at  time  t,  our  equation  becomes 


, 

dz  ' 

which  is  similar  in  form  to  the  equation  for  the  perturbation 
in  radius  vector. 

108.     To  integrate  the  equation  for  the  perturbation  in 
radius  vector. 

The  equation  is  (Art.  105) 

tf(rBr)  ,  /*,  .  ,        dR     .   -fdR,. 

*   '    —  }jir^  +  ^(™r)=a-j--  +  %n    -j-dt- 
dt          r  ^  da          J  de 

Let  us  consider  a  term  in  R  of  the  form 
Pcos{(pn-qri)t+Q}> 

where  P  is  a  function  of  the  mean  distances,  excentricities, 
and  inclinations,  and  Q  of  the  longitudes  of  the  perihelia, 
nodes,  and  epochs  :  then  uniting  this  term  with  the  non- 
periodic  part  of  R,  which  we  have  denoted  by  F,  we  have 

R  =  F  +  Pcos  {(pn  -  qn)  t  +  Q}; 

dR     dF    dP       ,, 
therefore  =       +    -  cos  {(pn  - 
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since  F  does  not  contain  e ;  therefore 

nP^ 

dR  -,.  de          .,  i\  ,      ™         i 

— —  dt  =  ;cos  {(pn  —  qn)t  +  Ql  +  mg, 

de          pn-qn 

where  g  is  an  arbitrary  constant.  It  may  appear  superflu- 
ous to  introduce  this  quantity,  since  an  arbitrary  constant  G 
has  already  been  added  in  Art.  105  ;  but  its  introduction  is 
merely  equivalent  to  a  small  change  in  the  value  of  (7,  and 
any  alteration  which  does  not  interfere  with  the  first  ap- 
proximation, obtained  by  neglecting  the  disturbing  force,  is 
of  course  permissible.  Thus  g  is  a  purely  arbitrary  quantity 
which  might  have  been  omitted,  but  which  we  shall  find  it 
convenient  to  retain,  leaving  its  value  to  be  assigned  here- 
after. 

fj  7?          r  fj  7?  d  Ji1 

Hence         a  -=-  +  2n  I  -j-  dt  =  2m  a  +  a  -j- 

da          J  de  *          da 


a°f+ 2«  f  cos  {(pn  -  qn)  t  +  Q} 

da     pn—qn  ) 

7  Tjl 

=  2m 'g  +  a  -v-  +  P±  cos  {(pn  -  qn')  t  -t-  Q}, 


-r,  US.L  d€ 

suppose,  where       P1  =  a^  +  -—^ 
Again  (see  Art.  13), 

r  =  a  \  1  +  ~  e2  —  e  cos  (nt  -f  e  —  -cr) 
L 


—  -  e2  cos  2  (nt  +  e  —  w)  —  . . .  I , 


therefore       -3  =  -^  jl  +  3e  cos  (n^  +  e-w)  + 

=  ?i2  {1  +  3e  cos  (n^  +  e  —  -sr)  +  . . .}, 

since  nzo?  =  /^. 

8—2 
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Hence,  by  substitution,  the  equation  for  the  perturbation 
in  radius  vector  becomes 


PiCoaKpn-qn^t+Q] 

\A/U  W\,V 

—  ri* r&r  {30 cos  (nt  +  €  —  &)  +  ...}.  ' 

109.  This  equation  must  be  solved  by  successive  ap- 
proximation, as  in  the  Lunar  Theory.  By  omitting  all 
small  quantities,  we  obtain  a  first  approximation  to  the  value 
of  rSr;  this  being  substituted  in  the  second  member,  and 
small  quantities  of  orders  higher  than  the  first  neglected,  we 
obtain  a  second  approximation,  which  will  be  correct  to  the 
first  order.  In  like  manner,  a  third,  and  higher  approxima- 
tions may  be  obtained. 

On  referring  to  Art.  59,  it  will  be  seen  that  small  quanti- 
ties of  the  second  order  being  neglected, 

1 

~2m    °' 

hence,  neglecting  all  small  quantities,  the  equation  of  the 
preceding  Article  becomes 

—^-P  +  n\  rSr  =  Zih'g  +  ^  a  -^  +  Pt  cos  {(pn  -  qri)  t  +  Q}. 


The  integral  of  this  equation  is 


+    2  —  /  —  '-  —7x2  c°s  {(pn  -  qn')  t  +  Q} 
n2  —  (pn  -  qri)2 

+  Acos(nt-B), 
where  A  and  B  are  arbitrary  constants.     Since,  if  all  small 
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quantities  be  neglected,  r  =  a,  we  have  as  a  first  approxima- 
tion, 


1   /0    ,     ,  m'     dC0\ 
=  —   %mg  +  —  a  -j-0 
n*a\      *       2       da  J 


H — r~2 — 7— 1 ^27 cos  ((Pn  ~ 

a  (n2  —  (pn  —  qny] 

+  -cos(nt-B). 
a 

110.     We  may,  however,  omit  the  last  term:   for,  con- 
sidering this  only,  the  radius  vector  of  the  planet  becomes 

{A  \ 

1  —  e  cos  (nt  +  e  —  or)  +  -,  cos  (n£  —  B)  +  . . .  j- 

A 

—  a  [1  —  {e  cos  (e  —  OT) j  cos  .5}  cos  wi 

tt 

-4 

+  {e  sin  (e  —  w)  +  -^  sin  5}  sin  ^  +  . . .] 
fl 


=  a{L  —  ^cos 
if  el  cos  (e  —  VT^)  =  e  cos  (e  —  -or)  —  ^  cos  J5, 

tfj  sin  (e  —  orj  =  e  sin  (e  —  OT)  +  -^  sin  $, 

a 

from  which  ^  and  ^  may  be  determined. 

Now  since  the  ellipse  upon  which  our  approximations  are 
based,  has  been  obtained  by  neglecting  the  disturbing  force, 
we  may  in  the  elliptic  formulae  replace  e  and  TS  by  et  and  OTJ 
respectively,  since  they  differ  by  quantities  of  the  order  of  the 
disturbing  force.  If  this  be  done,  our  first  approximation 
becomes 

_!/      ,       m'    dC\ 
n£a\  2      da) 

*    ^>v\ cos  i(Pn  -  ^')  t+Q}- 
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111.  In  order  to  obtain  a  second  approximation,  this 
value  must  be  substituted  for  Sr  in  the  right-hand  member  of 
the  equation  of  Art.  108.  Also  since  the  square  of  the  dis- 
turbing force  is  neglected,  we  may  write  e^  and  TZ>I  for  e  and  -cr 
in  this  equation.  We  will  write  for  brevity 

$r  =  L  +  P2  cos  {(pn  -  qri)  t  +  Q}.' 

Substituting  this  in  the  equation  of  Art.  108,  and  omit- 
ting those  terms  which  have  produced  the  first  approxima- 
tion*, we  have 
d*  .  rSr       o      . 


=  —  %n*ael  cos  (nt  +  e  —  -c^)  [L  +  P2  cos  {(pn  —  qri)  t  -f  Q}~\ 

=  —  Sn9ae1  L  cos  (nt  +  e  —  -c^) 
g 

-  £  wX  P2  cos  [{(p  +  1)  n  -  qri} 

-  1  wX  P2  cos  [{(p  -l)n-  qri} 

112.     On  the  form  of  this  equation,  we  have  an  import- 
ant remark  to  make.     In  consequence  of  the  term 
—  S^agj  L  cos  (nt  +  e  —  wj, 

its  integral  will  contain  the  term 
g 

-  aejUt  L  cos  (nt  +  e  —  -orj. 
28 

Thus  we  are  met  by  a  difficulty  :  our  equations  have  been 
formed  on  the  hypothesis  that  the  square  of  £r  is  small 
enough  to  be  omitted,  whereas  here,  we  have  a  term  capable 
of  indefinite  increase.  This  term,  then,  if  retained,  would 
ultimately  vitiate  the  whole  approximation.  The  difficulty 
might,  as  in  the  Lunar  Theory,  be  obviated  by  writing 

*  These  terms  are  omitted  for  the  sake  of  brevity  :  in  order,  therefore,  to 
obtain  the  complete  second  approximation,  we  must  add  to  the  integral  of  the 
above  equation  the  result  of  the  first  approximation. 
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cn  for  n  in  the  elliptic  formulae,  which  amounts  to  supposing 
the  perihelion  to  be  in  motion.  Its  motion  is  however  better 
found  by  the  method  of  the  variation  of  elements.  Indeed  it 
may  be  shewn  that  such  ternis  as  those  we  are  considering 
lead  to  the  formulas  which  have  already  been  obtained,  for  the 
secular  variation  of  the  elements*.  We  may  accordingly 
altogether  neglect  such  terms  if  we  suppose  the  elements  of 
the  ellipse  on  which  our  approximations  are  based  to  have 
been  previously  corrected  for  their  secular  variations. 

113.     With  this  understanding,  the  complete  integral  of 
the  equation  of  Art.  Ill  will  be 

3 


. 


g 


2n?-{(p+l)n-qn'}* 

cos  [{(p  + 1)  n  —  qn}  t  +  Q  +  e  — 

SvPno     P 
II  UC1  JTa 


cos  [{ (p  - 1)  n  -  qn"}  t  +  Q  -  e  +  trj 
+  A  cos  (nt  -  B). 

If  this  be  added  to  the  result  of  the  first  approximation, 
we  obtain  for  a  second  approximation 

1  /  '  m'     dC0\  ,  P, 

rcr  —  —,   2m  q  +  -^  a  -y-5    +  -5 — 7 —       — T\I 
?i'  \  2       da  J      n2  —  (pn  —  qn  )2 

cos  }(p*  —  ^w')  t  +  Q} 

q  /M2/7/>      P 


cos  [{(p  +  1)  n  -  qn}  t  +  Q  +  €  -  «rj 


cos  [{(jp  -  1)  n  -  2»'}  t+Q-e  +  vt] 
+  A  cos  (ni  —  J5). 

*  It  is  thus  that  the  Secular  Variations  are  first  obtained  in  the  Mecanique 
Celeste.     See  Pont^eoulant's  Systeme  du  Monde,  Supplement  au  Livre  II. 


120  PLANET  AEY  THEORY. 

The  arbitrary  constants  may  be  disposed  of  as  in  the  first 
approximation.  In  order  to  obtain  a  second  approximation 
to  the  value  of  Sr,  it  is  only  necessary  to  multiply  the  right- 
hand  member  of  the  above  equation  by 

-{l+el  cos  (nt  +e  —  OTJ)), 
neglecting  e*. 

114.     To  calculate  the  perturbations  in  longitude. 
We  have  (Art.  106) 

7^/1     ^d.r^r      £   dr     0    fdR  ,       0    [t 

h8v  =  2  —  j-  --  8r  -Y  —  2a     -,-  at  —  3n  I 
dt  dt         J  da  JJ 


tdR 
-7- 
de 


Taking  for  simplicity,  the  first  approximation  to  the  value 
of  rSr,  which  has  been  obtained  by  neglecting  the  first  power 
of  the  excentricity,  we  have 

s       1  /0    ,     ,  m      dC 
r$r  =  -«   2m  q  +-^-  a  - 


«  -  7  ,      , 

f|H  2        c^a  /      ri*-(pn  -  qn) 

cosKpn- 
therefore 

0  cZ  .  rSr          2P  (pn  —  on)    .    ,  ,  /N  .      ^^ 

2  —  -j—  =  --  2-3^  -  ^-7(2  sm  {(on  -  «r)  t  +  Q)  : 
<ft  nz-(pn  —  qny 

also,  neglecting  the  first  power  of  the  excentricity, 

:    "  "  "       '     S4=°- 

As  in  Art.  108,  writing  ^  m'C^  for  F,  we  have 
cZR     1      ,dCQ     JP 
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fdR,.          [fdR 
therefore  —  2a  I  -r-  dt  —  on  1 1  -^- 


-  /  sin  {( pn  —  qri)  t  +  Q], 
pn  —  qri  \      aa     pn  —  qri  ' 

where  /  is  an  arbitrary  constant. 
Hence  by  substitution,  we  have 


da 


pn  —  qn       (pn  —  qri) 


115.  This  expression  is  open  to  the  objection  of  contain- 
ing a  term  proportional  to  the  time,  which  being  capable  of 
indefinite  increase,  would  ultimately  vitiate  the  whole  ap- 
proximation. Here,  then,  we  see  the  advantage  of  having 
a  quantity  g  which  may  be  determined  at  pleasure  :  we  will 
so  determine  it  that  the  objectionable  term  shall  vanish. 
This  condition  gives 

1     dC 
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We  may  also  omit  the  constant/,  and  consider  it  as  con- 
tained in  the  epoch.  We  have  then,  writing  for  h  its  value 
no?  VO-  ~~  g2)>  an(^  neglecting  e2,  r  \  * 


no?  [pn  —  qn'      (pn  —  qn')* 


116.  Before  proceeding  to  obtain  the  perturbations  in 
latitude,  we  will  make  a  few  remarks  on  the  forms  of  the  ex- 
pressions for  Sr  and  &0.  If  we  confine  ourselves  to  the  results 
of  the  first  approximation,  it  will  be  seen  on  substituting  the 
value  of  Pj,  that  pn  —  qn  and  n*  —  (pn  —  qn')2  occur  as  divi- 
sors, and  that  the  expression  for  &0  contains  besides,  the 
divisor  (pn  —  qn')*.  The  second  of  these  may  be  written 

{(I  -p)n  +  qn'}  {(I  +  p)n-  qn'}. 
If  then  either 
(i)  pn  —  qn',     (ii)  (1  —  p)  n  -f  qn,     or  (iii)  (1  +p)  n  —  qri, 

be  very  small,  the  corresponding  terms  in  Sr  and  80,  though 
of  a  high  order,  may  yet  be  sensible.  This  is  especially  the 
case  with,  the  first,  since  as  we  have  remarked,  its  square 
occurs  in  the  expression  for  S0.  These  are  instances  of  what 
in  the  preceding  Chapter  have  been  characterised  as  long 
inequalities. 

The  period  of  the  term  P  cos  {(pn  —  qn)  t  +  Q],  which  has 
given  rise  to  these  inequalities,  is 

2"7T 

,;il'  pn  —  qn'  ' 

in  the  case  of  (i),  this  is  very  large,  and  in  that  of  (ii)  or  (iii) 

it  is  very  nearly  equal  to  —  ,  since  pn  —  qn'  is  nearly  equal 
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to  ±  n.  Hence  it  appears  that  terms  in  R  whose  period  is 
either  very  large,  or  nearly  equal  to  that  of  the  planet,  may 
give  rise  to  important  inequalities  in  the  radius  vector  and 
longitude.  Their  actual  importance  will  of  course  depend  in 
part  upon  the  order  of  the  principal  part  of  P  with  respect  to 
the  excentricities  and  inclinations,  i.e.  (see  Art.  50)  upon 
p~q. 

117.     To  integrate  the  equation  for  the  perturbation  in 
latitude. 

The  equation  is  (Art.  107) 


the  position  of  the  plane  of  the  orbit  of  the  disturbed  planet 
at  the  epoch  being  taken  for  the  fixed  plane  of  reference. 

Differentiating  the  expression  for  R  in  Art.  44,  with  respect 
to  z,  we  obtain 


or,  putting  z  equal  to  0,  and  substituting 

a  tan  if  sin  (n't  +  e  —  IT) 
for  z  (see  Art.  42), 

ma  tan  i  sin  (n't  +  e'-XV)  £  Z>0-  \  +  Dl  cos  <£  +  ...}  . 

(2  a  ) 

This  expression,  after  reduction,  consists  of  terms  of  the 
form 

P  sin  {(pn  —  qn)  t  +  Q}, 

where  p  and  q  are  positive  integers,  and  either  may  be  zero. 
Considering  one  such  term,  our  equation  becomes 

^^  +     (rBX)  =  Psin  {(pn  -  qn')  (  +  Q}. 
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Now  as  in  Art.  108, 

^  =  n2  {  1  +  3  e  c  o  s  (nt  +  e  -  OT  )  +  ...}; 

hence,  neglecting  the  product  eBX,  we  have  for  a  first   ap- 
proximation 


+  rf  .  rSX  =  Psin  {(  pn  -  qn')  t+Q}. 


The  integral  of  this  equation  is 

p 

rSX  =  —z  —  -  -  —  sin  {(  pn  -  an')  t  +  Q]  +  A  cos  (nt  -  B). 
n  —  (pn  —  qn) 

If  instead  of  taking  for  the  fixed  plane  of  reference,  the 
plane  of  the  orbit  of  the  disturbed  planet  at  epoch,  we  take  a 
plane  slightly  inclined  to  this,  we  may  omit  the  arbitrary 
term.  For,  denoting  the  planet's  latitude  with  respect  to  this 
plane  by  X,  we  have  approximately 

X  =  tan  i  sin  (nt  +  e  —  II), 

and  it  may  be  shewn  as  in  Art.  110,  that  omitting  the  term 
in  question  is  only  equivalent  to  changing  slightly  the  values 
of  i  and  fl. 


CHAPTER  VIII. 

ON  THE   EFFECTS  WHICH  A  EESISTING   MEDIUM  WOULD 
PRODUCE   IN  THE  MOTIONS  OF  THE  PLANETS. 

118.  IN  the  preceding  Chapters,  we  have  supposed  the 
planetary  motions  to  take  place  in  free  space,  and  the  results 
of  calculations  based  upon  this  hypothesis  manifest  a  very 
close  agreement  with  observation.  There  is,  however,  a 
remarkable  circumstance  connected  with  Encke's  comet  which 
seems  to  indicate  the  possibility  of  the  existence  of  a  very  rare 
medium,  too  rare  indeed  to  cause  any  sensible  resistance  to 
the  motions  of  the  planets,  but  which,  as  we  shall  presently 
see,  may  yet  influence  the  motions  of  comets,  in  consequence 
of  the  extreme  smallness  of  the  masses  of  these  bodies.  It 
has  been  observed  that  the  comet  above  referred  to  (which 
describes  an  elliptic  orbit  in  a  period  of  about  3J  years,)  has 
since  its  appearance  in  1786,  been  moving  round  the  Sun 
with  an  increasing  mean  motion.  Encke  attributes  this  to 
the  resistance  of  a  medium  pervading  space.  We  shall 
therefore  proceed  to  examine  the  effects  which  such  a  medium 
would  produce  upon  the  elements  of  a  planet's  orbit,  assuming, 
in  accordance  with  the  usual  theory,  that  the  resistance 
varies  as  the  product  of  the  density  of  the  medium  and  the 
square  of  the  velocity  of  the  planet.  We  shall  neglect,  in 
the  present  investigation,  all  forces  except  the  Sun's  attrac- 
tion and  the  resistance  of  the  medium;  consequently  the 
planet  may  be  supposed  to  move  wholly  in  one  plane. 
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119.  Let  r,  6  be  the  radius  vector  and  longitude  of  the 
planet,  s  the  length  of  an  arc  of  its  actual  orbit  measured  from 
some  fixed  point  to  its  position  at  time  t,  and  p  the  density  of 
the  medium.  Then  if  k  be  a  constant,  we  may  represent  the 

resistance  on  the  planet  by  Jcp  (-*-] ,  and  the  equations  of 
motion  will  be 


dV        /^#\2_  _  P     7    fds\zdr 

~d?~r(di)  ~~  ~?~  ^  (dt)  ~ds' 

1  cZ  /  2cZ^\_      ,    fds\2  dO 
rJt\dt)~~Kp\dt)rTs' 


7/3 

If  r2  -n=h,  these  may  be  written 


d? 


These   equations   are  the   same   in  form  with  those   of 
Art.  20,  and  may  be  treated  in  a  similar  manner,  —hp-fr  -r. 

taking  the  place  of  -j-,  and  —fy^jfT  tnat   °f  7™-     (See 
Art.  24.)     We  have  from  equation  (2) 

ds 


120.     To  obtain  a  formula  for  calculating  the  mean  dis- 
tance. 

We  might  proceed  as  in  Art.  25,  but  we  shall  here  employ 
the  method  of  Art.  26.     We  have 

<fs_      pdr     ,    fds\z 

^2 ^   ^7o  ""  ^  I  J+  I   ' 


.RESISTING  MEDIUM.  127 

and  by  a  known  formula  of  elliptic  motion 


dt          r       a' 
Differentiating  the  latter,  we  obtain 

9  ds  d*s  _  _  2/i  dr      /A  da 
dt~df~  ~*»  df-cfcf*' 

(fe 

and  multiplying  the  former  by  2  -77  , 

etc 

ds 


therefore 

a"  dt  r  \dtJ 


dt  ~        p     \dtj  ' 


121.    To  oitaire  a  formula  for  calculating  the  eccentricity. 
We  have,  as  in  Art,  27, 


~dt 


,  drds 
Differentiating  as  if  r  were  constant,  and  writing  -  &-^-  ^ 

fr 
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Now  from  equation  (3) 


/dr\2  ,  tf_MV 
(dt)  +r*~  tf 


,        _ 
"  r       h2 


,1       c 
therefore 


u?ede        07    ds  (u,     //,2(1—  e2)) 
^  -=-  =  -  2&p  T.  |-  -  ^—  ^  —  ^  f  ; 
/i2  a^  r  at  (r  h*       j  ' 


ds 


This  -result  may  also  be  obtained  by  differentiating  the 
formula  h*  =  pa  (1  —  e2),  and  substituting  the  expressions  for 

dh       ,  da         .     •   ,   oo 

-^-  and  -j-  ,  as  in  Art.  28. 
at         at 


122.     To  obtain  a  formula  for  calculating  the  longitude  of 
perihelion. 

We  have,  as  in  Art.  29, 


Differentiating  as  if  r  and  0  were  constant,  and  writing 
,    dr  ds  f     d?r 
'^dtdt**  df 

dr          o  //i        ^dtff      ,    dr  ds 

-COSec*(e-*)irt-kPll-tJt 

dh         .   ,  ds 
or  since  = 


RESISTING  MEDIUM.  129 


but  from  equation  (4)  of  Art.  22, 
dr  . 


therefore  =  -sin  (0-  w) 

dt  e  '  dt 

123.     Jb  oZto'ra  a  formula  for  calculating  the  longitude 
of  the  epoch. 

We  have  (see  Art.  13), 

0  =  w+/(w«  +  e-«r,  e). 
Differentiating,  the  elements  being  considered  variable, 


dt       dt      ndt  dt 

and  differentiating  as  if  the  elements  were  invariable, 

d0df 


Equating  the  two  values  of  -y-  ,  we  have 


=(\ 
dt      \dtj' 

-y- 
ut 

dm     1  dO  f  dn     de     dvr\      dfde_() 
df+ndt(   dt+Tt~~~dt)  +  ~aedt~    ' 

dO     h      df     dO         (ft      1\    .    IA       WA        K. 
TT=-a,    -j-^-j-^a  (TZ  +  -    sm(^-w)  (Art.17); 
d<     r      cZe      de         \^      f/ 


But 

at     \ 

.,       r  de         .dn 

therefore  •3i*s  —  t-si-r\  +  —   7    /    ^ 

dt  dt      V         h  J  dt 

nar2 
"HT 

C.  P.  T. 
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As  in  Art.  87,  we  may  omit  the  term  t  -=-,  if  we  bear  in 
mind   that   the   mean   longitude  will  then   be   denoted   by 

ndt  +  e.     Thus,  on  substituting  the  values  of  — =-  and  - 

at  at 

we  obtain 

de  _     2kp   .    fa  (        nr2 

-=-_sm(0-<*)    1-T 


dt  ' 

124.  We  shall  now  express  the  results  of  the  preceding 
Articles  in  a  form  convenient  for  application ;  and,  for  sim- 
plicity, terms  of  the  order  of  the  square  of  the  disturbing 
force  will  be  neglected. 

If  u  denote  the  excentric  anomaly,  we  have 

r—  a  (1  —  ecosu) — , (1), 

nt  +  e  —  OT  =  u  —  e  sin  u (2). 

fds^ 


W 
a  \1  — 


e  cos  u 

/j,  1  +  e  cos  u 
a  1  —  e  cos  u  ' 

„         ds  //I  +  e  cos  u\  l-i 

therefore     -7-  =  na  A  /   r  -  ,  since  wa-  =  y/^- 
rf^          Y  VI  -  e  cos  M/ 

And  to  transform  the  independent  variable  from  t  to  u, 
we  have  from  equation  (2) 

dt      1  ,,  N     1  cZe  ,  1  ^«r 

-7-  =5  -  (1  —  e  COS  U)  ---  j-  +  -  -y    , 

du     n  n  du     n  du 
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of  which  the  two  last  terms,  being  of  the  order  of  the  dis- 
turbing force,  may  be  omitted,  since  the  terms  in  which  the 
substitution  is  to  be  made  are  themselves  of  the  first  order. 

To  obtain  an  expression  for  sin  (0  —  vr),  which  occurs  in 
the  formula  for  the  longitudes  of  perihelion  and  of  the 
epoch,  we  have  from  equation  (1),  omitting  terms  of  the 
order  of  the  disturbing  force, 

dr  .        du 

dt  =  ae&mu^t' 

and  from  the  equation 

£  =  ^{l+ecos(0-w)}     (Art.  22), 
dr        e   . 


Equating  the  two  values  of  -7-  , 

ctt 

.    ,~        x      ha   .       du 
sin  (V  —  vr)  =•  —  sin  u  -y-  . 

fA  tlti 

Hence,   by  substitution,   the  formulae  of  the  preceding 
Articles  become 

da         07     2 /i  \      //l  +  ecoswN 

j—  =  —  zkpa*  (l+e  cos  u)  A  /  I- , 

du  '  V  U-ecosW 


-j~ 
du 


—  e  cos  u 


—  e  cos  U, 
9—2 
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125.  These  formulae  are  sufficient  to  determine  the 
elements  of  the  orbit  at  any  time,  and  being  perfectly  gene- 
ral, are  applicable  as  well  to  the  motion  of  cornets,  as  to 
that  of  planets,  but  before  we  can  integrate  them,  we  shall 
require  a  knowledge  of  the  form  of  p.  Now  the  analogy  of 
the  terrestrial  atmosphere  would  lead  us  to  suppose  that  if 
the  Sun  be  surrounded  by  an  ethereal  medium,  its  density 
decreases  as  the  distance  from  the  Sun  increases.  Moreover, 
the  researches  of  Professor  Encke  on  the  comet  which  bears 
his  name,  seem  to  indicate  the  law  of  the  inverse  square. 
We  will,  however,  merely  assume  p  to  be  such  a  function 

of  r.  that  when  multiplied  by  A  /( - ]  ,  and  deve- 

7  V  \l  —  ecosuj 

loped  in  a  series  of  cosines  of  u  and  its  multiples,  it  takes  the 
form 

A  +  Be  cos  u  +  Ce*  cos  Zu+  ... 
Thus  our  formulae  become 
da 


=  -  2ka2  { A  +  (A  +  B)  e  cos  u  +  ...}, 
-T-  =  —  2Jca  I A  cos  u  +  —  (1  +  cos  2u)  +  . . .  j-  , 


e  -j^  =  —  2ka  \  A  sin  u  +  —  sin  2w  +  ....  [•  , 


~7~~   =  —  AKU,  -\  J±  SILL  U  1 — 

du  2 


-T-  =  ka  (Ae  sin  u  +  . . .). 


126.     Supposing  the  orbit  nearly  circular,  to  examine  the 
effects  of  the  medium  upon  the  elements  of  the  orbit. 

Since  the  orbit  is  nearly  circular,  we  shall  neglect  squares 
and  higher  powers  of  e ;  thus  the  preceding  formulas  give 


on  integration 
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a  =  const.  -  2ka?  {Au  +  (A  +  B)  e  sin  u}, 

(  .              Be  f       sin  2w\) 
e  —  const.  —  2ka  <  A  sin  u  +  -~  (u-\ —  > , 

I  *    \  *        J) 

eer  =  const.  +  2&a  \  A  cos  u  +  -7-  cos  2w  !• , 

=  const.  —  JcaAe  cos  w. 

Hence  in  an  entire  revolution  of  the  planet,  the  mean 
distance  is  diminished  by  4t7rka*A,  and  the  excentricity  by 
%7rkaBe,  while  the  longitudes  of  perihelion,  and  of  the  epoch, 

remain  unchanged.     Also  from  the  formula  n  =  — V  ,  it  ap- 

a2 

pears  that  the  mean  motion  is,  in  an  entire  revolution,  in- 
creased by  QirknaA. 

127.  We  have  already  remarked  that  no  traces  of  a 
resisting  medium  have  yet  been  discovered  in  the  motion  of 
the  planets :  but,  since  k  varies  inversely  as  the  mass  of  the 
body  acted  upon,  the  formulas  of  Art.  124  shew  that  such 
a  medium,  though  too  rare  to  influence  the  planets,  might  yet 
sensibly  affect  the  motions  of  comets,  in  consequence  of  the 
extreme  smallness  of  their  masses. 


PBOBLEMS. 


1.  SUPPOSING  in  the  Problem  of  the  Three  Bodies  the 
relative  orbit  of  two  of  the  bodies  to  be  a  circle  described 
uniformly,  obtain  equations  for  determining  the  motion  of  the 
third  body ;  and  transform  the  system  of  co-ordinates,  so  that 
the  plane  of  the  circular  orbit  being  that  of  xy,  the  axis  of  x 
shall  always  pass  through  the  two  bodies  in  that  plane. 

2.  Shew  that  the  plane  of  the  orbit  of  a  planet  revolves 
about  the  planet's  radius  vector  as  an  instantaneous  axis*. 

8.  A  particle  is  describing  an  orbit  round  a  centre  of 
force  which  is  any  function  of  the  distance,  and  is  acted  upon 
by  a  disturbing  force  which  is  always  perpendicular  to  the 
plane  of  the  instantaneous  orbit,  and  inversely  proportional 
to  the  distance  of  the  body  from  the  centre  of  the  principal 
force.  Prove  that  the  plane  of  the  instantaneous  orbit  re- 
volves uniformly  round  its  instantaneous  axis. 

4.  A  particle,  acted  on  by  a  force  varying  as  the  inverse 
square  of  the  distance  and  always  tending  to  a  fixed  point, 
suffers  slight  disturbances :  prove  that  generally  there  is  a 
conic  with  the  fixed  point  for  focus,  with  which  the  body's 
actual  path  has  a  closer  contact  than  with  the  instantaneous 
conic. 

*  In  this  and  the  following  problem,  the  plane  of  the  orbit  must  be 
supposed  to  have  no  angular  velocity  about  a  normal  to  itself.  See  note  to 
Art.  19. 
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5.  Find  when  the  curvature  of  the  instantaneous  orbit 
of  a  body,  acted  on  by  disturbing  forces,  is  the  same  as  that 
of  the  actual  orbit  ;  and  shew  that  this  is  always  the  case 
when  the  only  disturbing  force  arises  from  the  action  of  a 
resisting  medium. 

6.  If  R  be  expressed  on  the  one  hand  as  a  function  of 
rx,  0l9  and  z  (Art.  11),  and  on  the  other  as  a  function  of  r, 
6,  i,  and  O,  6  being  measured  on  the  plane  of  the  orbit  from 
the  node,  prove  that 

dR_dR 
det~d&' 

and  obtain  a  formula  for  calculating  the  inclination. 

7.  If  R  be  expressed  as  a  function  of  t  and  the  usual 
elements,  obtain  the  formula} 

dtl         1     dR     di     cot0 


dt      h  sin  i  di  '   dt  ~~          h          di  ' 

where  6  is  measured  on  the  plane  of  reference  as  far  as  the 
node,  and  thence  on  that  of  the  orbit,  and 


n 

dt 


0    .  2  idfl\ 
—  2  sm'rj-jr   . 
2  dt) 


8.     The  central  force  being  ^  +  ^,  obtain  the  following 
equation  for  the  apsidal  motion, 


dt  ~        e*J/j,  r" 

a,  e  and  to-  being  elements  of  the  instantaneous  ellipse. 

9.     If  the  central  force  on  a  planet  be  ^  +  //,  shew  that 
the  planet  may  be  supposed  at  any  instant  to  be  moving  in  an 
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ellipse,  of  which  the  mean  distance  is  inversely  proportional 

to  the  planet's   distance   from   the    Sun,   provided  — ,  =  2pa, 

where  p  and  a  are  two  simultaneous  values  of  r  and  a.  Shew 
also  that  the  latus  rectum  of  this  ellipse  is  constant. 

10.  When  the  disturbing  function  R  is  independent  of  6, 

/»    -i  «    de       ,  dixr 

find  expressions  for  -j-  and  —rr  . 
at          at 

If  R  —  — ,  these  expressions  give  variable  values  for  e 

and  ts-,  whereas  the  motion  of  the  body  actually  takes  place 
in  a  fixed  ellipse:  shew,  this,  and  explain  the  apparent 
paradox. 

11.  A  planet  describes  an  orbit  under  the  action  of  a 
force  ^  tending  to   the   S.un,  p  not  being  quite  constant : 

obtain  the  following  equations  for  the  variations  of  the  ex- 
centricity  and  longitude  of  perihelion : 

=  —  cos  (6  —  -or), 


"GT  -       //i  \ 

/j,e  -T-  =  —  sin  (6  —  OT). 

If  d/ju  be  always  positive,  what  in  a  whole  revolution  is 
the  nature  of  its  effects  upon  the  excentricity  and  position  of 
the  major  axis? 

12.  Prove  that  for  any  closed  central  orbit,  the  element 
introduced  in  the  equation  of  vis  viva  is  subject  to  periodic 
variations  only,  when  the  disturbing  force  is  due  to  forces 
tending  to  centres  and  functions  of  the  distances  from  those 
centres. 
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13.     If  the  equation  of  the  Moon's  orbit  be  reduced  to 
the  form 

tfu 


shew  that  the  excentricity  and  longitude  of  perihelion  may 
be  found  from  the  equations 

~  =  -/sin  (6  -  «r),     e-^  =/cos  (0  -  w). 

Apply  these  equations  to  find  e  and  w,  when  /  is  a  small 
disturbing  force,  depending  only  upon  the  Moon's  distance 
from  the  Earth. 

14  Assuming  the  differential  equation  for  s  in  the  Lunar 
Theory  to  be 

d*s  (33  ) 

—  +  s  =  -  m*s  \  -  +  ^  cos  2  (6  -  m6)  I 


shew  that  if  7  be  the  longitude  of  the  Moon's  node, 


0       m*  t1  ~  cos  2  (m^~7)  ~  coa  2  (^  ~  'y) 

+  cos  2(0  -w0)}. 

^ 
From  the  above  expression  for  ~  ,  find  the  ratio  of  the 

mean  motion  of  the  node  to  that  of  the  Moon,  taking  into 
account  terms  of  the  order  m4. 

15.  If  two  planets  disturbing  one  another  were  revolving 
in  periods  of  350  and  201  days,  what  form  of  terms  in  the 
disturbing  function  would  demand  examination  ? 

16.  The  periods  of  Venus  and  the  Earth  are  2247  and 
365*256  days  respectively;   find  approximately  the  period  of 
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the  long  inequality  arising  from  their  mutual  perturbations, 
the  important  term  in  the  disturbing  function  11  being  of  the 
form 

PeV2  cos  {13  (nt  +  e)  -  8  (n't  +  e')  -  3w  -  2«r'}. 

17.  The  radius  vector  of  a  planet  is  affected  with  a  small 
periodical  inequality;  shew  that  its  effect  may  be  represented 
by  continued  and  periodical  alterations  of  the  excentricity  and 

PT 

longitude  of  perihelion,  the  period  of  either  being    p— —„, 

where  P  is  the  period  of  the  planet  and  T  that  of  the  in- 
equality. 

18.  If  in  addition  to  the  force  of  the  Sun  on  a  planet 
there  be  a  small  force  tending  towards  the  Sun,  and  varying 
inversely  as   the   mth   power  of  the  distance  of  the  planet 
from  the  Sun,  prove  that  the  perihelion  of  the  orbit  will  have 
a  progressive  or  regressive  motion  according  as  m  is  greater 
or  less  than  2. 

Can  you  explain  this  result  by  reasoning  similar  to  that 
used  in  Airy's  Gravitation? 

19.  It  has  been  found  by  comparing  theory  with   ob- 
servation that  the  perihelion  of  Mercury  progresses  at  a  rate 
greater  by  a  than  that  due  to  the  attraction  of  known  bodies  : 
shew  that  this  increment  would  be  accounted  for  if  the  law  of 

force  tending  to  the  Sun  were  -g  +  -4,  and  if  //  =  ac4  A/->  the 

orbit  being  supposed  to  be  nearly  a  circle,  and  the  mean  dis- 
tance to  be  c. 

20.  The  central  force  acting  on  a  body  being 

w, 
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shew  to  terms  inclusive  of  /*'  and  the  square  of  the  excen- 
tricity,  that  the  motion  is  in  an  ellipse  revolving  uniformly 
about  the  focus. 

21.  Shew  by  means  of  the  formula 

da  _  2na?  dR 
~dt~     fju     de' 

that  the  chief  perturbation  of  the  axis  major  of  the  Moon's 
orbit  may  be  expressed  by  the  equation 

a=a\l  +  0     ,  H  —  ^cos 
(        2n  (n  —  n  ) 

where  n  and  n  are  the  mean  motions  of  the  Moon  and  Sun 
respectively. 

22.  A  satellite  revolving  in  an  ellipse  of  small  excen- 
tricity  is  disturbed  by  another  satellite  revolving  about  the 
same  primary;   find  approximately  the  variation  of  the  mean 
distance  and  the  motion  of  the  apse,  corresponding  to  the 
terms 

ri* 

--  r~  [1  +  3  cos  {2  (n  -  n)  t  -f  e  —  e'J] 

TD 

in  the  function  R,  having  given 
dR       -dsr 


dt  ~~    //,     de  '       dt  fie          de  ' 

23.  Prove  that,  neglecting  periodical  variations,  the  ex- 
centricity  of  any  orbit  can   always   be   represented  by  the 
diagonal  of  a  parallelogram,  whose  sides  are  constant,  and 
angle  varies  uniformly. 

24,  Given  the  equations 

tan2  i  =  N*  +  N22  +  2N&  cos  (hf  +  Sx  -  8a), 


tan  II  =      _ 

N  cos  A    +  S    +       cos  a1 
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explain  the  nature  of  the  motion  of  the  node,  when  the  mini- 
mum inclination  is  zero. 

25.  Prove  that  as  far  as  secular  variations  only  are  con- 
cerned the  function  F  is  constant. 

26.  Considering  only  secular  variations,  obtain  the  fol- 
lowing equations : 

^A    /  /         y  \A/UJ  \  y^-y          r*  /  .       .         Q   •  CL\L  \  ^M 

2,   m  v  a  e  -7-  1  =  C,     Z  im  A/a  tan  i  -7-  1  =  0. 
\  cic  /  \  dt  J 

27.  If  the  squares  of  the  masses  of  two  mutually  dis- 
turbing planets  were  to  each  other  inversely  as  their  mean 
distances,  shew  that  the  nodes  would  oscillate  through  equal 
angles. 

28.  If  Mt  m,  m'  be  the  masses  of  three  bodies  mutually 
attracting  according  to  the  law  of  gravity,  M  being  much 
larger  than  m  or  m',  and  if  v,  v  be  the  velocities  of  m,  m  at 
distances  r,  r'  from  the  centre  of  M ,  supposed  fixed,  shew  that 
the  equation  of  vis  viva  for  this  case  may  be  assumed  to  be 

2        '  /2  ,  oTM-fm      m  ,  m>      m'\ 

mv  +  m  v  +  21T  5 h  ^— , r  }  =  0, 

\2a      r      2a      r  J 

2a  and  2a'  being  the  major  axes  of  the  instantaneous  ellipses 
of  m  and  m. 

29.  Infer  from  the  foregoing  equation  by  the  method  of 
the  variation  of  parameters  the  ratio  of  simultaneous  changes 
in  the   mean  distances  and  mean  motions  of  two  planets 
mutually  disturbing. 

30.  If  r  be  the  true  radius  vector,  6^  the  projected  lon- 
gitude, and  X  the  latitude  of  a  planet,  obtain  the  following 
equation  of  motion : 

fZV  ,     (dO^        fd\\*  .f*_dR 

-i  .1  "~~  /  COS    fo  1     -i     I    *"""  /I     -f     I    "T"     .>  -~*     t     * 

dt  \dtJ         \dtj      r'      dr 
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31.  Obtain  the  following  equation  between  the  pertur- 
bations of  a  planet  in  longitude  and  radius  vector,  whatever 
be  the  law  of  force,  provided  it  be  central  and  a  function  of 
the  distance  only,  and  provided  such  a  function  as  R  can  be 
found : 

.  dR 


rdt'          *}    dt  "dr 


dr 

where  F  denotes  the  central  force,  and  h  twice  the  sectorial 
area  described  by  the  undisturbed  planet  round  the  Sun. 

32.  If  the  orbits  of  two  planets  which  disturb  each  other 
be  very  nearly  circular,  shew  that  the  inequalities  of  the 
radius  vector  may  be  immediately  deduced  from  those  of  the 
longitude  by  means  of  the  equation 


•"•    dt        la     n      73  _  n 

~h 


dt 
33.     Integrate  the  equation 

^M  +  n*  .  rSr  =  5  {P  cos  (pnt  +  Q)}, 
determining  the   arbitrary  constants    so    that   Sr  =  0,  and 

7        ^ 

—  -,  —  =  0,  when  t  =  0  :  and  shew  that  for  small  values  of  t, 
at 


the  case  of  p  =  1  being  included. 
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34.     A  planet  moves  in  a  resisting  medium  of  which  the 
resistance 


apply  the  equation 


.f*P_2a_         o 

r3  J    dt  dr 

to  obtain  the  following,  in  which  e*  is  neglected  : 

*        *  t  rgr  .  3g  cos  wf  +  e  -  w 


*  sin  /rt  +  e  -  OT     =  0. 


35.  The  co-ordinates  of  the  position  at  any  time  t  of  a 
disturbed  planet  being  #+&c,  y  +  Sy,  z+  Sz,  reckoned  from 
the  Sun's  centre  as  a  fixed  origin,  and  referred  to  the  plane 
of  motion  at  a  given  epoch  ;  and  r  being  the  heliocentric 
distance,  x,  y  the  co-ordinates  of  the  position  which  the 
planet  would  have  had  at  the  time  t,  if  the  disturbance  had 
ceased  at  the  given  epoch  ;  obtain  the  following  equations  for 
determining  &P,  By,  Sz  to  the  first  order  of  the  disturbing 
force  : 


If-/  I/        \      f\  *JWUI        fy          I  UVJL.W  f^ 

dt>    -r..  --4)%--/  a^K^r=o, 


d  rA        dz         ' 

in  which  JJL  is  the  sum  of  the  masses  of  the  Sun  and  planet, 
and  R'  is  put  for 

-  TO-  (*-  +  y"  +  /•)-* + »'  {(*  -  *J  +  (y  -  y'T  +  0  -  *')"!  -  *, 

m'  being  the  mass,  and  a;',  yr,  ^'  the  heliocentric  co-ordinates 
of  the  disturbing  planet. 
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36.  Shew  that  the  effect  of  a  resisting  medium  on  the 
instantaneous  orbit  of  a  planet,  would  be  to  make  the  apsidal 
line  regrede  or  progrede,  according  as  the  planet  moved  from 
perihelion  to  aphelion,  or  from  aphelion  to  perihelion. 

37.  Two  small  planets  P,  Q,  very  near  each  other,  re- 
volve about  the  Sun  in  orbits  very  nearly  circular,  and  make 
two  revolutions  about  each  other  while  they  make  one  revo- 
lution about  the  Sun.     Compare  the  sum  of  their  masses  with 
the  mass  of  the  Sun. 

If  the  line  PQ  move  parallel  to  itself,  what  inference  do 
you  draw? 

38.  If  the  motion  of  a  planet  round  the  Sun  be  disturbed 
by  the  action  of  another  planet,  the  latter  being  supposed  to 
describe  a  circular  orbit  of  radius  a!  with  uniform  velocity  n, 
obtain  the  following  exact  equation  : 


dt  n 

2/A     2mV  2m' 

= — nr~  cos  m  ^ rT> 

f        «  (a 2  —  2a  r  cos  ft>  +  r  )2 

where  r  is  the  radius  vector  of  the  disturbed  planet,  rlt  0lf  z 
its  co-ordinates  referred  to  a  fixed  plane,  and  &>  the  inclina- 
tion of  the  radii  vectores  of  the  disturbed  and  disturbing 
planets  to  each  other. 

39.  Prove  that,  if  the  periodic  times  of  a  disturbed  and 
disturbing  planet  are  not  commensurable,  the  secular  changes 
of  the  orbit  of  the  disturbed  planet  are  the  same  as  they 
would  be  if  the  mass  of  the  disturbing  planet  were  distributed 
over  its  orbit,  in  such  a  manner,  that  the  part  of  the  mass 
distributed  over  each  portion  of  the  orbit  should  be  propor- 
tional to  the  time  which  the  planet  actually  takes  to  describe 
that  portion. 
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Additional  Examples. 

40.     If  the  disturbing  force  be  radial  and  equal  to  //r, 
prove  that  the  mean  distance  at  any  time  is 

a  —  2a  —  \e  cos  (nt  +  e  —  w)  +  1-  cos  2  (nt  4-  e  —  «•)  k 

,  1872.) 


41.  A  particle  moves  in  a  plane  subject  to  a  force  directed 
to  a  fixed  point  in  the  plane,  and  depending  on  the  distance 
from  the  point  :  prove  these  equations  : 

tfde     hdb__dR 
a3  dt  +  6*  dt~      da' 

h2  de      h  da  _     dR 
T3~di~tfTt~  ~~db9 

de      ,        .  da 


where  a  and  b  are  the  intercepts  of  the  direction  of  motion  on 
fixed  axes  through  the  fixed  point,  h  is  an  absolute  constant, 
E  is  the  potential  of  the  force  and  e  a  quantity  such  that,  v 
being  the  velocity  of  the  particle  at  time  £,  v  (t  +  e)  is  its 
distance  measured  along  the  direction  of  motion  from  the 
axis  of  x.  (M.  T.  1873.) 

42.  A  large  number  of  meteors  are  distributed  uni- 
formly in  the  circumference  of  a  circle  of  radius  b  concentric 
with  the  sun  and  in  the  plane  of  the  planet's  orbit  ;  shew 
that  the  perturbation  in  the  mean  distance  of  the  planet  is 
given  by 


-    - 
A6     \ 

M  being  the  whole  mass  of  the  meteors,  a,  e,  n,  e,  w  being 
elements  of  the  planet's  orbit,  and  b3  :  a3,  I/2  :  /i.2,  and  e2  being 
neglected.  (M.  T.  1880.) 
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43.  If  a,  I  be  the  semiaxes  of  the  instantaneous  ellipse 
of  a  planet's  orbit,  and  the  only  disturbing  force  be  a 
tangential  one  whose  acceleration  is/,  prove  that 


w  dt  dt 

(M.  T.  1881.) 

44.  If  H  be  proportional  to  the  excess  of  the  moment 
of  momentum  of  the  orbital  motion  of  a  disturbed  planet 
about  the  Sun  in  a  circular  orbit  at  distance  a  above  the 
value  of  the  same  thing  in  the  actual  elliptic  orbit ;  then 
shew,  from  the  equations  for  the  variations  of  major-axis 

and   excentricity,  that,  k  being  a  constant,  ~-=-  =  —  k  -j- . 

dt  dw 

Also  assuming  that,  when  the  inclination  of  the  orbit  to  the 
plane  of  reference  is  zero,  the  variational  equation  for  -&  is 

dvr  _  na  Jl  —  ez  dR 
dt  ae         de  ' 


,        ,,  dv         dR 

Shew  that  -7  -  =  k  -yyy-  . 

dt         dn. 


(M.  T.  1882.) 


45.  Explain  how  Elliptic  Integrals  and  Bessel's  Func- 
tions arise  in  the  development  of  the  disturbing  function. 

(M.  T.  1879.) 

46.  If  m,  m   denote  the  masses  of  two  planets,  n,  ri, 
their  mean  motions,  and  a,  a,  their  mean  distances  from  the 
Sun,  and  if 

m'  —^  C  cos  {(in  —  in)  t  +  7], 


m—zC'  cos  {(i'n1  —  in)  t  +  7}, 
a 

C.  P.  T.  10 
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represent   any  two   corresponding  terms   in   the   respective 
developments  of  the  portions 

r  rr 

m'  -72  cos  (r,  r'}  and  ra  -,  cos  (r,  r'), 


of  the  two  disturbing  functions,  prove  that 


(Smith's  Prizes,  1880.) 

47.  Shew  that  the  use  of  the  ordinary  disturbing  func- 
tion is  attended  with  certain  inconveniences  in  the  case 
where  we  have  to  determine  the  perturbations  of  a  superior 
planet  produced  by  the  action  of  an  inferior  one.  Also 
shew  that  these  inconveniences  may  be  avoided  by  referring 
the  motion  of  the  superior  planet  to  the  common  centre  of 
gravity  of  the  Sun  and  inferior  planet  as  origin,  and  by  taking 

a  quantity  proportional  to  -^  +  -  —  -  —  —  as  the  disturbing 

-t*>         p  f 

function,  where  3/is  the  Sun's  mass,  p  that  of  the  disturbing 
planet,  and  R,  p,  and  r  are  the  distances  of  the  disturbed  planet 
from  the  Sun,  from  the  disturbing  planet,  and  from  the  centre 
of  gravity  of  these  two  bodies  respectively.  How  would  you 
develope  the  above  disturbing  function?  (8.  P.  1880.) 


48.  Considering  the  secular  variations  of  the  excentrici- 
ties  and  perihelion  longitudes  of  two  mutually  disturbing 
planets  ;  prove  that  in  the  case  where  m  Ja  =  m'  Ja,  e  and  e 
have  the  same  greatest  and  least  values  E^  and  E2,  and  -or  ~  -sr' 

.ft  2  _    77»  2v 

is  never  greater  than  tan"1  1  —  ^-=-=^  J  . 

(M.  T.  1874.) 
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49.  If  there  were  only  two  planets  moving  about  the 
Sun  and  at  any  instant  the  orbit  of  one  were  circular,  the 
eccentricity  of  its  orbit  would  always  be  given  by 

c  .  cos  (at  +  /3) 

where  c,  a,  0  are  constants,  and  the  longitude  of  its  peri- 
helion would  increase  uniformly  with  the  time. 

(M.  T.  1876.) 

50.  If  7  be  the  angle  between  the  orbits  of  two  mutually 
disturbing  planets,  prove  that  approximately, 

/       2  ,      *  AT  :  *  ,  4mm'     aa  .  sin2 
m  >Ja.e  -}-  m  »Ja  .  e  + 


m 

is  constant.     (Pontecoulant,  Systtme  du  Monde,  Tome  I.  pp. 
458—460.) 

51.  If  I  be  the  mean  anomaly,  v  the  true  anomaly,  g  the 
angular  distance  of  the  perihelion  from  the  ascending  node, 
h  the  longitude  of  the  ascending  node  in  the  fixed  plane  of 
reference, 


-e2)}    H=     ^a  (1  -  e2)  cos  *, 
the  other  letters  having  their  usual  meanings,  and  if  R  denote 

the  ordinary  disturbing  function  plus  the  term  -j~,  prove 

Ziu 

that  x,  y,  zy  the  co-ordinates  of  the  planet  referred  to  the  fixed 
plane  as  the  plane  of  xy  and  to  the  line  from  which  the 
longitudes  are  measured  as  the  axis  of  x,  are  given  by  the 
equations 

Gcc  =  Gr  cos  (v  +  g)  cos  h  —  Hr  sin  (v  +  g)  sin  h, 
Gy  —  Gr  cos  (v  +  g)  sin  h  +  Hr  sin  (v  +  g)  cos  A, 


10—2 
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and  taking  I,  g,  h,  L,  G,  H,  as  the  six  elements,  prove  that 


^  =  dR     dH  _dE 

dt  ~  dl  y    dt  ~~  dg  '     dt~dh' 

(M.  T.  1878.) 

52.  If  M,  m  be  the  masses  of  Sun  and  disturbed  planet, 
H  the  disturbing  function  due  to  the  action  of  another  planet; 
•G7,  v  the  longitudes  of  the  perihelion  and  ascending  node,  e 
the  time  of  passage  through  perihelion  ;  a,  e,  i,  the  semiaxis 
major,  excentricity  and  inclination  of  the  orbit  to  a  fixed 
plane  ;  further  if 

M+m 
^  --  toT> 

K  =  *jM+m  Ja  (1  -  6*), 

X  =  jM+m.  Ja  (1  —  e2)  (1  -  cos  i)  : 

then  the  pairs  of  elements  (e,  //),   (K,  «r),   (v,   X)    will  be 
canonical,  i.e. 


de  _  (JO     dp  _  _  dto     o 
dt  ~  dJL '    dt~    ~  de  ' 


(8.  P.  1881.; 


APPENDIX. 

ON  THE  FORM  OF  THE  EQUATIONS  OF  ART.  39. 

1.  ON  referring  to  Art.   39,  it  will   be   seen  that  the 
formulae  which  have  been  obtained  for  calculating  the  elements 
of  the  orbit  involve  only  partial  differential  coefficients  of  R 
with  respect  to  these  elements,  multiplied  by  functions  of  the 
elements  which  do  not  contain  the  time  explicitly.     As  it  is 
to  this  circumstance  that  these  formulas  mainly  owe  their 
advantage  (since  it  renders  them  fit  for  use  as  soon  as  the 
partial  differential  coefficients  have  been  calculated),  it  will 
be  interesting  to  shew  d  priori  that   whatever  system   of 
elements  be  adopted,  the  formulae  for  their  calculation  can 
always  be  arranged  in  this  way. 

2.  If  the  motion  of  the  planet  be  referred  to  three  rect- 
angular axes  originating  in  the  centre  of  gravity  of  the  Sun, 
we  have  the  equations  of  motion  (see  Art.  9), 

d?x  ,  fjix  _  dR  ,  . 

^~         ' 


dR 

*  "s* 


tfz      iz      dR 
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Let  a,  b,  c,  d,  e,f  be  the  six  elements  introduced  by  in- 
tegrating these  equations  when  R  =  0,  and  for  ~  ,   ~ ,    -^ 

etc      at      at 

write  a?',  #',  z:  then  a?',  ;?/'  and  zf  can  be  expressed  as  functions 
of  t  and  the  elements ;  hence 

dx'  _  /dx\      dx  da     dx  db 
'dt~\dt)~da'didbdi  + ' 

where  in  f-g-J  the  elements  are  supposed  constant. 
If  in  equation  (1)  we  put  R  equal  to  0,  we  have 


therefore 


therefore  d^  =  d^  di+db  Ji 


and  similar  equations  hold  for  -=-  and  ^-  . 

Now  since  H  is  a  function  of  #,  ^  an(^  ^> 

dE^_dRdx     dRdy     dfidz 
,.   ;        da      c?a;  c?a      dy  da     dz  da 

_  /dx  dx'     dy  dy      dz  dz'\  da 
\da  da     da  da      da  da)  dt 

(dx  dx      dy  dy'     dz  dz\  db 

t     ( [     <7_       <7       I     j    

\da  db      da  db      da  db)  dt 

(dx  dx      dy  dy      dz_  dz'\  dc 
\da  dc     da  dc      da  dc)  dt 
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3.  We  may  eliminate  -^-  from  this  expression:  for,  sup- 
posing x,  y,  and  z  expressed  as  functions  of  t  and  the  elements, 
we  have 

dx  _  /dx\      dx  da     dx  db 

di  =  \dt)+da'dt+'db  dt  +  '"' 

but  by  the  principles  of  the  method   of  the   Variation   of 
Parameters 

dx      fdx 

~di' 

dx  da     dx  db     dx  dc 

therefore  -r-  -y-  +  -^  -,-  +  -j-  -j-  +  ...  =  0. 

da  dt      do  dt      dc  dt 

0.    .,    ,       dy  da     dy  db     dy  dc  - 

Similarly,    -f-  -j-  +  -^  -j-  +  -/-  -j-  +  ...  =  0, 
da  dt      db  dt     dc  dt 

dz  da     dz  db      dz  dc,  _ 

Multiplying  these  equations   by  —,-  ,    s    ,  &c.,  and  adding, 

we  obtain 

_  fdx  dx   dy  dy   dz  dz'\  da 
\da  da   da  da   da  da)  dt 


da   db  da   do  da)  at 

'dx  dx   dy  dy   dz  dz'\  dc 
+  \dc~da+'ac'd^  +  'dcda)dt 


If  this  expression  be  subtracted  from  that  for  -7-  in  Art,  2, 
the  latter  may  be  written 

dR  , ,  db  dc 
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,         r     7-1      dxdx      dxdx'     dydy      dydy 

where  [a,  b]  =  -y-  -y, 77  -7-  +  T  ~w  —  ~?i  •? 

L      J     da  d&      d6  da     da  d6      db  da 


.  i      . 

da  db      db  da  ' 


dz  dz'     dz  dz 

Similarly, 

dR     „     -.  da  .         .,  dc 


4.  By  successive  elimination  between  these  equations, 

we    can   obtain   expressions   for  -y-,    -y-  ,  &c.,  in  terms   of 

ctt      ctt 

-y-  ,  -y=-,  &c.,  [a,  Z>],  [a,  c],  &c.r   if,  then,  we  can  shew  that 
da     do 

[a,  5],  [a,  c],  &c.,  are  independent  of  the  time  explicitly,  it  will 

follow  that  this  is  also  the  case  with  the  coefficients  of  -y-  , 

da 

dR    p      .     ,,  .        «     da   db    p 

-^-,  &c..  in  the  expressions  for  -y-  ,  -y-  ,  &c. 
a6  a^    a^ 

5.  To  shew  that  [a,  b]  is  independent  of  the  time  ex- 
plicitly. 

Let  V  —  -  ;  then  the  equations  of  motion  give 


dt        dxy      \dt        dy'      \dt        dz' 

Now  differentiating  with  respect  to  t  only  so  far  as  it 
occurs  explicitly, 


d  r     ,  -,  _  dx  d  fdx'\      dx  d  fdx\ 
di La'    *~dadt (db)  +~dbdt (da) 

dx  d  (dx'\      dx  d  fdx\ 
~  dbdt(da)  ~~^adt  (db) 
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dx  d  fdx'\  dx'  dx 
dadb\dt)~db'da 
dx  d  /dx\  dx  dxf 
db  da  \dt)  da  db 


_____ 

~  dadb\dx)     dbda\dx) 

dyd_  fdV\  _  dy_d^  /dV\ 

dadb\dy)     dbda\dy) 

dz_  d^  fdV\  _dzd_(dV\ 

dadb\dz)     dbda\dz) 
_dxd(dV\      dy-d_fdV\      dzd^fd 
~dad~x\db)  +  dad   (db)  +  ~da~dz\db 


dady 

^^     _dyd_  ( 
db  dx\da   db  dy  \  da   db  dz  \  da 

_  d?V_   d*V  = 

da  db   db  da 

Hence  [a,  b]  does  not  contain  the  time  explicitly.  The 
same  is  of  course  true  of  [a,  c],  [b,  c],  &c.  It  follows,  then, 
that  whatever  system  of  elements  be  adopted,  we  can  always 
express  their  differential  coefficients  in  terms  of  the  partial 
differential  coefficients  of  R  with  respect  to  them,  multiplied 
by  functions  of  the  elements  which  do  not  involve  the  time 
explicitly. 

COKOLLARY.    We  now  see  that  in  the  expression  for  [a,  b], 
viz.  :  — 

(dx  dx  __  dx  dx'\      (dy  dtf__dy  dy'\      fdz^  dz  __  (h  dz\ 
\da~ab~~db~da~)  +  \da~db  ~db  ~da)  +  (da~db     dbda)' 
we  may  replace  x,  y,  z,  x,  y'}  z',  by  a?0,  ya,  z0,  aJ0',  y0',  <,  their 
values  at  the  origin  of  the  time,  without  altering  the  value  of 
this  expression.    We  may  next  adopt  instead  of  the  constants 
a,  b,  c,  d,  e,  /,  the  same  quantities  xQ,y^z^  #0',  y0',  £0'.     In 
this  case  we  have  evidently 
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fo,<]  =  i,  [y0,2/o']  =  i,  [*.,*.>!, 

[x0f,  O  =  -  1,  [y0',  yj  =  -  1,  [<,  *J  =  -  1, 

and  all  the  analogous  expressions  different  from  these  will 
vanish.  The  final  equation  of  Art.  3  of  this  Appendix  will 
give 

dxdE          ddR          dz      dR 


=_  =  _  =  _ 

dt  ~  dx0*  dt  ~  dyQ'  dt  dz0' 
and  no  other  combinations  of  the  arbitrary  constants  will  give 
simpler  expressions  than  these  for  the  variations  of  the 
.constants  adopted.  From  them  the  variations  of  the  usual 
elliptic  elements  might  be  derived  by  proper  analytical  trans- 
formations. The  group  of  constants  x0,  y0>  ZQ>  &0',  yd,  ZQ  may 
be  termed  canonical  elements. 

6.  From  the  formula  of  Art.  3  of  this  Appendix,  which 
is  due  to  Lagrange,  those  of  Chapter  II.  may  be  deduced  :  for 
this  we  refer  to  Pontecoulant's  Systeme  du  Monde,  Tome  I. 
p.  542. 

Thus,  if  we  denote  by  u  the  excentric  anomaly  of  the 
planet,  by  nt  +  I  its  mean  longitude  at  the  time  t,  by  e  the 
excentricity  of  the  orbit,  and  by  a  the  semi-major  axis,  con- 
nected with  the  constant  n  by  the  equation  n*a*  =  /-t,  we  shall 
have  by  the  formulas  of  elliptic  motion 
nt  +  l  =  u  —  esm  u. 

Next,  if  we  denote  by  X,   F,  the  two  rectangular  co- 
ordinates of  m  referred  to  the  plane  of  the  orbit,  the  major- 
axis  and  latus-rectum  being  the  co-ordinate  axes,  and  the 
origin  consequently  at  the  focus  of  the  curve,  we  shall  have 
—  ae  Y=al  —  e2  sin  u. 


Lastly,  if  we  denote  by  v  the  angle  which  the  major-axis 
of  the  orbit  makes  with  the  node  or  line  of  intersection  of  the 
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plane  of  this  orbit  and  the  fixed  plane  of  xy,  by  ft  the  angle 
which  the  node  makes  with  the  axis  of  x,  and  by  i  the  in- 
clination of  the  orbit  to  the  plane  of  xy  ;  we  shall  then  have 
by  transformation  of  co-ordinates, 

x  =  (Xcos  v  —  Ysm  v)  cos  ft  —  (X  sin  v  +  Ycos  v)  cos  i  .  sin  ft, 
y  =  (Xcos  v  —  Tsin  v)  sin  II  +  (JTsin  v  +  Ycos  v)  cos  i  .  cos  ft, 
z  =  (X  sin  v  +  Ycos  v)  sin  i. 

These  values  of  x,  y,  z  are  very  convenient  for  use  in  de- 
termining the  values  of  the  expressions  [a,  6],  [a,  c],  [5,  c],  &c., 
since  the  constants  in  the  expressions  for  x,  y,  z,  are  separated 
into  two  groups  ;  the  variables  X,  Y  only  containing  a,  e,  I, 
which  depend  on  the  form  of  the  orbit  and  on  the  position  of 
the  planet  on  that  curve  at  a  given  instant,  being  independent 
of  the  three  constants  v,  ft,  i,  which  fix  the  position  of  the 
major-axis  and  of  the  plane  of  the  orbit. 

By  the  aid  of  these  values  and  of  their  partial  differential 

coefficients  substituted  in  forms  such  as 

fdx  dx      dx  dx\      (dy  dy  __  dy  dy'\      fdz^  d^_^2  dz\ 
\dallb  "~  db  da)      \da  lib  ~~  db  ~da)      \da  db      db  da) 

we  shall  determine  the  values  of  the  fifteen  symbols  [a,  Z], 

[a,  e],  &c.,  and  by  an  easy  calculation  obtain  the  following 

values  : 

r     -n         an     r      n         s/afl-e2)     r      -,  ae 

M]=  -y,  [«,.]-        4^,  M-.^^ 

cos  i  Ja(l  -  e1)  ae  cos  i 


p,  ft]  =  sin  $     a  (1  -  e2). 

The  values  of  the  nine  other  symbols  [a,  e],  [a,  i],  [e,  I], 
[e,  ft],  [I,  i],  [I,  ft],  [ft,  v],  [I,  v],  [i,  v]  are  equal  to  zero. 

After  substituting  these  values  in  the  general  formula  at 
the  end  of  Art.  3  of  this  Appendix  we  shall  find 


156  PLANETARY  THEORY. 

dR  _     andl  _Ja(l-f)dv     cos  ija  (1  -  e2)  dtl 
da~   ~~2di  2a        ~di~  2a         ~  ~dt  ' 

dR  ae        dv        ae  cos  i    d£l 


dR  __  an  da 
~dl=~2~dt' 

dR  _Ja(l-e2)da  ae        de 

~dv~     ~2a        'dt  ~  Ja(l^/}  dt 

dR     cos  i  Ja(l  —  ez)da         aecosi     de       .     . 


.     . 
-j  —  sin  i*Ja(\.  —  e)  ji 

-<& 

dR       .    .    —  -  .  ^^ 


Whence,  by  an    easy   elimination,  and    observing    that 
7i2a3  =  /i,  we  deduce 


an  (1  -  e2)  dR 

tfyU.  C?^ 

2ofn  dR      a 


_      _ 
da  pe         de  ' 

dR  an  cos  ^ 


de 


p  sn          —  e 


sin  i  \/(l-e2) 
an  cos  i        (Z-R  an  dR 


fju  sn  -)  ^      /^  sin  i  V(l  -  e2)  ^^  ' 

formulae  from  which  we  can  readily  obtain  those  of  Art.  39  of 
this  work  when  we  take  into  consideration  that 

I  =  €  —  -ST,   I/  =  -57  —  O. 

For  by  Arts.  11  and  13,  R  may  be  expressed  as 
4>(a,e,i,r,  0-Q,,ty, 
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or  as  either  of  $  (a,  e,  i,  nt  +  e  -  «r,  w  -  H,  H)  =  E2, 


<f>(a,e,  i,nt+l,  v,  H)  = 
and  consequently 


dR 


_ 

dm  dl        dv  '      de      ~d 


and  so  dR,_dRz      dR,  _dR      dR. 

dl  ~~fa>     ~dp~~d^  +  ~de~' 

dR,     dR,     dR^     dR2 

d£l  ~  dii  +  dm       de  ' 
by  aid  of  which  the  needful  changes  may  be  effected. 

NOTES  TO  ARTS.   31   AND   33. 

7.  The  process  of  Arts.  81  and  33  may  perhaps  present 
some  difficulty  to  the  student,  arising  from  the  fact  that  the 
method  of  determining  the  position  of  the  planet  by  its  polar 
co-ordinates  on  the  plane  of  the  orbit,  i.e.  on  a  plane  moving 
with  the  planet,  is  not  applicable  to  the  geometrical  purpose 
of  denning  the  position  of  a  point  on  the  curve  of  reference, 
a  curve  which,  although  it  passes  through  the  planet  at  some 
given  instant,  has  nothing  to  do  with  the  planet's  subsequent 
motion. 

Of  course  the  position  of  any  point  may  be  determined 
by  its  polar  co-ordinates  on  a  plane  passing  through  it  and 
some  fixed  point  or  origin,  the  inclination  of  this  plane  to 
some  fixed  plane  through  the  origin,  and  the  longitude  of 
its  node.  But  inasmuch  as  an  infinite  number  of  planes 
can  be  drawn  through  two  given  points,  some  further  con- 
dition is  necessary  to  regulate  the  motion  of  that  on  which 
the  polar  co-ordinates  are  measured.  If  we  are  concerned 
with  the  motion  of  the  planet,  the  condition  is  that  this 
plane  shall  contain  the  direction  of  the  planet's  motion  at 
the  instant  under  consideration,  in  which  case  it  becomes 
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the  plane  of  the  orbit :  if  we  are  concerned  with  any  other 
point,  then  some  different  condition  is  requisite.  In  apply- 
ing this  method  to  the  curve  of  reference,  we  are,  of  course, 
at  liberty  to  assign  any  condition  we  please,  but  it  will  be 
convenient  that  it  be  such  that  the  co-ordinates  of  that  point 
of  the  curve  of  reference  which  coincides  with  the  planet 
shall  be  identical  with  those  of  the  planet  itself:  this  object 
is  attained  by  the  two  conditions  imposed  in  Arts.  31  and  33. 


ON   THE  EQUATIONS   OF  MOTION   OF  A  DISTURBED  PLANET. 

8.  The  principles  of  the  conservation  of  areas  and  of  the 
vis  viva  being  applicable  to  the  motion  of  the  Sun  and  any 
number  of  mutually  disturbing  planets,  four  first  integrals 
of  the  equations  of  motion  may  be  found.  Although  these 
integrals  are  of  little  use  in  determining  the  motions  of  the 
several  bodies,  which  cannot  be  completely  found  except  by 
methods  of  approximation,  we  shall  shew  how  they  may  be 
obtained  in  the  case  of  the  Sun  and  two  planets. 

The  differential  equations  of  motion  of  the  two  planets 
relative  to  the  Sun  are 


(1), 


M+ 
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where      .       p  =  {(at  -  xf  f  (y'  -  y)2  +  (zf 
From  these  we  find 


=  mm    -- 


and 


f\  '      *-  -1    \    /       f 

ro)mi»  (-3-p3)  (ya 


Hence 


my  +  mv 


Integrating,  we  have 
dz 


dz        ,  dz 


dy        ,  dy 


with  two  similar  integrals. 

These  equations  may  be  written  thus 

dz  ,  dy'\\ 
-rr  -  z  '-%  U- 
dt  dtj) 


.,  f     /    dz        dy\        ,  f  , 
M  4  m  (y  -T-  -  z  -£  }  +  m  [y 
{     V  dt         dtj          V 


with  two  similar  equations. 

For  a  first  approximation,  in  consequence  of  the  magni- 
tude of  the  Sun's  mass  by  comparison  with  that  of  any  of  the 
planets,  we  may  neglect  the  term  multiplied  by  mm',  and 
we  have 
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m>\y-j7  — 

dx        dz 


f   dy       dx 
m\x-~-  — -y-77 


Now 

is  the  projection  of  the  area 


dx        .  dz 


dx 


=  c, 


described  in  the  instantaneous  elliptic  orbit,  on  the  plane  of 
#y,  to  which  that  orbit  is  inclined  at  an  angle  i.     Hence 

(*  §  -  y  J)  =  J(M+  m)a(l-  e')  (1  +  tan2  i)'  *, 
and  we  have 


m          +  m)a(I-  (?)  (1  +  tan2  1)  " 

+  m  J(M+  m')  a  (1  -  e*)  (1  +  tanV)~*  =  C3, 
or,  retaining  only  terms  affected  by  M, 
m  Ja(I-e*)  (1  +  tan2  *)  ~* 

+  m'  Ja  (1  -  e'2)  (1  +  tan2  i')~*  =  constant. 

If  we  multiply  the  differential  equations  (1)  of  motion  of 
the  planet  m  by  the  factors 

dx        ,  dxr 
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and  multiply  also  the  differential  equations  (2)  of  motion  of 
the  planet  m  by  similar  factors,  and  add  the  six  products,  we 
shall  find 

_     fdx  6?x     dy  $y     dz  $?z 
\dt  dt*      dt  dtf      dt  o 

,(dz' 


O 

—  2 


dx        ,  dx 
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+  m+m 


(   tfx  ^    ,tfx'\ 

r  [*»-53  +  w^  -75- 
V     dt*  dfj 


dy        ,  dy' 

0™3T;     '  Tt 
-  2  -T  -  •  —  r        - 
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dr     m  dr 
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dy        ,du'\* 


=}=0 
-  ccY  +(y'-  2/)'  +  (z  -  z?> 

or   ntegrating 

f/cfoy     fdy\*  ,  /<M2)   ,      ,  f  /fl^N1     /cZyV  ,  fds/\*\ 

m|U)  +(J)  -|-U)i+m  1U)  +(i)  +  U); 
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an  equation  which  may  be  written 
ir    (fdx\*      (dy\* 

Mmdi)  +  uD  + 

dy      dy\z      fdz 


=  const., 

an    equation    containing    the    principle    of   the    vis    viva. 
Pontdcoulant,  Tome  L,  p.  210. 


EXAMPLE  OF  NUMERICAL  '  CALCULATION. 

9.  By  way  of  illustration,  we  will  here  shew  how  to 
calculate  the  inequality  in  the  earth's  radius  vector  due  to 
the  term 

m'<72  cos  2  {(n  -  ri)  t  +  e  -  e'{ 

in  R,  when  Jupiter  is  the  disturbing  planet. 

Following  the  method  of  Chapter  VIL,  and  taking  the 
result  of  the  first  approximation,  we  have  (Art.  110) 


or,  writing  for  Pl   its   value,  the   coefficient   of  this  term 
becomes 


m 


or  since  ri*as  =//,=  M  +m,  if  M  denote  the  Sun's  mass,  we 
may  write,  to  the  first  order  of  the  disturbing  force, 


*C,\ 

M  (3»  -  2ri)  (n  -  2n')        da       n  -  n  j  ' 
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a  more  convenient  expression  ;  for  since  n  occurs  in  the  form 
of  a  ratio,  it  is  no  longer  necessary  to  suppose  it  expressed  in 
circular  measure. 

The  values  of  C2  and  -^-2  must,  of  course,  be  found  by 

the  methods  of  Chapter  II.  :  by  slightly  changing  our  nota- 
tion, however,  we  shall  be  able  to  avail  ourselves  of  the 
calculation  of  M.  Bouvard,  as  given  in  the  third  volume  of 
the  Mecanique  Celeste,  and  the  third  volume  of  Ponte- 
coulant's  Theorie  Analytique  du  Systtme  du  Monde. 

If  c2  =  a  Cz,  a.  =  —  ,  we  have 
a 


2     a"    da      ada~~a'*dx  ' 
thus  the  coefficient  becomes 


M  (3n 

or,  taking  the  Sun's  mass  and  the  Earth's  mean  distance  as 
the  units  of  mass  and  of  distance, 


2 

^  ~ 


We  take  from  the  Systems  du  Monde,  Tome  in.*, 


logw  a  =1-2838993,  logw  cz  =  2-4501721, 

dr 

loglo  a  ^2  =  27581084; 

*  In  the  second  edition  of  the  first  volume,  published  in  1856,  more 
recent  values  of  m',  n  and  n'  are  given,  but  our  object  being  simply  to  illus- 
trate the  method  of  calculation,  the  above  values  have  been  retained. 

11—2 
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and  with  the  aid  of  a  table  of  logarithms  the  numerical  value 
of  the  coefficient  may  easily  be  found.  Thus,  in  consequence 
of  the  term  in  R  we  have  been  considering,  we  find  in  Sr  the 
term 

-  '0000092122  cos  2  {(n  -  ri)  t  +  e  -  e'}. 

ON  THE  METHODS  OF  CALCULATING  THE  MASSES  OF 
THE  PLANETS. 

10.  There  are  in  general  two  methods  of  determining  the 
masses  of  the  planets ;  either  by  observations  on  a  satellite, 
when  the  planet  is  accompanied  by  a  satellite ;  or  by  compar- 
ing the  inequalities  produced  in  their  motion  by  their  mutual 
action,  as   deduced   from   observation,  with   the  same  ine- 
qualities calculated  from  theory.     The  secular  variations  are 
best  adapted  to  give  the  most  exact  results ;  but  these  are  not 
yet  known  with  sufficient  accuracy  to  allow  of  this  use.     We 
are  therefore  obliged  to  recur  to  the  periodic  variations,  and, 
by  combining  a  vast  number  of  observations,  gather  from 
them  the  most  probable  results  *. 

11.  When  the  planet  is  accompanied  by  a  satellite  the 
formula  for  calculating  its  mass  may  be  obtained  as  follows  : 

Let  M,  ra,  m'  be  the  masses  of  the  Sun,  the  planet,  and 
the  satellite :  P,  P'  the  periodic  times  of  the  planet  about 
the  Sun,  and  the  satellite  about  the  planet ;  a,  a  the  mean 
distances  of  the  planet  from  the  Sun,  and  the  satellite  from 
the  planet.  Then  we  have 

2?raf  27ra'f 

JL     ==    ~ ~ ~~  ,  JL        ==  7~    I 

ro+ra'     PV3 
therefore  y^T  =  FV  ' 

*  Pont6coulant'   Systeme  du  Monde,  Tome  m.  p.  340. 
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or  approximately, 

m+m'_PV3 
M     "  PV  ' 

This  equation  gives  the  mass  of  the  planet  when  that  of 
its  satellite  is  known.  If  the  latter  be  neglected,  the  formula 
becomes 

m      PV3 


12.  In  the  case  of  the  Earth,  this  method  is  not  suf- 
ficiently exact,  but  the  following  may  be  employed.  The 
attraction  of  the  Earth  on  a  body  at  its  surface,  in  the 

parallel  of  which  the  square  of  the  sine  of  the  latitude  is  ~  , 

is  very  nearly  the  same  as  if  the  Earth  were  condensed  into 
its  centre.  (See  Pratt's  Figure  of  the  Earth,  Art.  89.)  Let 

then  sin2Z  =  ^,  g  =  ihe  Earth's  attraction  on  a  body  at  its 

o 

surface  in  latitude  I,  b  the  mean  radius  of  the  Earth,  E  the 
mass  of  the  Earth,  M  the  mass  of  the  Sun,  P  the  length  of 
the  year,  and  a  the  mean  radius  of  the  Earth's  orbit.  Then 

E 


E     gVP*      gPz  fb 
therefore  -^  =  7-5-=  =    -sr   - 

M     4?rV 


where  -  =  sine  of  Sun's  parallax. 

a 

For  the  methods  of  calculating  the  mass  of  the  Moon 
we  refer  to  Pont^coulant's  Systeme  du  Monde,  Tome  IV. 
p.  651. 


166  PLANETARY   THEORY. 


ON  THE  CONSTRUCTION   OF  ASTRONOMICAL  TABLES. 

13.  In  the  present  state  of  Astronomy,  the  values  of 
the  masses  and  elements,  upon  which  these  tables  depend, 
may  be  considered  as  approximately  known,  so  that  small 
corrections  only  are  necessary  :  we  propose,  in  this  note,  to 
shew  briefly  how  these  corrections  are  effected. 

A  series  of  observations  are  taken  in  KA.  and  N.P.D., 
separated  by  considerable  intervals  of  time  :  from  each  of 
these  the  geocentric  longitude  and  latitude  are  obtained,  and 
thence  the  heliocentric  longitude  and  latitude.  These  co- 
ordinates are  also  calculated  for  each  period  of  observation 
by  the  methods  of  the  preceding  Chapters,  employing  the 
existing  values  of  the  masses  and  elements.  The  differences 
between  the  observed  and  calculated  values  will  be  due  to 
errors  in  the  latter,  for  we  assume  .  that  our  observations  can 
be  depended  upon.  From  these  errors  in  longitude  and 
latitude,  the  corresponding  corrections  in  the  values  of  the 
masses  and  elements  are  obtained  as  follows  : 

Let  6  denote  the  longitude  of  the  planet  as  calculated 
from  theory,  6  +  $0  its  value  as  given  by  observation,  m,  m", 
&c.,  the  existing  values  of  the  masses  of  the  disturbing 
planets,  $m',  Sm",  &c.,  &n,  Se,  &sr,  8e  the  requisite  corrections 
to  the  values  of  the  masses  and  elements.  Then,  retaining 
only  the  first  power  of  the  excentricity  and  of  the  disturbing 
force,  we  have  by  the  methods  of  the  preceding  Chapters, 


where  m'P',  m"P",  &c.,  are  the  terms  due  to  the  perturba 
tions  of  the  planets  m',  m",  &c.     It  follows  that 


~ 
dn          de          cfa-  de 
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where  80,  being  the  difference  between  the  observed  and 
calculated  values  of  the  longitude,  is  known. 

A  similar  equation  will  result  from  each  observation  of 
the  series,  and  the  number  of  these  must  be  reduced  by  the 
method  of  least  squares  to  that  of  the  unknown  quantities: 
thus,  by  the  solution  of  these  equations,  the  required  cor- 
rections to  m',  m",  &c.,  n,  e>  OT,  e  will  be  obtained. 

In  like  manner,  from  a  comparison  between  the  observed 
and  calculated  values  of  the  latitude,  the  corrections  to  i  and 
fl  may  be  obtained. 

Tables  of  the  numerical  values  of  the  masses  and  elements 
will  be  found  in  Herschel's  Outlines  of  Astronomy,  and  in 
the  first  volume  of  Pontecoulant's  Systtme  du  Monde. 


THE    END. 
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